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1.  INTRODUCTION 


Recorded  Simulation  Control  Scheme  is  a method  of  generating  input 
signals  for  vehicle  shake  tests  performed  in  laboratory.  In  this  scheme  a 
computer  program  is  used  to  simulate  operation  of  the  vehicle  being  tested 
over  a terrain  of  interest.  The  simulation  results  create  axle  displacement 
records  which  are  stored  in  memory  of  a control  system.  At  the  time  of 
shake  test,  these  records  are  played  back  to  generate  axle  displacement  sig- 
nals. Thus  the  scheme  involves  two  distinct  steps  of  operation;  (a)  creating 
axle  displacement  records,  and  (b)  generating  shaker  input  signal,  as  shown 
in  Figure  1 . 

There  are  several  ways  of  implementing  the  recorded  simulation  con- 
trol scheme,  with  the  same  basic  approach  of  recording  the  computer  gener- 
ated axle  displacement  signals  and  playing  them  back  at  the  time  of  shake 
test.  The  implementation  scheme  selected  in  this  effort  is  shown  in  Figure  2, 
As  shown  in  the  figure,  the  axle  motion  records  are  generated  by  dynamic 
vehicle  simulation  computer  program.  Data  necessary  to  execute  the  pro- 
gram include: 

a.  Vehicle  data  (mass,  CG  location,  wheel  base,  etc.) 

b.  Operation  data  (terrain  profile  and  speed  of  operation). 

The  axle  displacement  records  generated  by  the  simulation  are  stored 
on  paper  tapes  and  transferred  to  magnetic  disc  of  a digital  controller  through 
adequate  software.  The  data  stored  on  the  magnetic  disc  are  retrieved  at  the 
time  of  running  the  shake  test,  and  the  digital  signals  are  converted  to  analog 
signals  using  a bank  of  digital  to  analog  signal  converters.  The  analog  sig- 
nals are  then  fed  to  the  shaker  through  appropriate  filters  and  amplifiers. 

The  records  stored  on  the  discs  typically  last  less  than  a minute  in 
real-time;  therefore,  the  same  records  are  repeated  until  the  testing  is 
completed.  Also,  the  simulation  program  is  executed  each  time  the  vehicle 
test  parameters  (terrain  profile,  speed,  load)  are  changed  and  the  axle  dis- 
placement records  on  the  discs  are  updated  to  show  the  changes. 
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Step  1;  Creating  Axle  Displacement  Records 


Step  2;  Generating  Shaker  Input  Signals 


Figure  1,  Recorded  Simulation  Control  Scheme 
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Figure  2,  Implementation  of  Recorded  Simulation  Control  Scheme 


This  manual  describes  both  steps  of  the  recorded  simulation  control 
scheme,,  Section  2 deals  with  creating  axle  displacement  records  from  the 
specifications  of  vehicle  and  terrain,  whereas  Section  3 discusses  the  con- 
trol system  which  converts  the  records  to  shaker  input  signals.  Instructions 
on  how  to  execute  both  the  steps  are  given  along  with  appropriate  illustrative 
examples. 


4 


2.  STEPl:  AXLE  DISPLACEMENT  RECORDS  FROM 
VEHICLE  AND  TERRAIN  SPECIFICATIONS 


This  chapter  deals  with  the  first  of  two  steps  involved  in  implementing 
the  recorded  simulation  control  scheme.  The  details  of  the  first  step  are 
illustrated  by  Figure  3.  As  shown  in  the  figure,  the  input  to  the  first  step  is 
in  form  of  specifications  for  the  vehicle  and  operation,  where  the  vehicle 
specifications  include  parameters  such  as  dimensions  of  the  vehicle,  spring 
constants  of  the  suspensions,  and  weights  of  the  components.  The  operation 
specifications,  on  the  other  hand,  include  details  of  the  terrain  profile  and 
speed  of  the  vehicle.  If  the  profile  of  the  terrain  is  not  known,  its  roughness 
parameter  is  sufficient  to  generate  a random  profile  using  an  associated 
program,  GRND,  described  in  Appendix  A.  3. 

The  vehicle  data  and  terrain  profile  are  fed  to  the  dynamic  vehicle 
simulation  program  to  obtain  simulation  results.  The  output  of  the  program 
is  in  many  forms,  as  shown  in  the  figure.  The  printout  summarizes  the 
input  parameters  and  the  dynamic  simulation  results.  The  simulation  results 
are  in  form  of  values  of  the  key  vehicle  variables  listed  as  a function  of  time. 
The  program  also  can  generate  plots  of  the  variables,  either  versus  time  or 
versus  horizontal  vehicle  displacement.  The  axle  displacement  results 
generated  by  the  simulation  program  are  punched  on  cards  and  they  can  sub- 
sequently be  transferred  to  paper  tapes  which  store  data  in  standard  ASCII 
format.  As  an  option  available  to  the  user,  the  output  of  some  of  the  vari- 
ables can  be  stored  on  a tape  for  future  statistical  processing  by  an  associated 
program,  SPEC,  described  in  Appendix  A.  3. 

The  rest  of  this  section  deals  with  the  dynamic  vehicle  simulation 
program. 

The  operation  of  the  program  is  described  in  Section  2.1.  Section  2.2 
deals  with  instructions  on  how  to  use  the  program  and  Section  2.3  illustrates 
the  use  with  an  example.  Details  of  the  subroutines,  program  nomenclature, 
associated  programs  (CRND  and  SPEC)  and  listing  of  the  program  are  in 
Appendix  A. 
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Vehicle 

Simulation 


Dynamic  Vehicle  Simulation  Program 


The  terrain-tire -vehicle  model  developed  in  reference  1 is  incorporated 
in  a computer  program  with  specifications  summarized  in  Table  1.  The 
program  simulates  operation  of  a vehicle,  shown  schematically  in  Figure  4, 
over  a given  terrain  at  a given  speed.  In  the  diagram  describing  the  vehicle 
model,  the  tire  models  are  not  shown  because  the  program  offers  a choice 
of  four  tire  models.  These  models  are  briefly  described  in  Table  2.  For 
details  of  the  vehicle  and  the  tire  models,  the  user  is  advised  to  read  ref- 
erence 1 . 


Table  1.  Specifications  of  the  Computer  Program 


Language: 

1 

FORTRAN  IV 

Compiler: 

Run  (Control  Data) 

Core  Used: 

70K  Octal 

Additional  Requirements; 

CALCOMP  Subroutines 

Plotter 

Card  Punch 

Paper  Tape  Punch-r 

For  a typical  vehicle  simulation  over  500  feet  of  terrain: 


! 

Using  Any  of  The 
First  Three  Tire 
Models 

Using 

Adaptive  Footprint 

Tire  Model 

Typical 

1 

1 

CPU  time 
on  CDC  6400 

100  sec 

1762  sec 

Typical  cost 
(government  rate) 

$18 

$58 

^Required  if  the  control  system  needs  the  axle  displacement  records  on 
paper  tapes. 

(a)  M-809  Truck 


(c)  The  Model  of  the  Vehicle 


Figure  4.  Formulation  of  the  Vehicle  Model 


The  computer  program  incorporates  the  terrain-tire -vehicle  model  to 
simulate  operation  of  the  vehicle  through  the  following  five  steps: 

a.  Vehicle  parameter  input 

b.  Profile  data  input 

c.  Initial  condition  calculation 

d.  Dynamic  simulation 

e.  Plotting  variables. 

These  steps  are  executed  by  the  program  through  a number  of  sub- 
routines. The  usage  of  subroutines,  each  of  which  performs  a specific  func- 
tion, makes  the  program  modular  in  structure  so  that  it  can  be  modified  just 
by  changing  the  relevant  subroutines  and  keeping  the  rest  of  the  program 
untouched.  The  manner  in  which  the  main  program  interacts  with  the  sub- 
routines is  shown  in  Figure  5.  How  this  interaction  performs  the  above 
five  steps  is  described  in  the  following; 

a.  Vehicle  Parameter  Input 

Execution  of  the  program  begins  by  calling  subroutine  PROGIO 
which  reads  the  input  data  for  the  vehicle  and  the  tire  models  through  three 
subroutines,  lOLAB,  lOLABD  and  IOLAB2.  The  data  is  in  form  of  data 
cards  with  an  alphanumeric  description  of  the  parameter,  its  value  and  the 
unit.  This  description  for  each  data  card  is  printed  as  soon  as  it  is  read. 
The  parameters  which  are  not  in  consistent  units  are  then  converted  in  the 
main  program  just  after  calling  PROGIO. 

b.  Profile  Data  Input 


Profile  data  consist  of  elevations  of  a random  profile,  with  an  rms 
value  of  1 foot,  generated  by  one  of  the  associated  programs,  GRND,  which 
is  described  in  Appendix  A.  This  profile  data  are  read  by  subroutine 
PROFILE.  Because  the  program  may  simulate  vehicle  operation  on  an  end- 
less profile,  all  the  profile  data  are  not  read  at  once.  Initially  only  as  much 
profile  is  read  as  the  program  dimensions  can  absorb.  Then  the  simulation 
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Subroutines 


Main  Program 


Figure  5.  Main  Program  - Subroutine  Interaction 
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Subroutine  s 


Main  Program 


\ 


POINTC 
(ITIRE=1, 2,3) 

ADAPT  - ADORG 
{ITIRE=4) 

SUSP 

BOGIE 


STEQU 

ADAPT 

(ITIRE=4) 

■STEQU  - RKDIF 


1 


\ 


IDYSYS 


/ 


\ 


RTREAD 

(ITIRE=2) 

FPRINT 

(ITIRE=3) 


hPROFILE 


/ 


^ CALCOMP  Subroutines 
**  S'/ stem  Subroutine 

(Not  in  the  Source  Deck^ 


PSTORE 

SECOND 

INITPLT 

>!t 

PLOT 


\ 


Scale 
* 

Axis 

Number 

* 

Line 

* 

Plot 


HVPLOTS 


/ 


FIN 


/ 


Step  d 


Step  e 


END 


Figure  5.  Main  Program  - Subroutine  Interaction 
Continued 
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is  carried  out  over  that  part  of  the  profile.  At  the  end  of  that  profile  part, 
another  profile  length  is  read  in  by  PROFILE,  and  the  process  is  repeated. 

The  profile  thus  read  in  is  filtered  by  subroutine  RTREAD  if  the 
rigid  tread  band  tire  model  is  being  used  and  by  subroutine  FPRINT  if  the 
fixed  footprint  model  is  being  used.  The  profile  is  then  converted  to  a pro- 
file with  the  user  specified  rms  value  (the  input  profile  has  an  rms  value  of 
1 foot). 


c.  Initial  Condition  Calculations 


From  the  data  of  the  terrain  on  which  the  vehicle  initially  rests, 
the  initial  conditions  of  all  the  state  variables  are  calculated  in  the  main 
program. 

d.  Dynamic  Simulation 

Integration  of  the  state  equations  to  determine  time  histories  of 
the  state  variables  is  coordinated  by  subroutine  DYSYS.  The  subroutine  calls 
subroutine  STEQU  in  the  beginning  to  initialize  values  of  the  state  variable 
derivatives.  Subroutine  ADAPT  is  called  next,  if  the  adaptive  tire  model  is 
being  used,  to  determine  the  initial  elevation  of  the  tire  center.  After  these 
initial  calls,  subroutine  DYSYS  executes  the  integration  process  by  calling 
subroutine  RKDIF  at  each  time  steps.  Subroutine  RKDIF  incorporates  a 
fourth  order  Runge-Kutta  algorithm  to  numerically  integrate  the  state  equa- 
tions. Through  integration  it  updates  the  values  of  the  state  variables  (dis- 
placements, velocities  of  the  vehicle  component)  at  each  time  step. 

The  state  equations  are  programmed  in  subroutine  STEQU,  which 
is  called  four  times  at  each  time  step  by  subroutine  RKDIF.  These  state 
equations  are  differential  equations  describing  the  heave -pitch  motions  of  the 
vehicle  and  its  components.  In  order  to  determine  updated  values  of  the 
derivatives,  subroutine  STEQU  calls  subroutine  POINTC,  if  any  of  the  first 
three  tire  models  (point  contact,  rigid  tread  band  or  fixed  footprint)  are 
being  used.  In  case  of  the  adaptive  footprint  tire  model,  however,  subrou- 
tine ADORG  (and  through  ADORG,  subroutine  ADAPT)  is  called  instead  of 
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POINTC.  These  subroutines  calculate  the  forces  exerted  by  the  tires. 

Forces  exerted  by  the  suspension  are  determined  by  calling  subroutine  SUSP, 
whereas  those  exerted  by  the  bogie  are  determined  from  subroutine  BOGIE. 
From  updated  values  of  forces  thus  obtained,  the  values  of  derivatives  are 
calculated  by  subroutine  STEQU,  from  which,  in  turn,  the  updated  values  of 
the  state  variables  are  calculated  by  subroutine  RKDIF. 

Subroutine  DYSYS  generates  outputs  from  the  results  of  numerical 
integration  in  the  following  manner 

(1)  The  values  of  the  key  vehicle  variables  are  printed  after  MM 
(user  selected)  time  steps. 

(2)  The  values  of  time;  and  front,  middle  and  rear  axle  displace- 
ments; are  punched  at  each  time  step.  The  punched  records 
can  be  converted  to  those  on  paper  tapes  according  to  the 
equipment  availability  at  the  user's  installation. 

(3)  If  plotting  of  certain  variables  is  desired,  subroutine  PSTORE 
is  called  to  store  the  variables  to  be  plotted.  The  plots  are 
then  produced  by  subroutine  VPLOTS. 

(4)  If  statistical  processing  option  is  to  be  exercised,  additional 
cards  are  required  by  subroutine  DYSYS  (as  described  in 
Appendix  A.  3)  for  storing  the  variables  on  magnetic  tape. 

Numerical  integration  is  halted  when  additional  profile  data  is  to 
be  read.  Subroutine  PROFILE  is  called  and  with  the  new  profile  data  the 
the  integration  is  resumed.  The  process  is  repeated  until  the  end  of  the 
profile  is  reached.  The  execution  is  then  moved  to  the  next  step.  During 
the  simulation,  subroutine  SECOND  keeps  track  of  the  execution  time. 
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e. 


Plotting  Variables 


If  the  user  requires  some  of  the  variables  to  be  plotted,  CALCOMP 
subroutines  INITPLT  and  PLOT  are  called  to  initialize  the  plotter.  Subroutine 
VPLOTS  is  then  called  to  create  6 by  5 inch  plots  of  the  variables  versus  time 
(or  displacement  as  requested  by  the  user).  The  plots  are  created  using 
CALCOMP  subroutines  SCALE,  AXIS,  NUMBER,  LINE  and  PLOT,  which  are 
called  by  subroutine  VPLOTS.  Plotting  is  termined  by  subroutine  DYSYS 
through  CALCOMP  subroutine  FIN. 

After  producing  the  plots,  the  program  is  terminated. 

2.2  User  Instructions 


2,2.1  Program  Input 

The  program  obtains  input  through  data  cards.  There  are  two 
main  categories  in  which  the  data  cards  can  be  divided: 

(1)  Vehicle -tire  parameters 

(2)  Terrain  data. 

The  vehicle-tire  parameters  category  includes  masses  of  the 
vehicle  components,  stiffness  of  the  suspensions,  tire  characteristics,  etc. 

The  terrain  data  is  in  form  of  ground  elevation  readings  at  a fixed  horizontal 
distance.  There  is  no  limit  to  the  length  of  the  profile  over  which  the  simu- 
lation can  be  carried  out.  Also,  if  profile  elevation  data  are  not  available, 
associated  program  GRND  (described  in  Appendix  A.  3)  can  be  executed  to 
generate  such  data. 

In  the  following  pages  the  data  cards  necessary  to  execute  the 
program  are  listed.  The  numerical  values  in  these  data  cards  are  to  be 
placed  within  the  'F',  'I'  or  'E'  format.  The  'A'  fields  are  used  to  provide 
parameter  name  and  units,  and  other  descriptive  legends  for  user  convenience. 
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DATA  CARDS 


General  Description  Data 

Card  No. 

Contents 

Format 

1-6 

Header  cards  (the  user  can  print  a 
heading  using  these  cards) 

80A1 

7 

Vehicle  type 

80A1 

8 

Tire  type 

80A1 

9 

Terrain  No. 

80A1 

10 

Tire  model  no.,  ITIRE 

ITIRE  = 1 Point  contact  model 

= 2 Rigid  tread  band  model 
= 3 Fixed  footprint  model 
— 4 Adaptive  footprint  model 

40A1,I5 

11 

Vehicle  speed,  V (mph) 

30A1,F10'4,  lOAl 

12 

Profile  rms,  y (in. ) 

(see  Note  1 ) 

30A1,F10«4,10A1 

Vehicle 

-Tire  Parameters  (See  Figure  4) 

Vehicle  Frame  Parameters 

13 

Header  card 

80A1 

14 

Mass  of  hull,  Mj^  (slugs) 

30A1,F10*4,10A1 

15 

Mass  of  bogie,  M^  (slugs) 

30A1,F10*4,10A1 

16 

Front  unsprung  mass,  each, 

M2  (slugs) 

30Al,Fl0-4,  lOAl 

17 

Middle  unsprung  mass,  each, 

M3  (slugs) 

30A1,F10*4,  lOAl 

18 

Rear  unsprung  mass,  each, 

M4  (slugs) 

30Al,Fl0-4,  lOAl 

19 

2 

Pitch  inertia  of  hull,  (slug  ft  ) 

30A1,F10'4,10A1 

20 

Pitch  inertia  of  bogie,  (slug  ft^) 

30Al,Fl0-4,  lOAl 

16 


21 

Horizontal  distance  between  front 
suspension  and  hull  CG,  L2  (ft) 

30Al,Fl0-4,  lOAl 

22 

Horizontal  distance  between  bogie 
connection  and  hull  CG,  (ft) 

30A1,  FlO-4,  lOAl 

23 

Height  of  hull  CG  from  hull  plane, 

GG  (ft) 

30A1,F10*4,  lOAl 

24 

Vertical  distance  between  bogie 
connection  and  hull  plane,  S34  (ft) 

30A1,F10*4,  lOAl 

25 

Horizontal  distance  between  bogie 
connection  and  front  stop,  (ft) 

30Al,Fl0-4,  lOAl 

26 

Horizontal  distance  between  bogie 
connection  and  rear  stop,  (ft) 

30Al,F10-4,  lOAl 

27 

Vertical  distance  between  bogie 
connection  and  bogie  CG,  (ft) 

30A1,F10*4,10A1 

28 

Half  bogie  length,  a (ft) 

30A1,F10*4,  lOAl 

29 

Stop  size,  h (ft) 
s 

30A1,F10*4,10A1 

30 

Stiffness  factor,  bogie  - stop, 

K , (see  Note  2) 
stop ' ' 

30A1,F10*4,  lOAl 

31 

Stiffness  factor,  bogie  - hull, 

H (see  Note  2) 

stop  ' ' 

30A1,F10‘4,10A1 

32 

Unloaded  clearance  between  frame 
stops  and  bogie,  ^^ga,p 

30A1,F10*4,  lOAl 

Suspension  Parameters  (See  Note  3) 

33-34 

Header  cards 

80A1 

Front  Suspension 

35 

Stiffness,  K2  (lb /ft) 

30A1,F10-4,10A1 

36 

Viscous  damping  constant  (jounce), 

(lb«  sec/ft)  (See  Figure  6) 

30Al,Fl0-4,  lOAl 

37 

Viscous  damping  constant  (rebound), 
B2vr  (See  Figure  6) 

30Al,Fl0-4,  lOAl 

38 

Jo\ince  travel,  Jn^  (ft) 

30A1,F10*4,  lOAl 

39 

Rebound  travel,  Rb2  (ft) 

30A1,  FlO-4,  lOAl 

17 


Forrfi  ' Force 


(a)  Dry  Friction  Damping  (b)  Shock  Absorber  Damping 


Figure  6.  Damping  Models 
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40 


Middle  and  Rear  Suspensions 
Header  card 

41  Stiffness,  (Ib/ft) 

42  Jounce  travel,  Jn^  (ft) 

43  Rebound  travel,  Rb^  (ft) 

44  Stop  stiffness  factor,  n (ft) 

(See  Note  2) 

Tire  Parameters 

45  Header  card 

46  Tire  stiffness,  K (Ib/ft) 

(See  Note  4) 

47  Tire  damping  constant,  b 
(lb*  sec/ft)  (See  Note  5) 

48  Tire  radius,  r (ft) 

49  Tire  footprint  length,  L (ft) 

(See  Note  4) 

50  Tire  distributed  (cubic)  stiffness, 
k”  (Ib/in.  ^)  (See  Note  6) 

51  Tire  distributed  (cubic)  damping, 
b"  (Ib’sec/ft^)  (See  Note  6) 

52  Tire  inflation  pressure,  p.  (psi) 

53  Tire  deflection  limit,  T^  (ft) 

54  Tire  stiffness  factor,  T 

(See  Note  2)  ® ^ 

Simulation  Data 


80A1 

30A1,F10*4,  lOAl 
30A1,  FlO-4,  lOAl 
30Al,Fl0-4,  lOAl 
30Al,Fl0-4,  lOAl 


80A1 

30A1,  F10*4,  lOAl 

30A1,F10*4,  lOAl 

30A1,F10*4,10A1 
30Al,Fl0-4,  lOAl 

30A1,F10*4,10A1 

30A1,F10‘4,10A1 

30Al,Fl0-4,  lOAl 
30A1,F10'4,  lOAl 
30A1,F10  4,  lOAl 


Printout  Card 

55  MM,  where  the  dynamic  results  are  40A1,I5 

printed  at  every  MM  time  steps 
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Plot  Card 


56 


{JPLT(K),  K=l,  40),IXAX,PLSTR,PLSTP  41Il,  9X,2F10*  3 


Column  Explanation 

1-40  Punch  1 in  the  Kth  column  if 

Kth  variable  is  to  be  plotted 
(See  Table  3).  Eight 
variables  maximum. 


41 


51-60 


61-70 


IXAX=0,  If  X axis  on  plots  is 
distance  travelled  (XREAR) 
IXAX=1,  If  X axis  is  time 

Distance  travelled  at  the  first 
point  of  the  plot  (as  measured 
by  rear  wheel  center  horizontal 
displacement) 

Distance  travelled  at  the  last 
point  of  plot  (as  measured  by 
rear  wheel  center  horizontal 
displacement) 


Terrain  Data 


57  (1)  Fixed  horizontal  distance  between  the  FlO'2,15 

data  points,  dx  (ft) 

(2)  Total  number  of  data  points,  JMAX 


58  Terrain  elevation  data,  in  ten  columns,  8E10*3 

and  eight  on  each  card  (ft), 

onwards  TER(  J),  J=1 , JMAX 


Note  1;  This  assumes  that  the  terrain  profile  used  in  the  illustrative 

simulation  has  been  generated  by  an  associated  program,  pro- 
gram GRND  described  in  Appendix  A.  3.  This  profile  has  an 
rms  (root  means  square)  elevation  value  of  1 foot.  A profile 
with  a different  rms  value,  but  with  a similar  spectral  density 
can  be  obtained  by  simply  multiplying  the  generated  profile  by 
the  rms  value  supplied  by  the  user  through  this  data  card.  If, 
however,  a terrain  profile  with  correct  rms  value  is  already 
available  to  the  user,  either  through  experimental  results  or 
through  some  other  terrain  generating  program,  this  data  card 
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Table  3.  List  of  Variables  Available  for  Plotting 


Variable 

Number 

Variable 

Program 

Symbol 

Dimension 

1 

Hull  pitch  angle 

THETA 

degree 

2 

Hull  pitch  velocity 

DTHETA 

degree/sec 

3 

Bogie  pitch  angle 

PHI 

degree 

4 

Bogie  pitch  velocity 

DPHI 

degree /sec 

5 

Hull  CG  displacement 

YO 

ft 

6 

Hull  CG  heave  velocity 

DYO 

ft/ sec 

7 

Front  axle  displacement 

YW2 

ft 

8 

Front  axle  heave  velocity 

DYW2 

f t / s e c 

9 

Middle  axle  displacement 

YW3 

ft 

10 

Middle  axle  heave  velocity 

DYW3 

f t / s e c 

11 

Rear  axle  displacement 

YW4 

ft 

12 

Rear  axle  heave  velocity 

DYW4 

f t / s e c 

13 

Front  tire  force  (vertical) 
(expressed  as  fraction  of  the 
equilibrium  force) 

FV2/WTR2 

Dime  n s io  nle  s s 

14 

Middle  tire  force  (vertical) 

FV3/2WTR3 

Dimensionless 

15 

Rear  tire  force  (vertical) 

FV4/2WTR4 

Dimensionless 

16 

Front  tire  force  (horizontal) 

FH2/WTR2 

Dime  nsionless 

17 

Middle  tire  force  (horizontal) 

FH3/2WTR3 

Dime  nsionless 

18 

Rear  tire  force  (horizontal) 

FH4/2WTR4 

Dimensionless 

19 

Front  suspension  force  (vertical) 

FY2 

lb 

20 

Middle  suspension  force  (vertical) 

FSY3 

lb 

21 

Rear  suspension  force  (vertical) 

FSY4 

lb 

22 

Front  ground  elevation  (effective) 

Y2 

ft 

23 

Middle  ground  elevation  (effective) 

Y3 

ft 

24 

Rear  ground  elevation  (effective) 

Y4 

ft 

25 

Bogie  CG  elevation 

YB 

ft 

26 

Hull  CG  acceleration 

D2Y0 

ft/sec^ 

27 

Front  suspension  displn. 

DEL2 

ft 

28 

Middle  suspension  displn. 

DEL3 

ft 

29 

Rear  suspension  displn. 

DEL4 

ft 

30 

Original  terrain  profile 
(as  seen  by  the  rear  tires) 

OTER(KPHS) 

ft 

31 

Effective  (filtered)  profile 
(as  seen  by  the  rear  tires) 

TER(KPHS) 

ft 

should  have  a value  of  12.0  inch,  so  that  the  terrain  rms  value 
is  not  readjusted  in  the  program. 


Note  2: 


Note  3: 


Note  4: 


Stiffness  factor  is  defined  as  a factor  by  which  the  normal 

stiffness  should  be  multiplied  to  simulate  a hard  contact.  Thus, 

when  the  suspension  exceeds  its  jounce  or  rebound  limits,  its 

stiffness  gets  effectively  multiplied  by  a factor  n,  as  shown  in 

Figure  7.  Similarly  when  the  tire  deflection  exceeds  the  limit, 

Tjg,  its  stiffness  gets  multiplied  by  factor  (see  Figure  7 (cont)). 

For  the  frame  stops,  the  situation  is  different,  because  unless 

the  gap  between  the  stop  and  the  bogie  becomes  zero,  there  is 

no  force  acting  on  the  bogie  and  hence  the  stiffness  factor,  when 

the  contact  does  occur,  has  to  be  keyed  to  some  other  stiffness. 

This  stiffness  chosen  is  the  suspension  stiffness  for  the  middle 

and  the  rear  suspensions.  Thus,  the  stiffness  of  the  stop  is 

Ko  K , as  shown  in  Figure  7.  In  extreme  case,  when  the 
J s x>o  o 

bogie  touches  the  hull,  the  effective  stiffness  is 


One  of  the  missing  parameters  in  the  suspensions  data  is  dry 
friction  coefficient,  which  is  extremely  difficult  to  model. 
Initial  experiments  at  TARADCOM  have  indicated  that  a force 
value  equal  to  5 percent  of  the  instantaneous  spring  force  which 
opposes  the  velocity  describes  the  dry  friction  approximately. 


Typical  design  curves  supplied  by  the  tire  manufacturer  are 
shown  in  Figure  8.  From  such  a curve,  knowing  the  equili- 
brium load,  W,  and  inflation  pressure,  p^  (PRESS),  equilibrium 
deflection,  YST,  can  be  determined.  Inverse  of  the  slope  of  the 
curve  at  point  A will  give  the  stiffness  value  k.  The  following 
equation  determines  footprint  length,  L. 

. . / -1  /r-YST\\ 

L = 2r  sin  I cos  ^ J j 

where  r is  radius  of  the  tire.  The  user  should  refer  to  reference 
2 for  details  on  tire  stiffness. 
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Figure  8.  Typical  Design  Curves  for  Tire 
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Note  5: 


Note  6: 


Data  on  tire  damping  are  relatively  scarce.  Typical  experimen- 
tal results  on  tire  damping  have  indicated  the  damping  ratio,  ^ , 
to  be  in  a 0.03-0.08  range.  For  the  illustrative  simulations  in 
this  report  a mean  value  of  ^ = 0.05  has  been  assumed.  From 
I,  the  damping  constant  b can  be  determined  from  the  following 
relation 


b =2 


I 


where  m is  the  mass  of  the  tire  in  slugs  and  k is  the  stiffness 
in  Ib/ft.  More  details  on  the  tire  damping  are  available  in 
reference  2 . 

The  value  of  the  distributed  (cubic)  stiffness  of  tire,  which  is 
used  in  the  adaptive  footprint  model,  k",  can  be  evaluated  by 
considering  deflection  at  the  equilibrium  condition; 


w - p.  A 


where  B = effective  width  of  the  tire  =*  A/L 

A = effective  footprint  area  at  equilibrium  obtained 
from  Figure 

e = cos-V---™: 

w :=  equilibrium  force  on  the  tire. 

The  cubic  damping  coefficient,  b",  is  evaluated  from  the  tire 
damping  coefficient,  b: 


2Br  e 
o 
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2.2.2  Program  Output 


The  program  generates  outputs  in  three  forms; 

a.  Printout 

b.  Plotter  output 

c.  Punch  output 

The  fourth  form  of  output,  i.  e.  , magnetic  tape  for  statistical 
processing  of  the  stored  variables,  is  optional  and  the  user  is  referred  to 
Appendix  A.  3 for  more  details. 

2. 2. 2.1  Printout 

a.  The  Input  Parameters 

All  the  vehicle  parameters  read  from  the  data  cards  are  listed 
with  the  appropriate  units. 

d.  Dynamic  Simulation  Results 

Values  of  the  following  variables  after  every  MM  (user  supplied) 
time  steps*  are  printed: 


NO. 

Serial  Number 

TIME 

Time  elapsed  from  beginning  of  the  simulation 

(sec) 

FV2 

Vertical  force  from  each  of  the  front  tires 

(lb) 

FV3 

Vertical  force  from  each  pair  of  the  middle  tires 

(lb) 

FV4 

The  above  for  rear  tires 

(lb) 

FH2 

Horizontal  force  from  each  of  the  front  tires 

(lb) 

FH3 

Horizontal  force  from  each  pair  of  the  middle 
tires 

(lb) 

FH4 

The  above  for  the  rear  tires 

(lb) 

YW2 

Front  axle  displacement 

(ft) 

YW3 

Middle  axle  displacement 

(ft) 

YW4 

Rear  axle  displacement 

(ft) 

*The  time  steps  are  such  that  the  vehicle  advance  by  a distance  dx  at  each 
step.  See  Data  Card  No.  57  for  the  value  of  dx. 
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YO 

Hull  CG  displacement 

(ft) 

THETA 

Hull  pitch  angle 

(degrees) 

PHI 

Bogie  pitch  angle  (forward  pitch  is  positive) 

(degrees) 

DEL2 

Front  suspension  deflection  (zero  when  unloaded) 

(ft) 

DEL3 

Middle  suspension  deflection  (zero  when 
unloaded) 

(ft) 

DEL4 

Rear  suspension  deflection  (zero  when  unloaded) 

(ft) 

Y2 

Terrain  elevation  under  front  wheel  center 

(ft) 

Y3 

Terrain  elevation  under  middle  wheel  center 

(ft) 

Y4 

Terrain  elevation  under  rear  wheel  center 

(ft) 

YB 

Bogie  CG  displacement 

(ft) 

FY2 

Front  suspension  force 

(lb) 

FSY3 

Middle  suspension  force 

(lb) 

FSY4 

Rear  suspension  force 

(lb) 

FFS 

Front  stop  force 

(lb) 

FBS 

Rear  stop  force 

(lb) 

D2Y0 

Hull  CG  acceleration 

(ft/sec^ ) 

D2YB 

Bogie  CG  acceleration 

(ft/ sec^ ) 

D2YW2 

Front  axle  acceleration 

(ft/sec^ ) 

D2YW3 

Middle  axle  acceleration 

(ft/sec^ ) 

D2YW4 

Rear  axle  acceleration 

(ft/ sec^) 

D2  THETA 

Hull  pitch  angular  acceleration 

(rad/sec^) 

D2PHI 

Bogie  pitch  angular  acceleration 

(rad/sec'^) 

At  the  end  of  the  simulation,  the  total  computation  time  needed  for 
the  simulation,  and  to  the  total  number  of  steps  executed  are  written.  A 
typical  printer  output  appears  in  Section  2.3. 

2. 2.  2.  2 Plotter  Output 

Plotter  output  depends  on  the  information  supplied  by  the  user 
on  card  no.  56.  Any  of  the  31  variables  summarized  in  Table  3 can  be 
plotted  against  distance  travelled  or  time.  Additional  variables,  up  to  a 
total  of  40  variables,  can  be  plotted  by  simply  adding  them  to  the  list  in 
subroutine  PSTORE.  Care  should,  however,  be  taken  to  make  the  additional 
variables  available  to  the  subroutine  PSTORE  through  COMMON  statements. 
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The  plots  are  self- scaled  and  of  6 x 5 inches  in  size.  They 
have  a variable  number  written  on  their  Y axis  which  is  keyed  to  the  number 
in  Table  3.  A typical  plot  is  shown  in  Section  2.  3. 

2.  2.  2.  3 Punched  Output 

The  program  generates  punched  output  of  the  time  and  axle 
displacement  data  for  use  in  creating  records  of  shaker  input  signals.  The 
variables  are  punched  in  the  following  order  and  format: 

Variables 

Two  sets  of: 

Time,  YW2  (front  axle  displacement), 

YW3  (middle  axle  displacement), 

and  YW4  (rear  axle  displacement) 

At  the  end  of  the  output,  the  program  punches  a series  of 
minus  ones  to  indicate  end  of  the  record  to  the  control  system  in  the  second 
step  (see  3.2.1).  The  punched  output  can  later  be  converted  to  those  on 
paper  tapes  according  to  equipment  availability  at  the  user's  installation. 

If  the  paper  tape  size  is  not  large  enough  to  include  the  complete  record, 
the  transfer  can  be  continued  on  the  second  paper  tape  by  punching  a series 
of  minus  ones  at  the  end  of  the  first  tape. 

2.2.3  Diagnostic  Message 

If  the  length  of  the  terrain  supplied  by  the  user  is  insufficient, 
the  program  terminates  with  a diagnostic  message  "TERRAIN  LENGTH 
INSUFFICIENT". 

In  such  a case,  the  terrain  length  has  to  be  increased  to  at 
least  the  length  of  the  truck. 


Format 

F7*4,  3F8*5, 
F7-4,  3F8-5 
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2.  3 Illustrative  Example 


The  computer  program  is  used  to  simulate  operation  of  a loaded  M-809 
truck  operating  at  18  mph  over  a 500  foot  long  track  of  1 inch  rms  roughness. 
The  program  employs  the  point  contact  model  to  simulate  the  operation  and 
the  terrain  elevation  data  used  in  the  simulation  have  been  obtained  from  the 
associated  program,  GRND,  described  in  Appendix  A.  3.  The  following 
lists  these  terrain  data  along  with  the  vehicle  and  tire  data  used  to  execute 
the  simulation.  The  output  of  the  program  is  in  three  forms;  a.  printout, 
b.  plots,  and  c.  punched  card  deck  (axle  displacement  records).  Of  these, 
the  printout  and  plots  are  included  in  this  section.  The  punched  deck  con- 
tains the  same  information  as  that  shown  in  the  plots. 
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GX28-7327  6U/M050 

FORTRAN  Coding  Form  Printed  in  U.S.A. 


GX28-7327-e  U/M060 


GX28'732/-6  U/M050 


Input:  Data  Cards  (Card  No.  58  and  onwards) 
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THIS  PKOGrAH  simulates  keave  pitch 
MOTION  OF  A three  axle  VEHICLE  OPERATImO 
OVrR  A SPECIFIED  PROFILE 
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DYNAMIC  SIMULATION 
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Output:  Printout  (first  page  of  dynamic  simulation  results) 
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DYNAMIC  simulation 


Output:  Printout  (last  page  of  dynamic  simulation  results) 
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ADDITIONAL  INFORMATION 


The  simulation  run,  described  in  this  section,  was  made  on  a CDC  6400 
machine  at  the  Smithsonian  Institution  Astrophysical  Observatory  in  Cam- 
bridge, Massachusetts.  The  control  cards  used  were: 


TACOM,  ABBA,  CMIOOK,  T600.  R.  FISH,  FOSTER-MILLER 
RUN(S) 

Note  — >SETCORE  (ZERO) 

LGO. 

CALPLT. 

7/8/9 


Source  Deck  (See  the  listing  in  Appendix  A.  4) 
7/8/9 

Data  Cards 
6/7/8/9 


The  run  used  70K  (octal)  core,  and  109  seconds  of  CPU  time.  The 
cost  for  the  r\m  was  about  $18  (government  rate). 
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3.  STEP  2:  SHAKER  INPUT  SIGNALS  FROM  AXLE 
DISPLACEMENT  RECORDS 


In  the  second  step  of  the  recorded  simulation  control  scheme,  the  axle 
displacement  records  created  by  the  vehicle  simulation  program  are  converted 
to  the  shaker  input  signals  through  two  tasks  shown  in  Figure  9.  The  tasks 
are: 


a.  Acquire  the  data  generated  by  the  first  step  and  store  them  on  an 
on  line  storage  medium. 

b.  Convert  the  records  on  the  storage  medium  to  time  synchronized 
shaker  input  signals. 

For  the  control  system  developed  in  this  effort,  the  first  step  involves 
reading  axle  displacement  data  from  the  paper  tapes  and  transferring  them  to 
a magnetic  disc  unit.  The  data  transfer  is  an  off-line  operation  and  needs  to 
be  performed  only  once  for  any  particular  simulation  data  set  (unless  the 
information  is  destroyed  accidentally).  In  the  data  set  read  by  the  control 
system,  there  are  four  elements:  time,  front  axle  displacement,  middle 
axle  displacement,  and  rear  axle  displacement. 

The  second  task  acquires  the  data  from  magnetic  disc,  converts  the 
displacement  values  to  analog  signals  using  D/A  converters  and  transmits 
them  at  a time  synchronized  rate.  The  time  synchronization  is  achieved  by 
comparing  real  time  signals  produced  by  a crystal  clock  with  the  stored  time 
values.  The  second  task  involves  repeated  playback  of  the  records  repeatedly 
because  the  displacement  records  are  of  a finite  length. 

The  following  describes  the  hardware -software  aspects  of  the  control 
system  and  illustrates  its  use  through  an  example.  Additional  details  of  the 
system  are  in  Appendix  B. 
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Axle 

Displacement 
Records  on 
Paper  Tape 

(a)  Data  File  Creation 


Figure  9.  Details  of  the  Second  Step  of  the  Recorded 
Simulation  Control  Scheme 
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3. 1 Control  System 


3.1.1  System  Hardware 

The  second  step  of  the  recorded  simulation  control  scheme  is 
implemented  using  an  EAI  590  hybrid  computing  system  which  incorporates 
a PACER  100  (16  bit  minicomputer)  interfaced  to  an  EAI  580  analog  com- 
puter, as  shown  in  Figure  10.  The  interface  is  such  that  there  is  a bidirec- 
tional communication  and  data  transfer  between  the  two  computers  during 
parallel  operation.  A magnetic  disc  system  capable  of  storing  two  million 
16  bit  words  is  connected  to  the  digital  processor  via  a direct  memory 
access  channel  (DMAC).  Input-output  by  the  system  is  accomplished  through 
a teletype,  a printer  and  a paper  tape  reader -punch.  Key  characteristics  of 
these  components  are  summarized  in  Table  4,  with  details  of  the  hardware 
and  operations  of  the  components  left  for  the  EAI  publications  listed  in  the 
end  (references  4,  5,  6 and  7). 

For  the  example  used  in  this  manual,  the  displacement  records 
to  be  stored  in  the  disc  are  18  seconds  long,  with  about  5000  time  steps.  The 
records  then  need  about  150K  16  bit  words  of  storage  space  on  the  disc. 

These  data  are  time  synchronized  by  a crystal  clock  prior  to  transmission. 
The  clock  has  a 16  bit  counter  register  and  the  value  stored  in  the  register 
is  compared  by  the  digital  processor  with  the  stored  signal  values.  When 
the  values  become  equal,  the  displacement  data  are  transferred  from  D/A 
converters  to  the  vehicle  shaker  input  electronics.  Signal  outputs  from  the 
D/A  converters  are  bipolar  (signed),  have  a full  scale  reading  of  10  volts, 
and  are  scaled  such  that  a full  scale  output  of  10  volts  equals  1 foot  of  axle 
displacement.  They  are  retrieved  from  the  D/A  output  trunks  on  channels 
one,  two  and  three  and  then  fed  to  the  shaker  servo  inputs  through  any 
required  signal  conditioners. 

In  addition  to  initiating  the  basic  hardware,  the  user  also  has  to 
make  several  hardware  patchboard  connections  on  the  analog  computer.  The 
clock/timer  requires  a connection  between  the  1 KHz  crystal  oscillator  out- 
put and  the  clock's  carry  in  input.  Also  the  control  program,  described  in 
the  next  section,  needs  the  1 millisecond  clock  output  to  be  connected  to  the 
clock  carry  in  input. 
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Figure  10.  Schematic  Diagram  of  the  Control  System  Hardware 


Table  4.  Control  System  Hardware 


PACKR 

100  Computing  System 

• 

16,  384  1 6-  bit  words  of  core  memory 

• 

Hardware  multiply  divide 

• 

Priority  interrupt  system 

• 

Memory  protect  system 

• 

power  fail/auto  restart  system 

• 

Five  universal  controller  positions 

• 

1 Microsecond  memory  cycle 

Direct  Memory  Access  Channel  (DMAC) 

• 

provides  for  I/O  concurrent  with  computation  on  cycle 
stealing  basis 

• 

1 Million  words/second 

Interval  Tiincr/Real  Time  Clock 

• 

provides  a high  resolution  16  bit  real-time  clock 

• 

Timer  update  signal  source  for  four  PACER  interval 
timers 

' Disc  Cartridge  Controller  and  One  Dual  Disk  Drive  on  DMAC 

• 

Provides  1.  1 million  words  of  storage  on  one  fixed  platter 

^ • 

Provides  1.  1 million  words  of  storage  on  one  removable 
platter 

Hybrid  Device  Interface  693 

• 

Parallel  processor  control 

• 

Dogical  control,  sensing  and  interrupts 

• 

Coefficient  device  setting 

• 

Data  monitoring/display 

• 

Eight  channels  D/A  multiplication  (15  bit) 

• 

15  Channels  A/D  conversion  (13  bit) 

Input/Output  Devices 

• 

High  speed  paper  tape  reader,  300  characters/scc. 

o 

High  speed  paper  tape  punch,  120  charactors/sec. 

• 

Teletypewriter  station  (console  device),  ASR  33 

• 

Line  printer,  80  column 

• 

CPU  control  front  panel 

Softwarr- 

System  s 

• 

MHDOS,  moving  head  disc  operating  system 

• 

FORTRAN  IV  compiler 

• 

Assembler 

• 

BTE,  basic  text  editor 

• 

CIG,  core  image  generator 

• 

HOI,  hybrid  operations  interpreter 

• 

COP,  control  o]>tions  processor 

• 

FlU,  file  interactive  utility 

• 

RTL,  run  time  library 

• 

HLR,  hybrid  linkage  routines 

• 

Applications  software  packages 
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System  Software 


Two  programs  are  required  to  operate  the  hardware  described 
in  the  previous  section.  The  first  program,  called  "File  Builder  Program,  " 
reads  the  paper  tapes  containing  the  axle  displacement  records,  packs  the 
data  in  a buffer  and  writes  them  on  a disc  file  as  shown  in  Figure  9(a)  (see 
Appendix  B.2.1  for  details).  The  second  program,  called  "Control  Program," 
retrieves  the  records  from  the  disc  storage  file,  loads  them  into  D/A  con- 
verters and  transmits  them  at  a real-time  synchronized  rate.  This  program, 
which  performs  Task  2 shown  in  Figure  9(b),  is  described  in  detail  in 
Appendix  B.2.2, 

The  file  builder  program  loads  data  in  blocks  of  1 0 time  steps 
into  the  buffer  before  transferring  them  to  a magnetic  disc  storage  file.  The 
loading  is  continued  until  the  end  of  the  record  is  indicated  by  the  buffer 
filled  with  minus  one.  The  data  file  on  the  disc  is  then  rewound  and  the  file 
is  scanned  for  bad  data  points.  A bad  data  point  is  defined  as  the  one  for 
which  the  time  value  is  outside  the  -1  to  20  range  or  an  axle  displacement 
value  is  outside  the  -1  to  +1  foot  range.  (This  definition  is  valid  for  the 
example  described  in  the  next  section  and  may  have  to  change  depending  on 
the  time  length  of  the  records  or  expected  values  of  the  axle  displacements.) 
The  file  is  then  scanned  for  such  bad  data  points,  10  time  steps  at  a time,  in 
the  b\iffer  and  if  a bad  point  is  found  the  10  steps  in  the  buffer  are  printed  on 
the  line  printer  with  an  error  message  requesting  the  user  to  correct  the 
data  point.  The  user  has  to  respond  with  correct  values  for  the  buffer  data, 
which  are  incorporated  by  the  program  through  backspacing  and  overwriting 
the  buffer  back  to  the  disc.  This  scanning  process  is  continued  until  the  end 
of  the  disc  file  is  reached. 

The  control  program  initializes  control  data  and  positions  them 
at  the  start  of  the  sample.  Then  it  waits  until  sense  switch  A,  shown  in 
Figure  11,  is  turned  on  to  begin  transmission.  Once  the  transmission  is 
initiated,  the  program  continuously  cycles  the  records  until  the  user  ter- 
minates the  program  by  turning  sense  switch  B on.  During  data  transmission 
the  program  retrieves  blocks  of  10  time  steps  from  the  disc  and  loads  them 
into  a buffer  in  core  memory.  The  next  operation  is  performed  in  two  steps. 


58 


Figure  11.  PACER  100  Control  Panel 
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In  the  first  step,  the  digital  values  are  loaded  into  a buffer  and  then  supplied  to 
the  D/A  registers.  In  the  second  step,  the  three  D/A  registers,  now  pre- 
loaded  with  the  front,  middle  and  rear  axle  displacement  values,  transfer 
the  signals  through  the  output  circuitry  at  the  correct  time.  After  the  trans- 
fer, the  D/A  registers  are  preloaded  with  the  next  time  step  values  and  the 
program  again  waits  for  the  correct  time.  The  data  transmission  is  continued 
for  the  10  time  steps  before  the  next  block  is  retrieved  from  the  disc.  If  the 
last  block  has  been  transmitted,  the  disc  file  is  rewound  and  the  record  is 
repeated. 


Time  synchronization  mentioned  above  is  achieved  by  a clock- 
timer  which  works  by  counting  down  1 KHz  pulses,  put  out  by  a hybrid  clock, 
from  an  initial  value  to  zero.  The  control  program  compares  the  counter 
value  to  the  corresponding  time  value  stored  with  the  displacement  records 
and  when  the  real-time  value  becomes  less  than  the  stored  value,  the  program 
allows  the  transmission  to  take  place. 

3.2  User  Instructions 


Instructions  for  using  the  control  system  in  terms  of  commands  and 
operations  required  to  execute  the  file  builder  program  and  the  control  pro- 
gram are  described  in  the  following.  Details  of  the  commands  used  in  com- 
munication with  the  operating  system  are,  however,  not  discussed  because 
they  can  readily  be  found  in  the  EAI  reference  manuals  (reference  4,  5,  6 
and  7).  Also  the  procedures  for  program  editing,  compilation  and  load 
module  generation  are  not  discussed  here,  but  are  presented  in  Appendix  B. 3. 

There  are  certain  assumptions  made  in  developing  the  software  and 
user  instructions  for  the  control  system: 

ao  The  time  values  to  be  stored  are  less  than  20  seconds  (i.e., 
record  length  is  < 20  seconds). 

b.  The  values  of  the  axle  displacements  are  between  +1  foot. 
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c.  The  number  of  time  steps  is  less  than  30,  000. 

d.  The  interval  between  the  time  steps  is  more  than  2 milliseconds. 


The  first  two  assumptions  are  used  to  detect  bad  data  points  when  the 
file  builder  program  scans  the  input  data.  The  third  assumption  ensures 
that  the  disc  storage  space  will  not  be  filled  up  and  the  final  assumption  pre- 
vents the  time  synchronizing  devices  from  losing  synchronization.  The  first 
two  restrictions  are  easy  to  change  by  appropriately  altering  the  file  builder 
program.  But  the  last  two  restrictions  cannot  be  changed  easily  because  of 
the  limitations  of  the  existing  hardware. 

3.2.1  File  Builder  Program  Operation 

The  file  builder  program  (described  in  Appendix  B.2.1)  reads 
the  axle  displacement  records  from  paper  tapes,  stores  them  on  disc,  prints 
them  (optional),  and  scans  them  for  bad  data  points.  Before  executing  the 
program,  however,  the  user  must  initialize  a disc  file,  for  the  program  out- 
put, on  logical  unit  21  via  a command  input  to  the  operating  system: 


Example 

where 


#N,  FILEl,  21,  3,  130 


#N 

FILEl 

21 

3 

130 


command  to  monitor  to  create  a disk  file 

name  given  to  file 

logical  unit  file  created  on 

file  type,  data  file  = 3 

words /record  for  data  file  in  octal 

carriage  return. 


The  user  must  then  load  the  core  image  file  of  the  program  from  disc  into 
core  memory; 


Example 


#L,  FILEBLD,  22 


where 

#L  = calls  loader  program 

FILEBLD  = name  of  file  builder  core  image  file 

22  = logical  unit  number  used  by  loader  in  octal.  (Do  not  use 

21,  because  it  is  already  used  in  the  previous  step.) 

Execution  of  the  loaded  program  requires  input  of  the  executions  address: 

Example  #G,  1000 

where 

#G  = command  to  link  to  execution  address 

1000  = standard  execution  address,  octal. 

Control  of  the  computer  is  now  transferred  to  the  file  builder  pro- 
gram and  a message  is  printed  by  the  program  to  notify  the  user  that  it  is 
ready  for  execution.  The  user  then  positions  the  paper  tape  in  the  reader, 
sets  the  data  (sense)  switches  to  implement  the  option  he  wishes,  and  hits  the 
console  carriage  return  to  initiate  execution.  The  message  and  the  data 
switch  options  are: 

"FILE  CREATION  PROGRAM" 

"SET  DATA  SWITCHES  A-D" 

"HIT  CARRIAGE  RETURN  TO  EXECUTE.  " 

User  response: 

DATA  SWITCH  A (ON)  = List  data  as  input 
DATA  SWITCH  B (ON)  = End  program  (at  any  point) 

DATA  SWITCH  C (ON)  = Pause  before  input  each  step 

DATA  SWITCH  D (ON)  = Skip  input  and  go  to  scan  (see  Figure  11). 
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During  execution  of  the  program,  the  axle  displacement  records 
are  read  from  paper  tape  and  stored  in  a biaffer  in  memory.  After  data  for 
10  time  steps  have  been  read,  the  buffer  data  are  transferred  to  the  disc  file 
(formats  for  data  input  and  storage  are  in  Appendix  B.4).  The  program  con- 
tinues the  process  until  it  encoionters  a time  value  of  -1*.  When  a -1  is 
detected,  the  program  asks  the  user  if  he  wishes  to  store  the  value  to  mark 
end  of  the  record  or  treat  it  as  the  physical  end  of  the  tape  and  read  the  next 
paper  tape  as  a continued  record: 

"-1  INPUT,  TYPE  0 TO  SKIP;  ELSE  1" 

User  response: 

1 which  means  store  -1  as  data,  because  it  is  end  of  the  record 

0 which  means  skip  -1,  because  the  record  is  continued  on  the 

next  paper  tape 

At  the  end  of  data  the  user  should  input  enough  records  starting 
with  minus  one  to  fill  out  the  last  buffer.  The  disc  then  encounters  -1  in  the 
last  time  entry,  initiating  the  program  to  write  an  end  of  file  mark  and  close 
the  disc  file.  After  the  file  is  closed  the  program  enters  the  scanning  phase, 
in  which  the  program  checks  the  created  file  for  time  values  outside  the  per- 
missible values  of  0 to  20  seconds,  and  axle  displacement  values  outside  the 
plus  or  minus  one  (foot)  range.  For  scanning  operation  the  data  are  read  into 
the  buffer  and  if  a bad  value  is  detected,  the  data  in  the  buffer  are  printed  on 
the  line  printer.  Also,  a message  is  printed  at  the  console  device  asking  for 
the  bad  data's  buffer  position  and  new  value: 

"INPUT  ID  # OF  BUFFER  TO  CHANGE  AND  NEW  VALUE,  12, 
E12.5" 

"ID.GE.41  STOPS" 


*Using  -1  indicator  this  way  allows  the  user  to  stop  the  tape  input  at  physical 
end  of  a tape  before  the  actual  end  of  the  record. 
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The  user  then  must  input  the  element  number  and  a new  value  to  correct  the 
bad  entry. 

Example  04+0. 12345E+00 

Which  means:  reset  element  No.  4 (rear  axle  displacement  of  the  first  time 
step;  the  element  numbers  are  assigned  sequentially  in  the  data  block)  to 
0.12345.  After  the  correction  is  made,  the  program  prints  the  message  again 
and  the  process  is  repeated  until  the  user  indicates  end  of  correction  by 
assigning  a value  more  than  or  equal  to  41  to  the  element  number.  The  pro- 
gram then  backspaces  the  file  by  one  block,  overwrites  the  bad  data  and 
resumes  scanning.  The  scanning  operation  is  completed  when  end  of  the  file 
is  reached. 

3.2.2  Control  Program  Operation 

The  control  program  obtains  the  displacement  data  from  the  disc 
and  generates  time  synchronized  shaker  input  signals  through  the  digital  to 
analog  converters.  This  program  requires  the  user  to  position  the  data  file 
to  be  transmitted  on  logical  unit  21  (octal)  and  connect  the  hybrid  clock  to  the 
digital  timer  input.  Positioning  of  the  disk  file  is  accomplished  by  a monitor 
command  input  at  the  console: 

Example 

where 

#P  = command  to  open  disk  file 
FILEl  = name  of  data  file  to  be  positioned 
21  = logical  unit  file  opened  on,  octal. 

The  clock/timer  requires  a jumper  plug  be  installed  from  the  1 kilohertz 
output  of  the  hybrid  clock  to  the  clock/timer  carry  in  input.  Then  the  clock 
mode  controls  are  set  to  run  with  a frequency  of  1 megahertz  (details  on 
operation  of  the  hybrid  clock  can  be  found  in  the  EAI  reference  manuals,  i.  e.  , 
references  4,  5,  6,  7 at  the  end. 


#P,  FILEl,  21 


:> 
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Execution  of  the  control  program  requires  loading  of  the  core 
image  file  of  the  program  from  the  disc  and  then  branching  to  the  execution 
address: 


Example 


#L,  CONTPG,  22 
#G,  1000 


(These  arguments  are  described  previously.)  Control  of  the  computer  is  now 
transferred  to  the  control  program,  which  prints  a message  to  the  console 
device  and  begins  testing  data  switch  A for  start  of  transmission: 


"PROGRAM  SHAKER" 
"LU=21  for  data  FILE" 


When  data  switch  A is  turned  on,  the  program  begins  transmission  of  the  three 
axle  displacement  signals,  which  are  brought  out  at  the  D/A  outputs  on  the 
analog  patch  board.  Once  transmission  begins  the  program  will  continue 
cycling  through  the  data  set  until  the  user  terminates  the  program  by  turning 
data  switch  B on,  at  which  point  the  program  prints  a message  on  the  console 
device  and  control  is  returned  to  the  monitor: 

"END  PROGRAM" 

3. 3 Illustrative  Example 

The  axle  displacement  records  generated  by  the  illustrative  simulation 
in  2.3  are  converted  to  shaker  input  signals  through  the  instructions  given  in 
the  previous  section.  These  axle  displacement  records  are  18  seconds  long 
and  include  4780  time  steps. 

The  signal  outputs  from  D/A  are  recorded  on  a strip  chart  recorder 
and  are  shown  in  Figures  12,  13  and  14  along  with  the  corresponding  plots  of 
the  displacements  obtained  from  the  simulation  run.  As  shown  in  the  figures, 
the  control  system  succeeds  in  providing  outputs  which  are  time  synchronized 
and  are  identical  to  the  digital  axle  displacement  records. 
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FRONT  AXLE  DISPLACEMENT  (ft. 


CONTROLLER  OUTPUT 


Figure  12.  Comparison  of  the  Simulation  and  the  Controller 
Outputs  for  the  Front  Axle  Displacement 
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R37.R  AXLE  DISPLACEMENT  (ft 


CONTROLLER  OUTPUT 


TIME  (si:c. ) 

Figure  14.  Comparison  of  the  Simulation  and  the  Controller 
Outputs  for  the  Rear  Axle  Displacement 


APPENDIX  A 


DETAILS  OF  THE  VEHICLE  SIMULATION  PROGRAM 


The  dynamic  simulation  program  used  for  generating  axle  displacement 
records  was  introduced  in  Section  2.  That  section  also  discussed  use  of  the 
program  and  illustrated  the  use  with  an  example.  In  Appendix  A,  further 
details  of  the  computer  program  are  provided. 

A.  1 Principal  Program  Nomenclature 

The  variables  used  in  the  vehicle  simulation  program  are  defined  in 
this  section  along  with  the  corresponding  symbols  used  in  the  analysis 
(reference  1).  All  program  variables  are  in  ft-lb-sec  \anits  except  where 
indicated  to  the  contrary.  (In  the  input-output  operations,  however,  the 
variables  may  not  be  in  the  standard  ft-lb-sec  units,  as  mentioned  in 
Section  2.2.) 

Program 


Variable  Name 

Symbol 

Explanation 

A 

a 

Half  of  bogie  length 

ABC 

b" 

Distributed  (cubic)  tire  damping  constant  in 
the  adaptive  footprint  model 

ADCT(I) 

Cose 

ADDT 

de 

Tire  sector  size 

ADDX(I) 

6x 

Distance  of  the  ith  tire  sector  periphery  point 
from  the  nearest  ground  point  (See  Figure  A- 3) 

ADSDT 

ADST(I) 

Sine 

ADTH(I) 

9 

Angular  position  of  the  ith  tire  sector 
(See  Figure  A-3) 

ADTT(I) 

Tan  9 

AKC 

k” 

Distributed  (cubic)  tire  stiffness  in  the 
adaptive  footprint  model 
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Program 
Variable  Name 

Symbol 

Explanation 

AL(I) 

Contact  length  for  the  ith  sector  (See 

Figure  A- 3) 

ALl(I) 

Contact  length  of  the  part  away  from  the  verti- 
cal for  the  ith  sector  (See  Figure  A-3) 

AL2(I) 

Contact  length  of  the  part  towards  the  vertical 
for  the  ith  sector  (See  Figure  A-3) 

APHI(I) 

<^i 

Ground  slope  corresponding  to  ith  sector 

BS 

B 

Tire  effective  width 

BT 

b 

Tire  damping  constant 

B2V 

B2v 

Shock  absorber  damping  constant 

B2VJ 

B,  . 

2vj 

Shock  absorber  damping  constant,  jounce 

B2VR 

®2vr 

Shock  absorber  damping  constant,  rebound 

CONL 

L 

Equilibrium  tire  footprint  length 

CPHI 

cos<^ 

CTHETA 

cos  © 

DDEL(I) 

5(0) 

Radial  revlative  velocity  for  the  ith  sector 

DDEL2 

'^2 

Suspension  relative  velocity  for  front 
suspension 

DDEL3 

K 

Suspension  relative  velocity  for  middle 
suspension 

DDEL4 

h 

Suspension  relative  velocity  for  rear 
suspension 

DELR(I) 

5(0) 

Radial  deflection  for  ith  sector  (see  Figure  A- 

DEL2 

•^2 

Suspension  displacement  for  front  suspension 

DEL3 

^3 

Suspension  displacement  for  middle 
suspension 

DEL4 

^4 

Suspension  displacement  for  rear  suspension 

DELTA(I) 

Vertical  deflection  corresponding  to  ith  sector 
(See  Figure  A- 5) 
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Program 


Variable  Name 

Symbol 

Explanation 

DER(M) 

Derivative  of  the  state  variables 

DPHI 

Bogie  pitch  velocity 

DTHETA 

© 

Hull  pitch  velocity 

DTIME 

Time  step,  dt 

DYE 

• 

Bogie  heave  velocity  (positive  upwards) 

DYW 

• 

Yw 

Axle  heave  velocity 

DYW2 

^w2 

Front  axle  heave  velocity 

DYW3 

Yw3 

Middle  axle  heave  velocity 

DYW4 

Yw4 

Rear  axle  heave  velocity 

DYX(M) 

Terrain  slope  for  the  Mth  ground  point 

DYO 

Hull  heave  velocity 

DY2 

Rate  of  change  of  ground  elevation  for  front 
wheel 

DY3 

Rate  of  change  of  ground  elevation  for  middle 
wheel 

DY4 

Rate  of  change  of  ground  elevation  for  rear 
wheel 

D2PHI 

Bogie  pitch  acceleration 

D2  THETA 

a • 

© 

Hull  pitch  acceleration 

D2YB 

f a 

Yb 

Bogie  heave  acceleration 

D2YW2 

Yw2 

Front  axle  heave  acceleration 

D2YW3 

Yw3 

Middle  axle  heave  acceleration 

D2YW4 

Yw4 

Rear  axle  heave  acceleration 

D2Y0 

• a 

^0 

Hull  heave  acceleration 

FBT(I) 

Damping  force,  ith  sector 

FBS 

^bs 

Rear  stop  force 

FB2F 

^b2f 

Front  suspension  dry  friction  force 

71 


Program 


Variable  Name 

Symbol 

Explanation 

FB2V 

bZv 

Front  suspension  shock  absorber  force 

FB3F 

^b3f 

Middle  suspension  dry  friction  force 

FB4F 

^b4f 

Rear  suspension  dry  friction  force 

FC2 

Vertical  force  from  each  of  the  front  tires, 
assuming  no  wheel  hop 

FC3 

Vertical  force  from  each  pair  of  the  middle 
tires,  assuming  no  wheel  hop 

FC4 

Vertical  force  from  each  pair  of  the  rear 
tires,  assuming  no  wheel  hop 

FFS 

^fs 

Front  stop  force 

FHG(I) 

Force  component  parallel  to  ground  for  the 
ith  sector  in  adaptive  footprint  model 

FH2 

^h2 

Horizontal  force  from  each  of  the  front  tires 

FH3 

^h3 

Horizontal  force  from  each  pair  of  the  middle 
tires 

FH4 

^h4 

Horizontal  force  from  each  pair  of  the  rear 
tires 

FKT(I) 

Carcass  stiffness  force  for  the  ith  sector 
(See  Figure  A-5> 

FK2 

^k2 

Front  suspension  spring  force 

FK3 

^k3 

Middle  suspension  spring  force 

FK4 

^k4 

Rear  suspension  spring  force 

FNG(I) 

Force  component  normal  to  ground  for  the  ith 
sector  in  adaptive  footprint  model 

FP(I) 

Pressure  force  for  the  ith  sector  (See 

Figure  A- 5)  g 

FSTOP2 

stop2 

Stop  force,  front  suspension 

FSTOP3 

^stop3 

Stop  force,  middle  suspension 

FSTOP4 

^stop4 

Stop  force,  rear  suspension 
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Program 


Variable  Name 

Symbol 

Explanation 

FSX3 

^sx3 

Horizontal  force  at  each  of  the  bogie  and 
middle  suspension  connections 

FSX4 

^sx4 

Horizontal  force  at  each  of  the  bogie  and  rear 
suspension  connections 

FSY3 

syi 

Vertical  force  through  middle  suspension 

FSY4 

^sy4 

Vertical  force  through  rear  suspension 

FTIME 

Simulation  time  limit 

FV2 

^v2 

Vertical  force  from  each  of  the  front  tires 

FV3 

^v3 

Vertical  force  from  each  pair  of  the  middle 
tires 

FV4 

^v4 

Vertical  force  from  each  pair  of  the  rear  tires 

FXX 

Horizontal  tire  force  in  adaptive  foot  model 

FX2 

^x2 

Horizontal  force  at  each  of  the  hull  and  front 
suspension  connections 

FX34 

^x34 

Horizontal  force  at  each  of  the  hull  and  bogie 
connections 

FYY 

F 

V 

Vertical  tire  force  in  adaptive  footprint  model 

FY2 

F 7 
y2 

Vertical  force  through  front  suspension 

FY34 

^y34 

Vertical  force  at  each  of  the  hull  and  bogie 
connections 

G(M) 

Difference  between  terrain  and  wheel  slopes 
for  the  Mth  point  in  the  rigid  tread  band  model 

GB 

Rear  stop  gap 

GF 

Front  stop  gap 

GG 

GG 

Hull  CG  elevation  from  hull  plane 

HS 

h 

s 

Stop  size 

HSTOP 

stop 

Stiffness  factor,  bogie-hull  (See  Note  2 in  2.2) 

IB 

'b 

Pitch  inertia  of  bogie 
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Program 
Variable  Name 

Symbol 

Explanation 

IDAMP(I) 

Ground  contact  index 

IDAMP(I)  = 0 if  ith  sector  out  of  ground  contact 
= 1 if  ith  sector  in  ground  contact 

IH 

'h 

Pitch  inertia  of  hull 

IPHS 

Total  number  of  ground  points  included  in  a 
simulation  phase 

ISCHEM 

Adaptive  model  scheme  number  (see  descrip- 
tion of  subroutine  ADAPT,  A.  2.  11) 

ITIRE 

Tire  model  number 

ITIRE  = 1 Point  contact 

= 2 Rigid  tread  band 
^ 3 Fixed  footprint 
= 4 Adaptive  footprint 

IXAX 

X axis  index 

IXAX  = 1 if  axis  is  TIME  in  plots 

= 0 if  axis  is  XREAR  in  plots 

IXT 

Point  numbers  on  plot 

JFIN 

Total  number  of  simulation  steps 

JGR 

Ground  point  number  corresponding  to  wheel 
center 

JL 

Number  of  ground  points  in  half  of  the  contact 
length 

JMAX 

Total  number  of  grotind  points  in  a simulation 
(e.g.,  5000) 

JNOD 

Ground  point  number  corresponding  to  a tire 
sector 

JN2 

Jn2 

Jounce  clearance,  front  suspension 

JN3 

Jnj 

Jounce  clearance,  middle  and  rear  suspensions 

JPHS 

Number  of  steps  in  a simulation  phase 

JPL 

Total  number  of  points  in  a plot 

JR 

Number  of  ground  points  in  terrain  length 
equal  to  tire  radius 

JREAR 

Serial  number  of  a simulation  step 
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Program 


Variable  Name 

Symbol 

Explanation 

JSTP 

Margin  at  ends  of  a simulation  phase 

K2 

Front  suspension  stiffness 

K3 

K3 

Middle  and  rear  suspension  stiffness 

KPHS 

Ground  point  serial  number  in  a simulation 
phase 

KSTOP 

K . 
stop 

Stop  stiffness  factor,  bogie-stop  (see  Note  2 
in  Section  2.2) 

KT 

k 

Tire  stiffness 

L,2 

^2 

Distance  from  front  suspension  to  hull  CG 

L3 

^3 

Distance  from  bogie  connection  to  hull  CG 

LIMP 

Limit  on  number  of  data  points  in  a plot 

LPHS 

Serial  number  of  a simulation  phase 

MB 

^b 

Bogie  mass 

MH 

Hull  mass 

MFF 

Total  number  of  ground  points  in  tire  contact 
length 

MM 

Number  of  simulation  steps  before  each 
printing 

MW 

Number  of  ground  points  between  rear  and 
front  wheel  centers 

M2 

^2 

Front  unsprung  mass,  each 

M3 

M3 

Middle  unsprung  mass,  each 

M4 

M4 

Rear  unsprung  mass,  each 

NADS 

Number  of  sectors  on  each  side  in  adaptive 
footprint  tire  model 

NAD2 

Total  number  of  tire  sectors  = 2NADS  + 1 

NBACK 

NADS  + 2 

NPLT(M) 

Variable  number  of  the  Mth  variable  to  be 
plotted 
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Variable  Name 

Symbol 

Explanation 

NPLTM 

Total  number  of  variables  to  be  plotted 

NSG(I) 

Difference  between  the  ground  point  numbers 
corresponding  to  the  ith  sector  and  the  wheel 
center 

NW 

Ground  points  between  rear  and  middle  wheel 
centers 

NXT 

Number  of  simulation  steps  skipped  between 
each  point  in  plots 

ODYX(M) 

Slope  of  the  Mth  point  in  unfiltered  terrain 
profile 

OTER(M) 

Elevation  of  the  Mth  point  in  unfiltered  terrain 
profile 

PHI 

<!> 

Pitch  angle  of  bogie  (forward  pitch  is  positive) 

PLT(M,  I) 

Value  of  Ith  point  in  the  Mth  curve  to  be  plotted 

PLSTP 

End  point  of  plot,  as  measured  by  rear  wheel 
center  horizontal  distance 

PLSTR 

First  point  of  plot,  as  measured  by  rear  wheel 
center  horizontal  distance 

PRESS 

Pi 

Tire  inflation  pressure 

RAD 

r 

Tire  radius 

RB2 

Rb2 

Rebound  clearance,  front  suspension 

RB3 

Rb3 

Rebound  clearance,  middle  and  rear  suspensions 

SB 

Sb 

Distance  between  bogie  connection  and  bogie  CG 

SBS 

Sbs 

Distance  between  bogie  connection  and  rear 
stop 

SDYX 

Sum  of  terrain  slopes  under  tire  contact  length 

SFS 

Sr 

fs 

Distance  between  bogie  connection  and  front 
stop 

SN 

n 

Stop  stiffness  factor  for  suspensions  (see 

Note  2 in  Section  2.2) 

SPEED 

V 

Vehicle  forward  speed 
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Variable  Name 

Symbol 

Explanation 

SPHI 

Sin0 

ST(M) 

Terrain  slope  for  the  Mth  ground  point  counted 
from  wheel  rear 

STGAP 

St 

gap 

Unloaded  gap  between  stop  and  bogie 

STER 

Sum  of  terrain  elevations  under  tire  contact 
length 

STHETA 

Sin  © 

SW(M) 

Wheel  slope  for  Mth  point  from  wheel  rear 

S3  4 

S34 

Distance  between  hull  plane  and  bogie 
connection 

TDEF2 

Tire  deflection,  front  tire 

TDEF3 

Tire  deflection,  middle  tires 

TDEF4 

Tire  deflection,  rear  tires 

TER(I) 

^o 

Elevation  of  the  ith  ground  point  in  filtered 
profile 

TERI 

Terrain  elevation  under  a wheel  center  at 
equilibrium 

TERI2 

Terrain  elevation  under  front  wheel  at 
equilibrium 

TERI3 

Terrain  elevation  under  middle  wheel  at 
equilibrium 

TERI4 

Terrain  elevation  under  rear  wheel  at 
equilibrium 

THETA 

© 

Pitch  angle  of  hull  (forward  pitch  is  positive) 

TIME 

t 

Simulation  time 

TRENT 

Total  tire  contact  length  ( = 2AR(I)) 

TRIM 

Tire  deflection  limit 

TSTIF 

T 

stif 

Stiffness  factor  for  tire  (see  Note  2 in 

Section  2,2) 

WEIGHT 

W 

Equilibrium  weight  on  a tire 
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Variable  Name 

Symbol 

Explanation 

WTR2 

Equilibrium  weight  on  front  tire 

WTR3 

Equilibrium  weight  on  middle  tire 

WTR4 

Equilibrium  weight  on  rear  tire 

XC{M) 

Distances  of  the  possible  contact  points  from 
wheel  center  vertical 

XDEL 

dx 

Distance  between  the  successive  ground  points 

XREAR 

X-coordinate  of  rear  wheel  center 

YB 

Yb 

Heave  displacement  of  bogie  CG 

YO 

^0 

Heave  displacement  of  hull  CG 

Y1 

Height  of  wheel  center  from  ground  reference 
plane 

Y2 

Terrain  elevation  at  front  wheel  center 

Y3 

Terrain  elevation  at  middle  wheel  center 

Y4 

Terrain  elevation  at  rear  wheel  center 

YHS 

The  larger  segment  in  partial  contact  case 
(see  Figure  A-4) 

YU 

Initial  height  of  the  wheel  center  from  ground 
reference  plane 

YU  2 

Initial  height  of  the  front  wheel  center  from 
ground  reference  plane 

YU  3 

Initial  height  of  the  middle  wheel  center  from 
ground  reference  plane 

YU  4 

Initial  height  of  the  rear  wheel  center  from 
ground  reference  plane 

YPLT(M) 

Value  of  the  Mth  plot  variable 

YRMS 

RMS  value  of  the  terrain  profile 

YSP2 

^sp2 

Initial  deflection  of  the  front  suspension 

YSP3 

Y o 
sp3 

Initial  deflection  of  the  middle  suspension 

YSP4 

Y A 
sp4 

Initial  deflection  of  the  rear  suspension 
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Variable  Name 

Symbol 

Explanation 

YST2 

Initial  deflection  of  the  front  tires  in  first 
three  tire  models 

YST3 

Initial  deflection  of  the  middle  tires  in  first 
three  tire  models 

YST4 

Initial  deflection  of  the  rear  tires  in  first 
three  tire  models 

YTER(M) 

Possible  wheel  center  elevations  for  rigid 
tread  band 

YW 

yi 

Displacement  of  wheel  center  (axle) 

YW2 

Yw2 

Displacement  of  front  axle  (shaker  input  signal) 

YW3 

Yw3 

Displacement  of  middle  axle  (shaker  input 
signal) 

YW4 

Yw4 

Displacement  of  rear  axle  (shaker  input  signal) 

ZHS 

The  smaller  segment  in  the  partial  contact 
case  (see  Figure  A-4) 
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A. E Subroutine  Description 


The  vehicle  simulation  program  includes  18  subroutines  interfaced  with 
the  main  program  as  shown  in  Figure  5.  In  addition  the  program  needs 
seven  CALCOMP  subroutines  (INITPLT,  AXIS,  SCALE,  NUMBER,  LINE, 
PLOT  and  FIN)  and  one  system  subroutine  (SECOND).  Each  of  these  sub- 
routines perform  specific  functions  and  therefore  the  program  has  a modular 
structure  in  which  changes  can  be  accomplished  just  by  altering  the  relevant 
subroutine  and  not  the  whole  program.  Table  A-1  describes  briefly  functions 
of  each  of  the  subroutines. 

In  order  to  understand  the  programming  details  of  the  subroutines 
dealing  with  the  tire  models  (POINTC,  RTREAD,  FPRINT,  APRINT,  ADORG 
and  ADAPT),  with  the  vehicle  model  (SUSP  and  BOGIE),  with  terrain  profile 
(PROFILE)  or  with  the  system  state  equations  (STEQU)  a study  of  reference  1 
is  a prerequisite.  Once  such  study  has  been  undertaken,  details  of  how  the 
models  are  incorporated  in  the  computer  program  can  be  obtained  from  the 
following.  Also  described  in  this  section  are  subroutines  DYSYS,  RKDIF, 
PSTORE  and  VPLOTS  which  deal  with  dynamic  simulation  and  plotting  part  of 
the  program.  Subroutines  lOLAB,  lOLABD,  IOLAB2  and  PROGIO  deal  with 
the  data  input  function  of  the  program  and  are  self  explanatory.  The  rest  of 
the  subroutines  are  described,  beginning  with  Main  Program. 

A. 2.1  Main  Program 

The  main  program  acts  as  an  interface  for  various  functions  of 
the  program.  The  first  function  of  the  program,  i.e.,  data  input,  is  carried 
out  by  the  main  program  through  subroutines  PROGIO  and  PROFILE.  The 
data  obtained  from  subroutine  PROGIO  are  converted  to  proper  units:  i.e., 
SPEED  is  converted  from  mph  to  ft/sec,  YRMS  from  inch  to  feet,  inflation 
pressure  PRESS  from  psi  to  psf  and  the  distributed  stiffness  AKC  from 
Ib/in.  ^ to  Ib/ft^. 

Index  IPHS  is  then  set  equal  to  the  dimension  of  array  TER  and 
index  LIMP  is  set  equal  to  the  Y dimension  of  array  PLT.  Index  ISCHEM, 
which  deals  with  sophistication  option  in  subroutine  ADAPT  (see  A. 2. 11)  is 
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Table  A-1,  A Summary  of  Subroutine  Functions 


No. 

Subroutine 

Primary  Function 

Group 

1 

lOLAB 

Reads  header  card  and  writes  it 

2 

lOLABD 

Reads  data  in  F format  and  writes  it 

Data 

^read/write 

3 

IOLAB2 

Reads  and  writes  header  card  with  a line 
skip 

4 

PROGIO 

Coordinates  data  input  and  vehicle 

parameter  printout  • 

Data 

input 

5 

PROFILE 

Reads  terrain  profile  and  processes  it 

6 

POINT  C 

Calculates  tire  forces  for  point  contact 
model 

\ 

7 

RTREAD 

Filters  profile  for  rigid  tread  band 
model 

8 

FPRINT 

Filters  profile  for  fixed  footprint  model 

^Tire 

9 

APRINT 

Divides  tire  in  sectors  for  adaptive 
footprint  model 

models 

10 

ADORG 

Coordinates  tire  force  calculation  for 
adaptive  footprint  model 

11 

ADAPT 

Calculates  tire  forces  for  adaptive 
footprint  model 

12 

SUSP 

Calculates  suspension  forces 

1 Vehicle 

13 

BOGIE 

Calculates  bogie  forces 

1 model 

14 

DYSYS 

Coordinates  dynamic  simulation 

15 

RKDIF 

Coordinates  Runge-Kutta  integration 

Dynamic 

algorithm 

^ simulation 

16 

STEQU 

Provides  derivatives  of  state  variables  > 

17 

PSTORE 

Stores  variables  for  plotting 

18 

VPLOTS 

Produces  plots  using  CALCOMP 
subroutines 

^►Plotting 

81 


normally  set  to  1*.  Index  LPHS  is  equal  to  1,  indicating,  thereby,  that 
subroutine  PROFILE  is  called  for  the  first  time. 

The  main  program  calculates  initial  conditions  for  the  dynamic 
simulation  from  geometric  relationships.  Then  the  simulation  is  carried  out 
through  subroutine  DYSYS,  and  at  the  end  of  the  simulation  the  main  program 
terminates  the  run. 

A. 2. 2 Subroutine  PROFILE 


Subroutine  PROFILE  reads  the  terrain  elevation  data  from  data 
cards  and  processes  them  according  to  the  requirements  of  the  tire  model 
being  used  in  the  particular  simulation  run.  The  reading  in  of  the  data  is 
accomplished  in  steps  or  "phases”  in  view  of  the  limited  core  capacity  used 
by  the  program.  The  size  of  the  phases  is  determined  by  dimensions  of 
arrays  TER,  OTER,  DYX  and  ODYX  which  store  various  aspects  of  the  ter- 
rain data.  These  arrays  are  dimensioned  500  so  that  program  can  fit  in  a 
reasonable  amount  of  core  space.  If  the  terrain  data  has  more  than  500 
points,  the  reading  is  temporarily  halted  at  end  of  the  first  500  points  and 
the  simulation  is  executed  until  it  uses  up  the  data  already  read.  Then 
subroutine  PROFILE  is  called  again  to  read  additional  data  and  the  process 
is  repeated  until  the  simulation  is  carried  out  over  the  complete  profile. 

This  phase  simulation  process  is  illustrated  in  Figure  A-1  and  described 
below. 


Subroutine  PROFILE  is  called  initially  from  the  main  program, 
and  subsequently  it  is  called  from  the  dynamic  simulation  routine  DYSYS. 
When  the  subroutine  is  called  for  the  first  time,  or  indicated  by  LPHS, 

Part  1 is  executed. 


*ISCHEM  is  set  to  2 if  more  accuracy  is  required.  However,  if  a cost 
saving  is  desired  while  still  employing  the  adaptive  footprint  model, 
ISCHEM  is  set  to  0. 
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Reading  the  Terrain  Profile 


Part  1 ; The  following  is  accomplished  in  part  1 . 


a.  Terrain  parameters  XDEL  and  JlvlAX  are  read. 

b.  The  numbers  of  profile  intervals  between  the  rear  and  the 
front  wheel  centers  (MW),  and  rear  and  the  middle  wheel 
centers  (NW)  are  determined.  The  numbers  are  then 
adjusted  to  the  nearest  integers  from  the  real  nximbers 
((25).(35)).>-^ 

c.  Simulation  parameter  FTIME  is  determined  from  JFIN, 
which  is  the  total  niimber  of  terrain  ground  points  minus 
the  number  of  points  under  the  truck.  Incremental  time 
(time  step)  DTIME  is  determined  from  the  speed  of  the 
vehicle  SPEED  and  the  interval  between  the  ground  points 
XDEL.  Finding  DTIME  this  way  ensures  that  the  vehicle 
advances  from  ground  point  to  ground  point  as  the  simula- 
tion progresses.  The  initial  time  is  set  to  zero. 

d.  Adaptive  footprint  initialization  subroutine  APRINT  is 
called  if  ITIRE  is  4. 

e.  The  three  more  sophisticated  tire  models  (rigid  tread  band, 
fixed  footprint  and  adaptive  footprint)  require  knowledge  of 
terrain  elevation  at  JSTP  ground  points  on  either  sides  of 
the  wheel  center  ground  point.  The  value  of  JSTP  is  cal- 
culated thus  ((45)-(51)): 

JSTP=  JL=  Integer  (CONL/Z’i'XDEL)  +1  for  fixed  footprint 
model 

JSTP  =»  JR=  Integer  (RAD/XDEL)  +1  for  rigid  tread  band 
and  adaptive  footprint  models 

JSTP  = 0 for  point  contact  model 


^Numbers  in  brackets  refer  to  the  corresponding  line  numbers  on  the  left 
side  of  the  listing  (see  A. 4). 
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The  value  of  JSTP  is  then  used  to  set  margines  on  the  ends 
of  a simulation  phase. 

f.  Because  the  profile  is  read  from  data  cards,  the  number  of 
elevation  points  read  in  each  simulation  phase  has  to  be  an 
exact  multiple  of  the  number  of  data  points  on  each  card 
(say  8).  This  is  ensured  by  ((54)-(59)). 

g.  The  number  of  points  to  be  skipped  while  storing  the  data 
for  plotting  (NXT)  is  determined  from  the  starting  and 
stopping  values  of  the  plot  and  size  of  the  plot  array 
((61)-(62)). 

h.  The  profile  is  read  and  then  converted  to  a profile  with 
correct  RMS  value  ((67)~(69))  (see  Note  1 in  2.2). 

i.  Index  JPHS,  which  determines  the  number  of  steps  the 
dynamic  simulation  is  allowed  to  progress  in  a particular 
simulation  phase,  is  calculated  by  (71).  This  value  for 
JPHS  ensures  that  at  the  end  of  the  simulation  phase,  the 
front  wheel  center  stops  before  the  last  data  point,  with 
JSTP  margine  left  for  computation  (see  Figure  A-1). 

Part  2 ; Part  2 is  executed  only  when  PROFILE  is  called  sub- 
sequently from  DYSYS.  The  major  tasks  of  part  2 are; 

a.  The  part  of  the  profile  under  the  vehicle  is  converted  to  the 
initial  part  of  the  next  phase  profile,  because  the  rear  part 
of  the  truck  has  not  yet  operated  on  it  ((91)-(92)),  This 
conversion  establishes  the  first  MWl  (=M'W+l+2  JSTP) 
points  of  the  next  profile,  as  shown  in  Figure  . 

b.  The  value  of  ITEST,  i.e.,  the  terrain  data  points  to  be 
read  in  this  phase,  is  adjusted  so  that  it  is  divisible  by  the 
number  of  data  points  on  each  data  card  ((78)-(83)). 
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c.  The  value  of  IPHS,  i.e.,  the  total  number  of  ground  points 
in  this  phase,  is  adjusted  so  that  the  number  of  new  profile 
data  read  in  does  not  exceed  the  total  data  points  left  to  be 
read  (85). 

d.  New  terrain  data  is  read  (94). 

e.  Finally  a new  value  of  JPHS,  i.e.,  the  number  of  steps  the 
simulation  is  to  progress  in  this  phase,  is  determined 

((98)-(99)). 

Part  3:  Part  3 is  executed  each  time  subroutine  PROFILE  is 
called,  and  it  accomplishes  the  following  tasks: 

a.  The  terrain  data  are  transferred  to  a new  array,  OTER, 
because  array  TER  stores  the  filtered  profile  for  the  rigid 
tread  band  or  the  fixed  footprint  tire  models  ((1  03)-(l  04)). 

b.  Terrain  slopes  at  each  data  point  is  determined  by  (109) 
and  stored  in  arrays  DYX  and  ODYX.  Array  DYX  stores 
the  filtered  slopes  in  the  next  step. 

c.  If  the  rigid  tread  band  tire  model  is  being  used,  subroutine 
RTREAD  is  called.  If  the  fixed  footprint  tire  model  is 
being  used,  subroutine  FPRINT  is  called.  These  sub- 
routines filter  the  terrain  profile,  and  store  new  elevations 
in  array  TER  and  the  slopes  in  array  DYX. 

A. 2. 3 Subroutine  RTREAD 


Subroutine  RTREAD  filters  the  profile  for  the  rigid  tread  band 
tire  model.  This  task  is  accomplished  by  the  following  steps: 

a.  The  diameter  of  the  rigid  tread  band  is  divided  in  2NS  equal  parts 
of  length  XSEG,  with  value  of  NS  not  less  than  4,  as  shown  in 
Figure  A-2,  ((11)-(15)). 
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b.  For  each  node  on  wheel  periphery,  similar  to  node  I in  the  figure, 
distance  from  the  wheel  center  vertical  is  calculated  (18)  and 
stored  as  XT(I). 

c.  Slope  of  the  tire  periphery  for  each  node  is  calculated  (21)  and 
stored  as  SW(I). 

d.  For  each  node,  the  nearest  groundpoint  towards  rear  is  determined 
and  its  serial  number  difference  from  the  wheel  center  ground  point 
(J)  is  calculated  and  stored  as  JPT(I)  ((24)-(27)).  The  value  of  J 

is  at  least  JR+1,  so  that  the  rear  most  ground  point  needed  for 
calculation  does  not  have  an  unknown  elevation. 

e.  The  slope  of  the  terrain  under  node  I,  ST(I),  is  estimated  by  linear 
interpolation  of  the  slope  values  of  the  surrounding  gro\md  points 
(37). 

f.  The  difference  between  the  terrain  slope  and  the  periphery  slope  is 
calculated  and  stored  as  G(I)  (40). 

g.  The  calculations  e.  and  f.  are  repeated  for  all  2NS+1  nodes 
((32). (54)). 

h.  Index  IS(I)  is  assigned  value  0 for  negative  G(I)  and  value  1 for 
positive  G(I)  ((42)-(44)). 

i.  If  the  sign  of  G(I)  changes  from  positive  to  negative  going  from 
node  I-l  to  node  I,  a contact  point  may  exist  between  the  two  nodes. 
For  all  such  cases,  distances  of  the  possible  contact  points  from 
the  wheel  center  vertical,  XC(K),  are  determined  by  linear 
interpolation  (52). 

j.  The  terrain  slopes  at  the  possible  contact  points  are  also  deter- 
mined from  linear  interpolation  and  stored  as  DYC(K)  (53). 
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k.  Terrain  elevations  at  the  possible  contact  points  are  found  by  linear 
interpolation  from  the  nearest  ground  point  on  the  rear  and  stored 
in  array  TPT  (61). 

l.  Elevation  of  the  wheel  center  for  each  possible  contact  point  is 
found  and  stored  as  YTER(I),  for  I = 1,  K (62), 

m.  The  largest  value  of  YTER(I),  from  the  values  calculated  in  the 
previous  step,  is  the  actual  wheel  center  elevation  and  the  corre- 
sponding contact  point  slope  is  DYC(I).  The  values  of  YTER(I)  and 
DYC(I)  thus  found  are  transferred  to  arrays  TER(J)  and  DYX(J) 
respectively  ((67)-(72)).  These  values  represent  the  filtered  ter- 
rain profile  elevation  (plus  a radius  bias)  and  the  filtered  profile 
slope  corresponding  to  the  ground  point  J. 

n.  The  steps  e.  to  m.  are  repeated  for  ground  points  from  JR+1  to 
IPHS-JR  ((30)-(74)). 

o.  The  elevations  for  the  ground  points  1 to  JR  and  IPHS-JR+1  to  IPHS 
are  assumed  to  be  equal  to  those  for  ground  points  JR+1  and  IPHS- 
JR  respectively.  The  corresponding  slopes  are  assumed  to  be 
zero  ((74)-(81)). 

p.  Finally  the  radius  bias  is  removed  from  the  entire  filtered  profile 
((83)-(84)). 

A. 2. 4 Subroutine  FPRINT 

The  profile  is  filtered  for  the  fixed  footprint  model  by  subroutine 
FPRINT.  This  is  accomplished  by  determining  averages  of  the  terrain  ele- 
vation and  the  terrain  slope  under  footprint  of  the  tire. 

To  initiate  the  filtering  process,  indices  JL  and  MFF  are  calcu- 
lated from  size  of  the  footprint  length  CONE.  The  limits  for  the  calculation 
are  set  from  ground  point  number  JL+1  to  number  IPHS-JL  so  that  there  are 
enough  margins  left  when  the  averages  at  the  end  points  of  a terrain  data 
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record  are  calculated.  Averages  for  terrain  elevation  and  slopes  are  then 
obtained  by  adding  the  ground  point  elevations  and  slopes  under  the  tire  foot- 
print and  dividing  the  sums  by  the  total  number  of  ground  points  under  the 
footprint  - MFF  ((13)-(22)).  The  averages  thus  obtained  are  stored  in  arrays 
TER  and  DYX.  Also,  just  as  for  the  rigid  tread  band  model,  the  filtered 
elevations  for  the  ground  points  1 to  JL  and  IPHS-JL+1  to  IPHS  are  assumed 
to  be  equal  to  those  for  points  JL+1  and  IPHS-JL  respectively.  The  slopes 
for  these  groundpoints  are  assumed  to  be  zero  ((23)-(29)). 

A. 2. 5 Subroutine  APRINT 

When  the  adaptive  footprint  tire  model  is  being  used,  subroutine 
APRINT  is  called  from  subroutine  PROFILE  to  divide  the  wheel  into  pie 
shaped  sectors  and  to  calculate  various  parameters  associated  with  the  sec- 
tors. In  dividing  the  wheel  into  sectors,  only  a part  of  the  wheel  most  likely 
to  be  in  contact  with  the  ground  (say  ±30°  from  the  vertical)  is  considered. 
The  number  of  divisions  are  selected  to  be  19  (=NAD2)  and  the  size  of  each 
sector  then  is  60° /19  = 3*15°.  The  sectors  are  arranged  so  that  the  center 
line  of  the  middle  sector  is  vertical. 

For  each  cf  the  NAD2  sectors,  sector  angle  ADTH(I),  i.e.,  the 
angle  made  by  the  sector  center  line  with  the  vertical,  is  calculated  (see 
Figure  A-3).  Other  parameters  calculated  are  sine  of  the  angle  (ADST(I)), 
cosine  of  the  angle  (ADCT(I)),  tangent  of  the  angle  (ADTT(I)),  number  of 
ground  points  in  the  sector  center  line  projection  (NSG(I))  and  ADDX(I) 

(shown  in  Figure  A-3).  These  parameters  are  constants  and  therefore  they 
are  not  recalculated  during  the  simulation, 

A. 2. 6 Subroutine  DYSYS 


The  dynamic  simulation  part  of  the  program  is  coordinated  by 
subroutine  DYSYS.  Dynamic  simulation  includes  integrating  the  differential 
equations  describing  the  vehicle  behavior  (state  equations)  and  providing  the 
results  in  appropriate  form  to  various  output  devices.  To  perform  integra- 
tion, DYSYS  calls  subroutine  RKDIF  (described  in  A. 2. 7)  which  incorporates 
a fourth  order  Runge-Kutta  numerical  integration  scheme.  The  differential 
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equations  needed  for  integration  are,  however,  stored  in  subroutine  STEQU, 
which  is  called  four  times  every  time  step  by  RKDIF  to  obtain  updated  values 
of  the  derivatives  of  the  variables.  As  the  integration  progresses,  the  re- 
sults are  supplied  to: 

a.  The  printer  (after  every  MM  time  steps) 

b.  Subroutine  PSTORE,  which  stores  values  of  the  variables  to  be 
plotted.  (The  plots  are  later  produced  by  subroutine  VPLOTS 
called  at  the  end  by  DYSYS.  ) 

c.  The  card  punch,  to  prepare  records  of  the  axle  displacements. 

Selected  variables  can  also  be  stored  on  magnetic  tape  by 
adding  a few  cards  as  described  in  A.  3.  These  variables  can  then  be 
processed  by  program  SPEC. 

With  this  introduction,  the  individual  steps  in  the  subroutine 
are  explained  in  the  following 

a.  Equilibrium  deflections  of  the  tires  (YST2,  YST3,  and  YST4)  and 
the  suspensions  (YSP2,  YSP3,  and  YSP4)  are  determined  in 
((29)-(38))  if  any  of  the  first  three  tire  models  are  being  used. 

If  the  adaptive  tire  model  is  being  used  the  procedure  is  as 
follows: 

(1)  Initial  height  of  the  tire  center,  Yl,  is  selected  to  be  equal  to 
two  times  the  radius  plus  the  terrain  elevation  under  the  tire 
center  (TER(JGR)).  (See  Figure  A-3.  ) 

(2)  Subroutine  ADAPT  is  called  to  determine  the  vertical  force, 
FYY,  on  the  tire. 

(3)  The  vertical  force  is  compared  with  the  equilibrium  weight  on 
the  tire,  i.e,,  WTR2,  WTR3  or  WTR4,  depending  on  the  posi- 
tion of  the  tire  in  the  vehicle. 
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(4)  If  FYY  is  less  than  the  weight,  the  tire  center  is  moved  down 
by  0.3  inches  and  the  procedure  is  repeated  from  (2). 

(5)  If  FYY  is  more  than  the  weight,  the  tire  center  is  moved  up 
by  0 . 3 inche  s . 

(6)  Subroutine  ADAPT  is  called  to  determine  FYY. 

(7)  If  FYY  is  less  than  the  equilibrium  weight,  the  tire  center  is 
moved  down  by  0.03  inches  and  the  step  (6)  is  repeated. 

(8)  If  FYY  is  more  than  the  equilibrium  weight,  the  wheel  center 
height  is  stored  as  YI12,  YI13  or  YI14  depending  on  whether 
the  tire  is  in  front,  middle  or  rear  of  the  vehicle. 

b.  System  subroutine  SECOND  is  called  (before  and  after  the  simula- 
tion) to  determine  the  wall  clock  time  required  to  execute  the 
simulation. 

c.  Subroutine  STEQU  is  called  to  initialize  derivatives  of  the  state 
variables  and  other  variables  involved  in  the  state  equations  (such 
as  forces  on  the  vehicle).  These  initial  values  are  then  printed  in 
the  simulation  output.  If  the  user  requires  plots  of  the  variables 
starting  from  the  equilibrium  values,  subroutine  PSTORE  is  also 
called  ((75)-(134)). 

d.  Subroutine  RKDIF  is  called  to  numerically  integrate  the  state 
variables  and  generate  values  of  the  variables  at  the  next  time 
step  (137). 

e.  The  variables  are  printed  (after  every  MM  time  steps)  ((126)-(1 31)). 
The  selected  variables  are  also  stored  for  plotting  ((160)-(169)). 

f.  The  values  of  time  (TIME)  and  axle  displacements  (YW2,  Y'W3  and 
YW4)  are  stored  on  an  intermediate  tape  (TAPE  8)  and  punched 
((122)  and  ((1 79) -(1  91 ))).  The  intermediate  tape  is  used  so  that 
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more  than  one  set  of  data  can  be  stored  on  each  card.  (The  data 
on  these  cards  can  later  be  transferred  to  another  mass  storage 
medium,  such  as  paper  tape  or  magnetic  disc.) 

g.  Steps  d.  to  f.  are  repeated  for  every  time  step  until  the  vehicle 
reaches  the  end  of  the  section  of  the  profile  already  read  (i.e.,  end 
of  a phase  of  the  simulation  is  reached). 

h.  Subroutine  PROFILE  is  called  to  read  additional  profile  data 
((153)-(155)). 

i.  Steps  d.  to  h.  are  repeated  until  end  of  the  profile  is  reached,  as 
indicated  by  TIME  > FTIME  (100). 

j.  If  plots  are  requested,  the  plotting  subroutines  are  initiated  by 
calling  INITPLT.  The  origin  is  fixed  by  calling  PLOT  two  times 
and  subroutine  VPLOTS  is  called  to  generate  the  plots.  The  plot- 
ting is  terminated  by  calling  FIN. 

A. 2. 7 Subroutine  RKDIF 


RKDIF  is  a numerical  integration  subroutine  which  calculates 
the  values  of  the  state  variables  at  time  t+dt,  given  the  values  at  time  t, 
using  a 4th  order  Runge-Kutta  method.  The  integration  scheme  is  sum- 
marized below. 

a.  The  iteration  procedure  starts  with  the  values  of  the  state  variables 
yf,  y2»  etc.,  at  time  t. 

yj(t)  i =*  1,  n 

b.  The  slopes  Dy^(t)  are  then  determined  from  y^(t)  by  calling  STEQU. 
Dy.(t)  = dy.(t)/dt 
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c.  The  values  at  time  t + -g-  are  then  determined, 

t 

Y^l  = Yi  + ^Yi  * dt/2 

d.  The  slopes  Dy^(t  +dt/2)  are  then  determined  by  calling  STEQU 
and  using  the  values  of  y^j^  found  in  c.  above. 

e.  The  values  y^2  time  t + dt/2  are  then  determined 

''a  = + “^11  • 

f.  The  slopes  Dy^2  then  determined  from  STEQU  using 

the  values  of  y^2  found  in  e.  above. 


g.  The  values  y^^  time  t + dt  are  then  determined 


yi3  = I'i  + °yi2  • 

h.  The  slopes  Dy^2  time  t + dt  are  then  determined  from  STEQU 
using  the  values  of  y^^  found  in  g.  above. 


i.  Finally,  the  values  of  the  state  variables  at  time  t + dt  are  found 
as  follows 


y.(t  dt)  = y.(t)  +(Dy.  + 2Dy.j 


+ 2DYi2  ^^i3^ 


A. 2. 8 Subroutine  STEQU 


This  subroutine  incorporates  the  differential  equations  describ- 
ing the  vehicle  motion  (state  equations)  in  terms  of  the  state  variables  and 
forces  on  the  vehicle  components.  In  order  to  update  values  of  the  derivatives, 
subroutine  STEQU  obtains  the  current  values  of  the  state  variables  from  sub- 
routine RKDIF.  From  these  values,  forces  acting  on  the  tires,  suspensions 
bogie  and  hull  are  obtained  by  calling  subroutines  POINTC  (for  the  first  three 
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tire  models),  ADORG  (for  the  adaptive  footprint  model),  SUSP  and  BOGIE. 

The  derivatives  of  the  state  variables  are  then  obtained  from  the  state  equa- 
tions (i.e.,  equation  1,  2,  4,  39,  41  and  43  in  Appendix  B of  reference  1). 

The  accelerations  are  calculated  directly  from  the  derivatives  of  the  state 
variables,  with  the  exception  of  the  bogie  heave  acceleration  which  is 
obtained  from  equation  7 in  Appendix  B of  reference  1. 

A. 2. 9 Subroutine  POINT C 

Subroutine  POINT C determines  the  vertical  and  horizontal  forces 
acting  on  each  tire  using  the  point  contact  tire  model.  For  the  rigid  tread 
band  model  (or  the  fixed  footprint  model)  the  profile  is  first  filtered  by  sub- 
routine RTREAD  (or  FPRINT)  and  then  the  filtered  profile  is  used  to  obtain 
the  forces  from  subroutine  POINTC.  The  task  is  performed  by  the  following 
steps: 

a.  Terrain  elevation  and  time  rate  of  change  of  elevation  under  each 
tire  is  determined  ((20)-(25)).  In  these  equations  index  KPHS 
represents  the  ground  point  under  the  rear  wheel  center.  Adding 
NW  and  MW  to  KPHS  gives  ground  points  under  the  middle  and  the 
front  tires  respectively. 

b.  The  vertical  forces  on  the  tires  are  determined  from  the  point 
contact  model  equations  (see  reference  1,  Appendix  A)  ((2  8)-(44)), 
Equilibrium  tire  deflections  YST2,  YST3  and  YST4  used  in  the 
equations  are  precalculated  in  subroutine  DYSYS. 

Co  The  horizontal  forces  on  the  tires  are  obtained  by  multiplying  the 
vertical  forces  by  the  appropriate  profile  slopes  ((47)-(49)). 

A„2.10  Subroutine  ADORG 


This  subroutine  coordinates  calculations  of  the  vertical  and 
horizontal  tire  forces  for  the  adaptive  footprint  tire  model.  The  force  cal- 
culations are  actually  carried  out  by  subroutine  ADAPT,  but  they  are  in 
terms  of  general  displacement,  velocity  and  force  variables.  These  general 
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variables  are  converted,  in  subroutine  ADORG,  to  specific  variables  for  the 
front,  middle  and  the  rear  tires,  as  summarized  in  Table  A-2.  Of  these 
variables,  Yl  (elevation  of  the  wheel  center  from  the  ground  reference  plane, 
see  Figure  A- 3),  JGR  (ground  point  under  the  wheel  center)  and  DYW  (heave 
velocity  of  the  wheel  center)  are  needed  by  subroutine  ADAPT  to  carry  out 
the  force  calculations.  Therefore,  subroutine  ADORG  assigns  these  general 
variable  names  to  the  specific  variables  for  the  front,  middle  and  rear  tires. 
Subroutine  ADAPT  is  called  next,  for  each  of  the  three  tire  positions  and  the 
resulting  general  force  variables,  FXX  and  FYY  are  converted  to  specific 
variable  according  to  Table  A-2. 


Table  A-2.  Conversion  of  General  Tire  Variables  to  Specific  Variables 


General  Variables 

Specific  Variables 

Front  Tires 

Middle  Tires 

Rear  Tires 

FYY 

FV2 

FV3/2>!' 

FV4/2* 

FXX 

FH2 

FH3/2>:' 

FH4/2* 

YW 

YW2 

YW3 

YW4 

DYW 

DYW2 

DYW3 

DYW4 

YU 

YI12 

YU  3 

YU4 

JGR 

KPHS+MW 

KPHS+NW 

KPHS 

'■'Forces  FV3,  FV4,  FH3  and  FH4  are  for 

tire  pairs . 

A. 2. 11  Subroutine  ADAPT 


Subroutine  ADAPT  determines  the  horizontal  and  the  vertical 
forces  acting  on  tires,  using  the  adaptive  footprint  tire  model.  There  are 
two  parts  to  the  calculations: 

a.  Geometric  calculations 

b.  Force  calculations 
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The  geometric  calculations  deal  with  obtaining  radial  deflections 
and  contact  lengths  for  the  tire  sectors  in  contact  with  the  ground.  The  force 
calculations  use  the  geometric  parameters  thus  determined  to  obtain  various 
force  components  for  the  tire  sectors  and  then  the  total  horizontal  and  verti- 
cal forces  at  the  tire  center. 

Before  commencing  the  calculations,  the  storage  arrays  are 
cleared  and  index  IBACK  is  set  to  zero.  Setting  IBACK  to  zero  indicates 
that  the  calculations  will  start  from  the  bottom  sector  (1=1)  and  progress 
towards  the  front  of  the  tire  (see  Figure  A-3).  After  all  the  sectors  in  contact 
with  the  ground  are  considered,  IBACK  is  set  to  1 and  the  calculations  pro- 
ceed along  the  rear  of  the  tire. 

Part  1 - Geometric  Calculations  ((47)-(109))  - There  are  three 
choices  of  accuracy  available  to  the  user  in  the  geometric  calculations.  The 
choice  is  indicated  by  selecting  a value  of  ISCHEM.  For  ISCHEM=0,  an 
algorithm  with  low  accuracy  is  used,  where  as  ISCHEM=2  indicates  choice  of 
an  accurate  but  expensive  algorithm.  Normally  ISCHEM=1  is  used,  which 
means  that  an  algorithm  with  an  in  between  accuracy  and  cost  will  be  used. 

ISCHEM— 0;  The  geometric  relations  used  in  obtaining  the 
radial  deflection,  DELR(I),  and  the  contact  length,  AL(I)  (see  Figure  A-3), 
for  sector  I are  approximate. 

DELR(I)  = DELTA(I)«  ADCT 

AL(I)  = (RAD-  DELR(I))>:=  ADDT 

These  calculations  also  presume  that  the  ground  contact 
is  continuous  and  includes  the  bottom  point  of  the  tire. 
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ISCHEM=gl ; More  complex  and  accurate  relations  are 
used  in  determining  the  radial  deflections  and  the  contact  lengths  of  the  tire 
sectors.  The  accuracy  is  also  improved  by  including,  in  the  contact  length 
calculations,  the  sectors  which  are  only  partly  in  ground  contact.  However, 
it  is  still  assumed  that  the  ground  contact  is  continuous  and  includes  the 
bottom  point  of  the  tire, 

ISCHEM=2 : This  scheme  is  identical  in  all  respects  to 
the  above  scheme,  except  that  the  ground  contact  may  now  be  discontinuous 
and  may  not  include  the  bottom  point  of  the  tire. 

The  steps  involved  in  the  geometric  calculations  are  then 
summarized  in  the  following: 

a.  Sector  number  I is  set  to  1 (bottom  sector), 

b.  Ground  point  number  corresponding  to  each  sector  is 
determined  and  stored  as  JNOD, 

c.  Slope  of  the  profile  under  the  sector  is  determined 
from  the  stored  data  ((50)-(52)),  If  the  sector  is  on 
the  forward  side,  the  slope  is  calculated  correspond- 
ing a ground  point  ahead  of  the  sector  groxond  point 

to  improve  accuracy, 

d.  Vertical  deflection  DELTA(I)  is  calculated  (55)  from 
the  geometric  relations  between  the  variables  shown 
in  Figure  A- 3, 

e.  Radial  rate  of  deformation  of  the  sector,  DDEL(I)  is 
determined  by  subtracting  the  radial  velocities  of  the 
two  ends  of  the  sector  (56), 
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f.  If  ISCHEM=0,  the  radial  deflection  and  the  contact 
length  are  found  from  the  approximate  model 
{(78)-(81))  and  the  calculations  are  moved  to  step  1. 
For  this  approximate  model,  APHI(I)  is  considered 
to  be  equal  to  ADTH(I),  and  the  sector  is  considered 
to  be  out  of  ground  contact  if  DEIjR(I)  is  negative 
(see  step  p. ). 

g.  For  any  of  the  following  conditions,  the  accurate 
geometric  calculation  algorithms  tend  to  become 
unstable  and  therefore  the  calculations  are  performed 
according  to  the  approximate  model. 

(1)  ||APHI(I)  - ADTH(I)|  - Tr/2|  <0*1  (radian) 

(2 ) II APHI(I)  - ADTH(I)  + ADDT  /2  | - it  /2  | <0*1 

(3)  ||APHI(I)  - ADTH(I)  - ADDT/2|  - Tr/2|  < 0«1 

h.  For  1=1,  i.e.,  the  bottom  sector,  DELR(I)  is  equal 
to  DELTA(I)  (61).  For  1=^1,  DELR(I)  is  ND  in  Figure 
A- 3 and  it  is  determined  by  geometric  relationships 
in  triangles  BNZ  and  DVZ  (63), 

i.  If  DELR(I)  is  positive,  IDAMP(I)  is  set  to  1 so  that 
the  sector  will  be  considered  in  damping  force  calcu- 
lations. If  DELR(I)  is  negative,  the  sector  is  con- 
sidered for  partial  contact,  as  described  in  step  m. 

If  DELR(I)  is  more  than  the  radius  of  the  tire,  the 
accurate  algorithm  is  considered  unstable  and  the 
approximate  algorithm  described  in  step  f.  is  used. 

j.  The  contact  lengths  ALl(I)  and  AL2(I)  (i.e.,  SN  and 
QN  in  Figure  A- 3)  are  determined  from  geometric 
relationships  in  triangle  OQS  ((70)-(7l)).  If  either 
ALl(I)  or  AL2(I)  is  more  than  DELR(I),  the  sector  is 
considered  to  be  in  partial  contact  as  described  in 
step  m. 
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k.  The  contact  length  AL(I)  is  found  by  adding  ALl(I) 
and  AL2(I). 

l.  The  value  of  I is  incremented  by  1,  and  if  I is  less 
than  the  total  number  of  sector,  NAD2,  the  procedure 
is  repeated  from  step  b. 

m.  The  sectors  which  have  exhibited: 

(1)  DELR(I)  negative 

(2)  Either  ALl(I)  or  AL2(I)  less  than  DELR(I) 

are  considered  for  partial  contact  ((97)-(l 09)).  For 
determining  the  partial  contact  lengths,  difference 
between  APHI(I)  and  ADTH(I)  is  stored  as  DIFPT. 

n.  If  DIFPT  is  negative  and  (1)  DELR(I)  is  negative,  the 
sector  is  considered  out  of  ground  contact  and  process 
moves  to  step  p.  (2)  DELR(I)  is  positive,  the  sector 
is  considered  to  be  in  full  contact  and  the  process 
moves  to  step  k. 

o.  Terms  ZHS  and  YHS  are  determined  from  the  geo- 
metric relations  shown  in  Figure  A-4.  If  ZHS  > YHS, 
the  sector  is  considered  to  be  completely  in  groxmd 
contact  and  the  process  moves  to  step  k.  Otherwise, 
AL(I)  is  found  by  adding  ZHS  and  YHS  and  dividing  the 
sum  by  Cos  (DIFPT).  If  AL(I)  is  found  to  be  negative, 
the  sector  is  considered  out  of  ground  contact  and  the 
process  moves  to  step  p. 

p.  If  there  is  no  sector-ground  contact,  DELR(I)  and 
AL(I)  are  set  to  zero.  In  case  of  ISCHEM=0  or  1, 

no  sector-ground  contact  means  that  there  is  no  more 
contact  in  front  of  the  tire  and  the  calculations  are 
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Figure  A-4.  Tire  Sector  in  Partial  Contact 
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repeated  for  the  rear  sectors  (I=:NBACK)  by  starting 
from  step  b.  ((93)-(94)).  If,  however,  the  rear  sec- 
tors have  also  been  considered,  as  indicated  by 
IBACK=1  (91),  the  calculations  for  Part  1 are  com- 
pleted and  Part  2 is  executed  as  described  below. 

For  ISCHEM=2,  no  such  end  of  contact  is  assumed 
and  all  the  sectors  are  considered  for  possible 
multiple  contacts.  Therefore,  the  process  is  trans- 
ferred to  step  1. 

Part  2 - Force  Calculations  ((1 1 3)-(l 31 ))  - After  the  geometric 
calculations  are  completed.  Part  2 of  the  subroutine  is  executed.  The  force 
calculations  considered  in  Part  2 involve  determining  the  pressure,  stiffness 
and  damping  forces  for  each  sector,  finding  their  components  perpendicular 
and  horizontal  to  the  ground,  eliminating  the  negative  perpendicular  forces 
and  finally  determining  the  inertial  vertical  and  horizontal  forces  (FYY  and 
FXX  in  Figure  A- 5),  The  steps  in  detail  are: 

a.  Forces  FXX,  FYY  and  contact  length  TLENT  are  set 
initially  to  zero.  These  parameters  are  determined  by 
performing  addition  of  the  force  components  and  the  sector 
contact  lengths  using  a DO  loop  ((1 1 6)-(l  31 )). 

b.  Pressure  force  for  a particular  sector,  FP(I)  (see  Figure 
A- 5)  is  determined  by  multiplying  the  tire  pressure  to  the 
sector  contact  area  (118). 

c.  Carcass  stiffness  force,  FKT(I),  is  determined  by  inte- 
grating the  incremental  stiffness  force  along  the  radial 
deflection  (121). 

d.  Damping  force,  FBT(I)  is  determined  from  the  distributed 
damping  constant  ABC  and  the  rate  of  sector  deformation 
calculated  in  step  e.  of  Part  1.  If  the  sector  is  not  in  the 
ground  contact,  IDAMP(I),  and  subsequently  FBT(I),  are 
zero  (125). 
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Figure  A- 5.  Tire  Forces  in  the  Adaptive  Footprint  Model 


104 


e.  Total  sector  force  normal  to  the  ground,  FNG(I),  is  found 
by  adding  the  components  of  the  pressure,  stiffness  and 
damping  forces.  Since  only  the  compressive  forces  can  be 
transmitted  between  the  tire  sector  and  the  ground, 
FNG(I)=0,  if  FNG(I)<0. 

f.  Total  sector  force  horizontal  to  the  ground,  FHG(I)t  is  found 
similarly. 

g.  Vertical  sector  force  is  found  from  components  of  the 
forces  vertical  and  horizontal  to  the  ground.  The  sector 
forces  are  added  to  give  FYY. 

h.  Horizontal  sector  forces  similarly  give  FXX. 

The  end  products  of  Part  2,  and  that  of  the  subroutine,  are  the 
vertical  force,  FYY,  and  the  horizontal  force,  FXX,  exerted  by  the  tire  on 
the  axle. 

A. 2. 12  Subroutine  SUSP 

Subroutine  SUSP  calculates  the  forces  in  the  suspensions 
according  to  the  following  steps; 

a.  Displacements  and  velocities  of  both  ends  of  the  suspensions  are 
required  to  calculate  the  forces.  Those  for  the  bottom  ends  are 
available  directly  from  the  state  variables  YW2,  YW3,  YW4  and 
DYW2,  DYWS,  DYW4.  The  displacements  and  velocities  for  the 
top  ends  are  calculated  from  equations  13,  16,  23,  24,  27,  28,  34 
and  37  in  Appendix  B of  reference  1. 

b.  Leaf  spring  stiffness  forces  are  calculated  from  equations  10,  11, 
20,  21,  31  and  32  of  the  reference  report  and  leaf  spring  dry 
friction  forces  are  calculated  from  equations  14,  15,  25,  26,  35 
and  36.  It  should  be  noted  that  the  dry  friction  factor,  FRSF,  is 
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assumed  to  be  0.  05,  which  means  that  the  dry  friction  force  is 
5 percent  of  the  instantaneous  sprint  force  and  opposes  the  velocity. 

c.  Stop  contact  forces  are  calculated  from  equations  18,  29  and  38. 

d.  Viscous  damper  force  for  the  front  suspension  is  calculated  from 
equation  17. 

e.  Finally  the  forces  are  added  to  give  the  total  suspension  forces 
according  to  equations  9,  19  and  30.  The  horizontal  forces  are 
obtained  from  equations  40,  42  and  44. 

A. 2. 13  Subroutine  BOGIE 


Subroutine  BOGIE  calculates  the  frame  stop  forces  and  the  forces 
acting  on  the  bogie. 

a.  The  stop  forces,  determined  from  equations  45,  46,  47,  and  48, 
are  assumed  to  act  on  the  unsprung  masses.  This  is  so  because 
the  bogie  in  the  model  represents  just  the  inertial  properties  of 
the  suspension  springs.  The  stiffness  and  the  damping  properties 
of  the  spring  are  represented  by  the  middle  and  the  rear  suspen- 
sions. The  stop  forces  act  on  the  ends  of  the  spring  and  therefore 
on  the  unsprung  masses  at  the  other  ends  of  the  suspensions  in  the 
model. 

b.  The  horizontal  and  the  vertical  forces  transmitted  by  the  bogie  to 
the  hull  are  calculated  using  equations  5 and  8,  respectively. 

A. 2. 14  Subroutine  VPLOTS 


Subroutine  VPLOTS  includes  the  CALCOMP  plotting  subroutines, 
which  are  executed  to  produce  plots  of  the  variables  stored  by  subroutine 
PSTORE  in  array  PLT.  The  data  for  the  Y axis  of  the  plots  are  obtained  from 
PLT  (IPLN,  I),  where  IPLN  is  the  number  of  the  plot.  These  data  are 
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transferred  to  array  A(I)  in  subroutine  VPLOTS.  Similarly  the  X axis  data 
stored  in  PLT(9,I)  are  transferred  to  array  B(I).  Number  YPN,  which 
represents  the  particular  variable  being  plotted  according  to  the  list  in 
Table  3,  is  written  on  the  Y axis  of  the  plot. 

The  rest  of  the  subroutine  consists  of  calling  the  standard 
CALCOMP  subroutines  SCALE,  AXIS,  NUMBER,  LINE  and  PLOT. 

A. 2. 15  PSTORE 

Subroutine  PSTORE  stores  the  variables  to  be  plotted,  which 
are  selected  by  the  user  by  punching  1 in  appropriate  columns  of  card  No.  56 
(see  2.2).  The  column  numbers  correspond  to  the  variable  numbers  in 
Table  3,  and  are  stored  in  array  NPLT(J),  J=l,  NPLTM,  where  NPLTM  is 
the  total  number  of  plots  desired. 

When  PSTORE  is  called,  the  system  variables  available  for 
plotting  are  converted  to  an  appropriate  member  in  array  YPLT(M),  where  M 
corresponds  to  the  variable  number  in  Table  3,  Then  if  the  Mth  variable  is 
to  be  plotted,  YPLT(M)  is  converted  to  PLT(J,IXT)  where  J=l,  NPLTM  and 
IXT  is  the  plot  point  number  supplied  by  DYSYS. 

Simultaneously  X axis  data  are  stored  in  array  PLT  (9,  IXT). 

The  X axis  can  either  be  XREAR  (the  X coordinate  of  rear  wheel  center)  or 
TIME.  The  user  makes  the  selection  by  punching  appropriately  in  column  41 
of  card  56,  and  thereby  setting  the  value  of  IXAX. 
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A.  3 Associated  Programs 


The  programs  described  in  this  appendix  can  be  used  in  conjunction 
with  the  vehicle  simulation  program  as  described  in  Figure  3.  Program 
GRND  can  be  used  to  generate  a profile  if  such  data  are  not  available, 
whereas,  program  SPEC  can  be  employed  to  determine  spectral  densities 
of  various  variables  associated  with  the  simulation. 

A. 3.1  Program  GRND 

A. 3. 1.1  Description 

This  program  generates  a ground  profile  in  terms  of  eleva- 
tions at  specific  intervals.  The  spectral  density  of  the  generated  profile  is 
given  by 


S (W  = A/n^;  < Q < pL 

Yq  Aj  A^ 

= 0 elsewhere 

where  is  a reduced  frequency  given  by: 

Q = W/v  = Ztt/A 

and  and  cut-off  wavelengths,  selected  to  be  X^^  = 57  feet  and  ^2~ 

2 inches. 


From  the  definition  of  power  spectral  density 

00 

= r s (n)  dfi 

00  ^0 


108 


The  details  of  the  computation  required  to  generate  the  pro- 
file is  available  in  Appendix  C of  reference  2 and  therefore  it  is  not  repeated 
here. 


The  parameters  XX,  XXl , XX2,  YY,  YYl  and  YY2  in  the 
program  are  random  numbers,  selected  from  random  number  tables.  These 
random  numbers  are  used  to  obtain  PHI(K)  in  reference  2)  and  THETA(K) 

(=  20/(i2j^*Ai2j^)).  The  discrete  reduced  frequencies  are  calculated 

such  that  the  log  scale  from  = 0.11  to  “ 35.4  is  divided  equally  in 

149  parts.  For  each  frequency,  the  profile  spectral  density  (A/Cl^)  is 

determined  and  stored  as  S(K). 

Subroutine  EVALFN  is  called  for  each  of  the  5000  ground 
points  to  calculate  terrain  elevation  using  the  150  discrete  frequencies 
obtained  earlier.  The  calculated  elevations  are  punched  on  data  cards  in  an 
8ElO*3  format.  The  mean  square  value  of  the  profile  (SIGMA2)  is  calculated 
by  numerical  integration  and  it  is  printed  to  help  confirm  validity  of  the 
generated  profile.  Also  printed  are  the  total  number  of  elevation  values 
which  lie  in  each  of  the  40  bands  of  elevation  magnitude,  spread  between  -Za 
{a  = RMS  value  of  the  profile  = ^/  ) to  +2a.  This  table  is  useful  in  deter- 

mining probability  distribution  of  the  profile. 

A,  3. 1 .2  Inputs 

There  are  essentially  only  two  inputs  to  the  program, 
assuming  that  the  profile  spectral  density  shape  is  not  changed. 

a.  XDEL,  i.e.,  the  interval  at  which  the  elevations  are  determined. 

It  is  selected  to  be  0.1  foot. 

b.  A,  i.e.,  the  ground  roughness  parameter.  It  is  selected  as  0.1105, 
so  that  the  RMS  value  of  the  profile  elevation  is  1 foot.  Since  the 
terrain-tire -vehicle  simulation  program  converts  the  generated 
profile  to  a profile  with  the  user  selected  RMS  value,  the  value  of 

A used  in  the  profile  generation  program  may  not  be  changed. 
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A. 3. 1.3  Outputs 


The  following  are  the  program  outputs. 

a.  A deck  of  cards  with  the  terrain  elevation  coordinates  in  SElOO 
fo  rmat , 

b.  Mean  square  value  of  the  profile  elevation. 

c.  The  distribution  of  profile  elevation  (G(M)). 

A.  3. 1 .4  Listing 

Printed  on  the  following  pages. 
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SKKI  I <Z.»S  CKI  •0Ml,l*0MEG4  CKI  »•» 

CONTINUE 

DO  2 NN=1,5I)00 

H=NN-1 

XtXOEL»N 

call  EVALFN (FCMMI , DFOX (MM» ,XI 
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9'/7C3103 
'i  If') 

33.00600000 
35.00300000 
51 .00000003 
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A.  3.2  Program  SPEC 


A. 3. 2.1  Description 

Program  SPEC  employs  a fast  fourier  transform  algorithm 
to  generate  spectral  densities  of  different  variables  from  their  time  histories, 
which  are  stored  on  a magnetic  tape  during  execution  of  the  vehicle  simulation 
program  (see  Section  A.  3.2.2  for  details).  These  time  histories  are  read 
from  the  tape  and  stored  in  array  F(I).  These  data  are  processed  in  the  fol- 
lowing ways  before  determining  the  spectral  densities: 

a.  The  number  of  discrete  values  in  each  time  history  is  reduced  to 

4096,  because  it  is  more  convenient  to  process  2^  points  using  the 

12 

fast  fourier  transform  algorithm.  Usually  the  last  4096  (=2  ) 

points  are  selected  so  that  effects  of  the  start  up  transients  are 
reduced  while  performing  the  fourier  analysis. 

b.  The  time  histories  are  processed  by  a cosine  taper  window  (see 
p.  325  of  reference  3)  to  reduce  effects  of  leakage  and  the 
processed  data  are  stored  in  array  A(I). 

After  processing  the  time  histories,  subroutine  FFT  is 
called  to  carry  out  the  fourier  transformation.  The  results  are  returned 
through  array  A(I)  in  form  of  complex  fourier  transforms.  The  transforms 
are  converted  to  spectral  densities  for  each  of  the  4096  frequencies  using 
equation  9.134  of  reference  3,  This  information  is  condensed  in  terms  of 
spectral  densities  for  a set  of  discrete  frequencies  FR(I)  which  are  evenly 
spread  on  a logarithmis  scale  from  wavelength  XDEL  (profile  elevation 
interval)  to  4096  XDEL.  For  this  condensing  operation  the  spectral  densities 
around  each  of  these  discrete  frequencies  are  averaged.  All  the  spectral 
densities  are  also  divided  by  a factor  0.875  to  compensate  for  the  effects  of 
the  cosine  taper  window  (reference  3,  p.  327). 

Finally,  an  approximate  value  of  RMS  for  the  variable  time 
history  is  calculated  by  adding  average  areas  under  the  spectral  density 
curve. 
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A.  3. 2 . 2 Input 


The  input  to  the  program  is  in  form  of  a tape  and  a data  card 


a.  Tape  - This  tape  stores  time  histories  of  the  variables  generated 
by  the  vehicle  simulation  program.  It  can  be  created  by  inserting 
statements,  such  as  the  following  statements,  at  the  two  places 
marked  by  commend  cards  in  subroutine  DYSYS  of  the  program. 


Place  1 


C STftTI'iriCAL  PR0C':SSIM3  OPTION  , SEE  DESCRIPTION  OF  PROGRAM  SPEC  IN 
C THE  USER  MANUAL  FOR  THE  CARDS  TO  QE  INSERTED  HERE  FOR 
C EXGERCISING  thF  OPTION 

C NCJR  IS  THE  TOTAL  NUMDER  OF  VARIABLES  SELECTED 
35  NCUR=10 

FV42:FVA/2. 

FH42=FH4/2. 


WRITE  (9,173)NCUR,  SPEED,  ITIRE,  YRMS 
WRITE  (9,  174)  YW2,YW3,  YW4,  FV2,  FV42,  FH2,  FH42,  D2YW2, 
D2YW4,  D2Y0 

173  FORMAT  (IX,  12,  ElO. 4,12,  Flo. 4) 

174  FORMAT  (IX,  10E12.4) 


Where  tape  9 is  used  as  a storage  tape  and  it  is  specified  as  so  in 
the  program. 


Place  2 


C STATISTICAL  PROCESSING  OPTION  , SEE  DESCRIPTION  OF  PROGRAM  SPEC  IN 
140  C THE  USER  MANUAL  FOR  THE  CAROS  TO  BE  INSERTED  HERE  FDR 

C EXCERCISING  THP  OPTION 
FV42=FV4/2. 

FH42=FH4/2. 


WRITE  (9,  174)  YW2,  YW3,  YW4,  FV2,  FV42,  FH2,  FH42,  D2YW2, 
D2YW4,  D2Y0 

b.  Data  Card  - The  total  number  of  points  in  each  of  the  time  histories 
is  needed  in  l6  format  (this  number  is  available  as  TOTAL  STEPS 
in  the  vehicle  simulation  program  output). 
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A.  3.2. 3 Output 


The  program  output  consists  of  (a)  the  spectral  density  values 
for  a set  of  reduced  frequencies,  (b)  the  corresponding  values  of  spectral 
densities  for  the  frequencies  in  Hertz,  (c)  log  values  of  all  the  spectral 
densities. 


The  printout  is  repeated  for  each  of  the  NCUR  variables. 

On  the  top  of  each  set  tire  model  no.,  speed,  RMS  value  of  the  terrain  profile, 
and  the  variable  serial  number  are  printed.  The  approximate  RMS  value  of 
each  variable  is  printed  at  the  bottom  of  the  set. 

A. 3. 2. 4 Listing 

The  following  is  a listing  of  the  program. 
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PROGPftM  TIRF ( INPUT, OUTPUT, PUNCH, TAPE5I INPUT .tapes: OUTPUT, TAPE/ t 
IPUNCH, TAPES ) 

C*»»»****TERRAIN  TIRE  VEHICLE  SIMULATION  PROGRAM  »•*******•»•* 


A.  4 Listing  of  the  Simulation  Program 


O 

V. 

a 


» o o 

» CD  (/) 

* (/) 


» Ql  O 

* >-  UJ  N 

* CD  -I  fw) 

» -1  f 

* o M a 

» bj  z;  js 

» fi,  I « 
» o 

» u UJ  ^ 

* UJ  ^ N. 

* > M-t 

» tu  o ^ 

* o u. 


ft 

•> 

•> 

n 

UI 

<r  X 

ft 

B 

3 

1-  j-  to 

(O 

X 

> 

UJ  I-  u 

X 

ft 

X to  • 

z 

UI 

k- 

fO 

^ X »Yj 

to 

z 

Z -ft 

z 

3 

(SI  • X 

X 

ft 

kH  O 

UJ 

>- 

O »o  (O 

ft 

X 

1-  -t 

O 

• k-  u. 

< 

X 

k-<  ft- 

k~ 

z 

l-i  (O 

ft  k- 

UJ 

UI 

1- 

X >-  <SJ 

X 

X 

X 

z o 

> 

to 

X 

z 

CL  ••  X 

to 

«* 

k^  O 

>■ 

o 

3 

fi 

w 

(SI  (SI  U 

ft 

X 

X 

z 

ft. 

B k«  *• 

to 

t-H 

(O  ft 

ft 

X 

O' 

o 

<\j 

• to 

(D 

ft 

ft 

k'  -ft 

X 

o 

o 

K 

X 

o > > 

to 

mi 

ft  o 

X 

3 

u 

(O 

V 

• fO 

X 

0. 

X Jt 

z 

z 

X 

Q 

rsj  u 

(0 

X 

3 

o 

kH  ft- 

* 

1- 

to 

» 

• 

Oft*' 

u. 

z 

ft 

ft 

US 

mJ  1- 

z 

o 

fo 

<VJ 

» to  fO 

CO 

> 

z 

z 

k-  to  CO 

o 

UJ 

z 

X 

B > >- 

ft 

CSJ 

< 

< 

X 

ft  X o 

-> 

-1 

o 

«• 

> 

V,  • to 

CD 

UJ 

UI 

V 

(O  z « 

ft 

X 

13 

M 

» 

(SI  fo  u. 

w 

ft 

z 

z 

o 

to  “>  ft 

< 

UJ 

-1 

o 

a ft  » 

to 

3 

3 

X 

o 

Ui  ftpft 

>• 

{S\ 

z 

> 

• to  (Si 

ft 

ft 

ft 

ft 

Z Z O 

ft 

H 

UJ 

fO 

o 

CD  >•  >“ 

w 

(SJ 

ftt 

tmt 

Z W 

X 

M 

UI 

K 

z 

«• 

> kU. 

(D 

o 

o 

ft  to  ^ 

z 

a 

z 

z 

•0 

o 

B CSJ  • 

ft 

nH 

o 

•JS 

13  -a  X 

X 

4 

T 

UI 

• 

V 

• ft  .f 

M 

Z 

l(H 

<(S 

0 ft  O 

(> 

M 

• 

X 

<VJ 

*> 

CD  to  X 

mJ 

to 

z 

ft 

« k-  O 

to 

z 

to 

UJ 

Z 

M 

>•  > u. 

ft 

UJ 

o 

X 

k X < 

k«< 

UI 

to 

< 

u 

X 

* • * 

fO 

o 

> 

o kH  ft 

ft 

z 

< 

«r 

«> 

ft 

J-  fO 

X 

<o 

ft 

B 

X • -ft 

(Si 

ft 

' X 

li 

kO 

>•  -I  X 

z 

ft 

z 

UJ 

k- 

ft 

«r  k»  o 

c 

ws 

1 * 

O 

ft 

X 

A 

# 

B UI  ft 

H 

M 

ft 

X 

ft  X 4- 

< 

k^ 

1 r 

to 

0 

O 

M 

••  B • 

Z 

Z 

w 

M 

a 

o 

(0  2 

z 

z 

< 

z 

w 

OJ  CO 

z 

-t  # (Si 

ft 

ft 

z 

ft 

ft 

o 

eo  ft  X 

ft 

UI 

x 

o 

a B 

X 

ft 

> »o  X ro  (SI 

u. 

US 

ft  (0  kft 

X 

K 

♦- 

kH 

z 

•» 

» 

* ^ ft 

Z 

z 

ft 

ft 

o 

ftp 

^ X o 

u 

ft 

mi 

h- 

Ui 

UI 

< 

ro  UJ  • 

>- 

ft 

w 

UI 

•H 

z 

z o < 

tf 

(SJ 

< 

« 

,J  k- 

X CD  k- 

>■  o z CD 

X 

X 

Ui 

O X ft 

o 

*H 

X 

z 

X CO 

«•  < 

Ui 

a •> » 

ft 

M 

ft 

M 

kM- 

k- 

O ft^ 

z 

z 

UJ 

<t 

fO 

X 

» (SI  ft 

(St 

ft 

(SI 

ft 

u 

o 

ft  to  o 

Jj 

ft 

X 

•> 

W -J  • 

z 

CD 

0 

s 

& 

ft 

O X (O 

to 

1- 

CO 

o 

UJ  CO 

«> 

a 

a 

>•  UJ  fO 

ft 

Z 0. 

z 

ft 

z 

< z o 

o 

u: 

k-  UJ 

M 

z 

• 

• a > 

X 

ft 

ft 

UJ 

ft 

o 

z u < 

a 

3 

k- 

CO 

k- 

X 

z 

<T 

(SI  • U 

ft 

I 

O f\j 

X 

X 

ft  ft  ft 

z 

V 

< 

111 

UJ 

• 

» 

^ M 

ft 

z 

k<« 

M z 

u> 

k-  to  k- 

o 

M 

H 

Z ui 

CD 

k- 

UJ 

B Z (SI 

kO 

ft 

to  a 

ft 

-I  r 

0 z S 

«» 

ft 

o 

< 

<3  X 

M 

X 

X 

• > > 

ft 

X 

ft 

V 

Z M 

z 

ft  z o 

o 

4 

(n 


O X 

K </l 

O tiJ 
U.  Ui  ^ 
X o 

M O M 
0^  U.  X 
O O UJ  hJ 
O 0^  CD  > 

a:  0. 

ft  Z Ui 
Z «X  X 

C/l  UJ  (J  k. 
M •«- 
X I 


• Z I 
CD  N*  « 

z r I 

3 » 
X Z < 
M O » 

• M C 
X V 
X z 

o 


o 

X ^ 

X z 

c ► 

o o 


C\J  CO  u.  u. 

> a V » 

V «•  UJ  ^ 

> K'  o M 
O Z Z X 
UJ  > o u. 
c/9  CJ  u.  • 

V O N to 
Z * Z CO 
O <NJ  o U. 
Z Z X * 
Z ^ X (O 
o o u. 
o o o u. 


X ft.  CVJ  O Z ^ liJ 
V o X M v •!- 

-J  »-  K Z ft.  V 
(.>  to  z < ^ ^ o 
X X 00  Z ft.  to  z 

UJ  » ftr  > ^ ^ Z 

> o v)  a X z o 

K ta  V N V • V 

z •>  Z Z Z I.  _ 

OOOOOk-Ot  ___  

ZUJZZZ(OZZXZZk*«4tl  .< 
zuixzc«jxzsczKMi>a 
ofLoooeLO>poooovzz 
ucoouu  »UUO<  vZZ 


I-  M < 4*  ' 

^ V w 

V liJ  ^ O V 
UJ  to  O CO  m 
z < X r fo 
M X O <H 

V X Sv  e > 

z z z 

>(/1000'>«C\Jt£>U9 


H O 
M M 
CO  o 
M O 
3 Z 
00, 
o 

< mi 


ti 

lD  • 

» CJ 

n 

LU  N, 
UJ  to 

a X 
CO  z 
II 

C t« 

UJ  to 
I Ui  z 
: a Of 

> CO 


X 

0. 


• O M O 

» 

ooooooooooo 


liJ  br 
..J 

z z 


^ o z 

^ fo  o 

^ » Z 

* o to 

c/)  • «-t  z M 

CO  (SI  to  UI  » o 

UJ  Z Z T*  • 

Z»  UJMBOlte 

zoXQoontp 

II  ^ k-l  tfl  lA  :-» 

to  « B > II  M o ^ 
(O  II  tOZ»M  — 

UJ  U CO  to  X X z 
Z^MMfi.MOUJ 
Z < M ^ B fcr 

(0  ft 

X z 
ft  H 


UJ  w 
3 H 
Z X 
M UJ 
X 
Z 13 
o to 
a M 


o 
z 
^ ft 

If 

CO 

X «i 
ft  < 

*i  u 


o o 


123 


INITIAL  CONDITIONS 
YN2=TER(1*NM) 


lA 


00 

K 

V 

o 


o 


T3 

to 


vO 


CJ  V 

V . 

— <VJ  * 


II  ^ 

I ^ •• 


z i/i 

(/)  II 
< <x 
M M 
M Ui 
I X 


V >•  V BL  »-  >- 


• • nr  o;  X v 

^7  O U)  UJ  UJ  *>•  (/) 

_l 

I •aooa<xn**||  C 

m O H M ••  II  M lo  ^ (j 

••OXZZIIQ^O^Q'  X-l 

oaaook»K»v  zu 


o l^  e 

(A  «jD  K 


in 

K 


124 


- 

t- 

«* 

u. 

a 

Z 

Ul 

O 

o 

UJ 

o 

ty 

o 

Ul 

Q 

f» 

u 

(A 

X 

V 

» 

o 

o 

UJ 

2 

z 

lu 

a 

t/1 

O' 

Z 

t-i 

LD 

X 

<x 

X 

» 

o 

M 

o 

% 

2 

0. 

X 

u. 

<s 

M 

r 

1 

t 

1 

1 

CO 

U-  X 

X 

X 

(/I 

M <Z 

< 

•■ 

X 

X 

CL 

tA 

(/)  M 

o 

CD 

>-  • 

L- 

1/1 

» .J 

C/1 

* 

r r a 

u 

X 

(/) 

2 M n 

a 

» 

U. 

DC  •-  U «> 

ID 

M 

(A 

> OC  ^ DC 

— 

CM  < » 4 

X 

-1 

•> 

CD  UJ  CO  UJ 

>• 

UJ 

UJ 

UJ 

UJ 

<VI 

fO 

*■0^  (A  X 

o 

o 

o 

o 

o 

•J 

(A 

-1  ->  UJ  X 

o 

o 

o 

o 

o 

» 

•> 

» 

2 

r 

2 

2 

» 

CM  — G.  <«• 

«» 

2 

CD  CM  # O 

a 

UJ 

UJ 

UJ 

Ui 

•J 

• 

• 4H  <J  O 

o 

UJ  2 

QC 

0/ 

QC 

• 

CM 

2 '-'  fO  ID 

ID 

ty  M 

M 

M 

M 

CM 

•J 

<0  Or  <t  • 

-* 

M ^ 

»- 

K 

W> 

» lU  » o 
fO  Q o H 

CO  • w 

hi  «x  w 

«.  £ m ^ 

*-*  C/1  a 
2 t-  CO  •■ 


•O  OJ  £ ■»-l 
^2  0'-' 
• M O ►- 

CM  K •■  J 

CD  o a CL 

. Q£  » < 2 

: liJ  fV  •• 

I (M  2 *2 


^ A A «» 

«0  OO  «D  qO 

M o o a o 
M (T  O'  O'  cr 


|w«^--w  V V V V 


^ Of  O'  O'  0^ 

«»  X 2 S Z 

. ^ ^ ^ X 

)«»C9M«>.C3T^(Mro  -^a 
'POOs,r>oK  • • • •N. 

>'^rr<\^’^f3nfctr*cf  * 
)CD  *■  ••rOO'UJUJUJUJX 
► • • • *0 

l_J2(M_|0'aJUJCUUI(M 
iOZ'^OQ^Q^^QiOi^ 

:<QfO<u'ii.u.u.i^a( 

:o2U.oq^mihmi-4u.< 


o * o.  &. 

uj  c/>  o o 

UJX«» 

X (Lq^cm  xcD(vjro^xcoojroorou.cDffi«»(otoc/> 

c>  covo  ^X222XMM^UO(/)C/)^(0<X^Z 

» oncT'VwriooooooooSQOOQnooo 
X cDcs  ^ ^cDCDSCDCDmmeocDeDCDCDcDODCQcaaaeDm 

'-'2  0 0 2 »»ooooooooooooooooooo 

; < Off  Ui  < 

.U.OOOXU.UCJUOUUUOUUOOUUUUOUO 

J ^ 


125 


rj 


u 

o 

CL 


OD 

to 

to 


z 

>- 

u. 


o 

oe 


3 

O 

oc 

a 


a 

O' 

1/1 

o. 


0.  «»  «»  (/)  <>  u,  < 

<x  «.«.(/)Zi-i. 

— oz— f>oujM»-’ 

fcp-oj  cviruojza  rozcozcj  Hv*<ccn^cOQ:jco- 
cno^Vcricn-)o:>->^30ir</)o«>vcnQ'ua3<<aH'^< 

0(y‘'-'Ccr‘OP*'-'r'QoQfT'^oooae'Cor3CCi 
CD  ^CDCDOS'^flOaDCQCQCDCDCO  •ffiCnfflCDeSCaCDCSCDalCDi 

ozooooooooooozoooooooooooa*-’ 


* » 
I O 
: O'  w 
I ^ ^ 

o -r 

O' 


*4  ^ ^ 
II  *-  tl 
5^  — 


UJ  UJ  « . . _ . . 

«-ar<<r<«'i<<<a<t<<o'<<<«<<<«T<r<<ujC7ujoif  oii.a.  ii 
OZUUUUOUOUUOUZOUUUUOUUUOUtfU.O'U.-JOMZ-J 


o a u*  z 

U Z Z Ui 


126 


3 

V) 


4 

a. 

4 

o 

UJ 

• 

Uv 

o 

4 

CO 

X 

4 

X 

» 

4 

Ul  z 

4 

to 

CkC  Z 

4 

X 

M » 

4 

M z 

4 

< 

IH  X 

4 

4 

l-T  X 

U-  X 

X 

4 

X to 

M < 

<r 

4 

» » 

M Z 

X 

4 

a to 

to 

M 

4 

UJ  o m 

M •- 

•* 

4 

hm  K to 

4 

fO  ► 

X Z o 

a 

4 

e X to 

z M e 

4 

UJ  o •u. 

» -J  u\ 

4 

z X ro  to 

or  M — 

CO  or 

4 

MO  ^ * 

< » X 

x < 

4 

^ o • ^ 

UJ  to  X 

CL  UJ 

4 

3 tJ  tvj  M 

« to  a 

3 tL 

4 

o < d to 

3 UJ  o 

• X 

4 

Oi  • • 

* at  • 

a • 

4 

ca  o »o  fo 

« a. 

4 

3 M z -i 

M • o 

to  a 

4 

to  » 

rH  o a 

3 43 

4 

to  * ai 

w (D  '0 

• ta 

4 

O UJ  M ^ 

or  < -* 

1-  • 

4 

z CO  z * 

Id  • X 

X O 

4 

M < -5  4* 

o u > 

M M 

4 

(OX  • z 

» X o 

• w 

4 

<0  tL  M • 

Id  < • 

N-  K 

4 

UJ  X ro 

z • ^ «I 

X -J 

4 

U Z > • z 

M (0  Ci  a 

z a 

4 

O M etc  CM  * 

kr  CD  00  o 

» • 

4 

Oi  O X <Vi 

U.  • ^ tfl 

(0  « 

4 

0.  M X ••  z 

• -J  -J 

X o 

4 

3 UJ  ^ 

UJ  Z UJ  CL 

0-  -H 

4 

a o X z X 

X O CO  Ul 

4 

M 

M O < 

k- 

m 

< 

o 

O 

M 

•L 

k-  ••  -J 

o 

CO 

Ul 

Ul  z 

0 

o a • 

• 

z 

0. 

e 

z 

M 

t> 

X 

a. 

» < o 

a. 

z 

o 

at 

M 

(\l 

• 

<r 

UJ  CL  ar  e 

«> 

M 

O' 

X 

z 

T- 

C3 

z «>  z 

o 

X 

o 

kf 

•J 

M 

«x 

M 

M 

MM  * U> 

CO 

M 

a 

M 

a. 

M 

M 

U ^ CD 

•• 

to 

M O 

I 

^ z 

to 

M 

X 

•> 

V » z 

O' 

0. 

a 

»-4  O 

O M 

z 

X 

n 

o »-  z 

td 

M 

z 

^ o 

ar  3 

M 

M 

tr 

V 

o 

M X N 

K 

N N. 

0. 

cr 

z 

X N X 

K 

UJ 

M 

a. 

o 

z 

UJ 

«> 

to 

<r  UJ  3 

o 

(O 

J 

M 

• 

K 

< 

o 

X 

I 

I 

Z fV  z 

z 

< 

a. 

to 

M 

M 

M 

M 

Ul 

>-  M O 

CL 

X > ai  <3 

a 

M 

< 

M U. 

> 13  O K M O 

0. 

X 

a. 

• . 

bm 

M 

O X 

V 

13 

V >k  V 

N 

X 

X 

^ to 

3 < 

CL 

X 

z 

m 

z z z 

z 

z 

Z 

(L 

X 

o 

M 

0. 

o 

a o c o 

o 

CO 

O O M 

a. 

ar 

3 

-J 

z 

UJ 

Z X X 

z 

X 

z 

Z 

(O 

UJ 

«« 

X 

4 

M 

M 

z 

O' 

Ul 

Ml 

M 

o 

o 

z 

M 

X 

3 

u. 

u. 

•J 

to 

A 

3 

< 

M 

o 

to 

M 

z 

Ul 

z 

M 

M 

r 

M 

Z 

o 

o 

M 

a 

M 

X 

UJ 

< 

fO  M 

z 

X 

«X 

O 'Tt 

M C3 

M 

M -J 

-1  ♦ 

Z 

fc-  + 

z 

< UJ 

u.  o 

t-  Ul 

< o 

o 

a c 

1 Ul 

z ^ 

O «d 

X 

X 

-1  z 

Z 

J z 

z 

o 


U) 


CD  UJ  C7 
I J o 


>r  *1^  Ui 
X *0  O Q. 
K 0^  O (/) 


Xftrzz^ 

0000>-00fl,u. 

OOOO  •‘OO  •*-! 


I u 

I UJ  *-<  — 

: o ^ ro 

I X (M  • 

« • o 


o u.  oe 

. Of  w 

: z 

: «t  < 

. a X r 
<x  a:  Qi: 

td  o o 
0^  u.  u. 


U.  V 

o — 
(X  ♦ 

Iij  fO 
CD  U 
Z ^ 
X OJ 
Z _l 


X X 
V 

A m 


^ < 1 
i-i  a:  ( 
ct  oc.  : 


4^09^ 

cv»  • X o 


<—  +01 


o 111  r 
QL  X I 
0.  h-  , 


I CVJ 


U. 
UJ  M 
Z H 

>-l  o 

X Id 

ac  z 

UJ 


fO 

j»_<  »^x  tvi^ 

X Z ^ (/) 

M b-  o ••  u; 

II  U.  < M UJ  h-  ^ 
Z X M Z to  X 

M q:  o II  UJ  «i>' 

_U,U-OUJZI-U- 
ZXi-»->Wu.ZZX»-* 


^ (O 
o ••  UJ 
Ui  W tr 
X (O  X 
tl  UJ  T-” 


O 

UJ  o 
UJ  Ul  J 
-I  B.  Uu 
«f  CO  * 
O V ur 
CO  -J  X 
UJ  M 
UJ  O >- 
Z X c 
11 


O UJ 

u.  o 


UJ  CO  CO 

c n ^ 

X » 


• -i  af 
I ->  n 


CL  CL 
!•  ^ 
to  to 


d UJ  c 

z r 1 


X ^ K » 
)-  O U.  i 


u. 

o o 


to  ' 

Id  n 

z « 

IM  (L 
CD  w 

(X  X 

« M 

X u. 


to 

o o 


X V M *3“ 


z cr  ry 

O UJ  UJ 

o • • 

W UJ  UJ 

X or  at 


UJ  ui 

oc  0£ 


li.  >•  u*  ^ K 


127 


MAKE  SURE  THAT  NUMBER  3F  BATA  READ  IS  DIVISIBLE  BY 
ITESTtIPHS 

95  XTEST=8.*FL0AT(ITEST/8»-FL0ATCIIESTI 
IFCABSIXTESD.LE.O.sailGO  TO  94 
ITEST=ITEST-1 


SUBROUTINE  PROFILE  74/F4  OPT  = l FTN  01710/78  ll.OO.SIt  PAGE 


128 


SUBROUTINE  PROFILE  74/71*  OPT  = l FTN  4.6»433a  01/10/78  11.00.54  POGE 


w 


Uj 

z Ui 

M O' 

M 

3 H 

O M 
Of  •• 

OD 

CM 
{/I  CM 

4- 
UJ  CM 
w *• 

M CM  O b. 

b.  » «X  Z 

CM  llJ  M 

Oj  CM  O'  CM  O'  U; 

«|  r-  Kr  CM  a.  3 

M O'  b.  Z Z 

b.  o o w a 

O ^ J ►-  -J  K 3 

CL  ^ ,J  Z 1-0 

(L0<tC<0  UJZ 
O u O O O o'  bJ 

j ro  .♦CM 

4 CM  CM  CM 

O 

O O 


129 


Ui 

o 

er 


OC 

QC 


Ui 

Z 


a 

Ui  z 

O'  z 

*-  z 


li.  X 

^ M « 

z 

a t/» 

o H-  • 

►-  » ^ 

^ ra  z o 

UJ  ^ M O 

a » ^ lA 

O fO  k.  <-<’ 
r -I  • X 

a (sj  (/)  o 

z «J  UJ  o 

< • z • 

ffi  M)  0.  ^ 
r • a 
a o a 

«f  • o i<> 

Ui  <M  ^ 

/v  Z • X 

^ a > 
» X o 
O ko  ««  * 
M X • ^ 
O • C/O  o 
M (\i  (S  O 
Of  X •>  irt 

0£  ^ Z QC 
O Z O Ui 
U » U H 
fO  • o 

Ui  z a 
z ^ 
M csj  Of  e 
l~  Z •>  o 
0 3 •h-  \fy 
« O <X>  CO  ^ 
uj  « *-4  * a* 
Of  CD  • >-  Ui 
K3«X»- 
ar  {/I  z V X. 

* UI  e 
UI  Of  T Z Z 
Z Ui  M M Of 
»-•  ♦-  » ♦“  o 

»-  -J  X X K 
3 M z Z Z 
o u o o 
Of  ^ z z 
<D  UJ  < z K 
3 ^ UJ  O O 
(/)  M or  (3  u 
u 
o 


cn  • 
I — 
a a 


a u 

U'  a 

t/l  3 


»-  X ^ 
X • o 

z -M  IT 

» CJ 

to  tn  k- 

X w Q. 

O to  3 

X M » 

(0-^0 
X o u> 
a UN  •- 
«l  V 
• O CO 
CO  • •* 
X 


«.  «•  (O 


♦ 

I 


- o : 

• >-  ' 

(O  (/)  X o t 


I «-  o CO 


<o 

2 


u w. »-  >-  X : 


X z z 

a o o 

V M H 

z CO  CO 
CO  o z z 
z r uj  Ui 
Of  z z z 

V o M •-< 

B o 


; o:  • X 

> ^ k-  a 

■ X ^ < 

> M • Of 

: u CO  II 
I k-  Z CD 
: II  w u 
I to  u.  to 
I z X 


u -J 
lij 
G 

X cvi  ►-« 

Ik  ^ CO 

rr  G Z ►- 
O'  to  <r  <i 
<r  z Of  w o 
<r  » w II  j 
• OJ  X M u. 
O ♦ l> 


CO  X 
Ui  I 

a G 
o < 
u a 

CO  » 

a 

X < 

(if  (kf 

X ^ 
a Of 
►-I  o 
a CO 


C^  II  M W k-l 
UJ  CM  •• 
cocoaok- 
X Z 3 o X 


G 

k-i 

to 


a 

o 


o — 

— ►-»  CO 

o — ' o 


J ^ CJ  X »- 


a V 


kxk>  ^»-i«-lfUi^ 
XX  I-*  ^ H G 

M ^ Z 

I O k-l  tM 


4 4-.  Ui 


.J  z 

3 (O 


X — k-  w k> 

^ ^ k»  Z 

u.  a o a c 

*-•3030 


I II 
CO  3 
4-  Z 

a a ^ 
3 *-<  ai 
••  II  ^ 
^ 3 

a a o 

3 3 0 


(\i  •»  I 
CO  »-•  • 
r I 


C7  II  < 
cvi  H I 
o a ( 
n c 3 ' 
X 3 3 • 


to 

k- 

z a 

H-*  3 
<T  3 

a — 

a X 

-ij  >“ 
K 3 
II 


^ k*  : 

3 to  : 

CJ 


••  G -J  •«. 

»-«  U.  M • W 
«-  Z k-  • 

k*  G CO 

CO  a II  II  uj 
II  o ->  ►-  K 
• U M 00  X 
M »|»  UJ 
W o 00  ^ u. 

G ••  M X M 

to 


• M 

M G 
«•  UI 
U.  O 


UJ  CJ  *■ 
• UJ  ** 

M • t 
kH  O ' 
••4  UI  ' 
to  O I 


M I U.  ^ 


^ o 


• a I- 

I O CO 

r -J  — 

• CO  » 

r N. 


G 
I G I 

I I •*» 


O Ui 

^ Q 


o X 
r U'  X 


M M a * cli 

X a 

II  <4; 


I M •-«  M CO  X 


k-  CO 

X •• 

II  — 
- X 
X >-* 
^ o 
O X 
X 3 


Ui  II  • 

3 »-•  » 

z 


O O I 

u o » 


I O M 
) X w 

>•  o 


G 

o o 


3 

o 

a 


UN 


130 


SUBROUTINE  RTRE»0  14/74  0^1=1  FTN  4.6+433B  01/10/78  11.00.54  PAGE 


f>J 


U. 

» « 

Ul  M 

a ^ 

X u 
I X 

• 


X o 

X 

» I 

«•  o 

(J  < 

& q: 


o 


</> 

z 

o 


4- 

o 


or 

Of 


u o 

o 

"J 

z 

X * 

(/> 

u. 

a 

z 

ft 

z 

o 

t/i  u 

tA 

UJ 

0. 

(/) 

ca  fi. 

* 

Z 

X 

0m 

X 

< CO  z 

< 

M 0m 

• 

•« 

X 

» 

M X • 

^ 

M 

X bH 

UJ 

X 

0m 

oc 

X 

z 

ffi  0. 

X M — 

•• 

-J 

w 

X 

•3 

•3 

X *3 

M ft' 

Z ft'  X 

• 

A/ 

Z 

C/l 

U.  U UJ 

ft. 

c/1  UJ 

UJ 

u 

X 

e 

O' 

X 

o 

z 

3 X Z 

M 0.  H 

tat 

Ul  X tar 

X 

X 

X 

M 

• 

UJ 

♦ 

•• 

• U 

M ft  Ul 

1 ^ o 

II  bJ 

Z Z X 

M 

X 

c 

UJ 

UJ  “3  O 

X 

^ O X 

a 3 X 

M ♦ II 

S 

X M •• 

M Z 

M 

9* 

3 

3 

II 

11 

Z 

li 

11 

«C  II 

^ 

M Z 

O «• 

UJ 

0m 

Z 

Z 

CJ 

0m 

0m 

•9  fO 

0m 

«k 

z ^ z 

o 

M M M 

«i>  M 

'J  ft.  n 

o 

X 

•9 

X 

X 

CM 

M 

CM  "3 

CM  3 Z 

^ U «» 

0£  ^ 

z «- 

X 

w 

X 

X 

X 

m0 

X 

X 

m0 

Ul  X X 3 

O IL  b> 

Lii  Z 

M X (kr 

a; 

X 

z 

z 

X Of 

-3 

X z 

> z X a 

o 

0.^0. 

►-  o 

z o UJ  o u 

ij 

X 

o 

e o 

X UJ  K O 

X 

UJ 

O O UJ  Ul  z 

o 

^ ^ 1- 

X u 

Of  « X 

a M 

X 

Q 

u 

u a o X 

a 

o 

X 

z a X z UJ 

ilJ  X 

j 

lA  Ik- 

e 

X z 

e 

X 

Al 

w 

z ^ 

^ bH 

fO 

Ul  o 

(M 

(M 

CM 

• CM 

o u 

X 

X 

X 

• X 

u 

u o 

u 

(J 

131 


iu 

O 

«x 

a 


* 

sO 

dr 


M T 

t-  it.  X 
V M <3 

* ►-  r 

a.  c/)  -> 

o *“  • 

fA  K a 
O ^ M O 
bj  » ^ 

O pa 

O -I  » X 
IT  ► ta  >- 
tSJ  ca  o 
UJ  o 
Z » 0^  * 
M fa  0.  -• 
tt  z » o 
a.  n o o 
^ • cn 

O M « «> 
O Z • X 
u.  n O > 

• V Q 
O rPt  « » 
UJ  X » 

X *(0  0 
M (M  CD 
U.  X *■  lA 

O'  ^ Z O’ 
O X O itJ 
U.  *0  H> 
fa  • o 
111  z c * 
z » < — 

M Pi  Of  a 

t~x  ^ 

H-  Z3  • I-  lf» 
Z O CD  CO  *-> 
M q:  M » ar 

0^  CC  ^ IJ 

ascDx*- 
iL  ca  z V ^ 
•>  UJ  o 
Ul  V X O'  X 
X iu  ^ oc. 
♦Ml-  •■  I-  C5 
H-  -J  X V K 
3 i-t  Z Z Z 

o u,  o o 

IK’  -I  Z Z 
CD  UJ  «et  z Z 
3 J liJ  O C 
ca  •-*  a:  o o 


0. 

o 


’ o ^ ">  at:  3 

: o I I uj  II 

. -J  Ml  ca  5*  "5 


3 n 
I t 

u.  u. 


(a  o 

z z 
z z 


X -1 

*-l  » 

b.  (M 
•M  M 
II  b. 
^ b. 
n z 


X * 

-I  a 
*a  M I 

11  M : 
Ml  T I 

z z : 

•a  *a  I 


zoo  Oi 

M « « U/ 

tH  Z X ^ 

o II  UJ  X UJ  (a 

• • M ^ O 3 il 

o o ca  (/)  z <• 

II  - C3  II  |l  *M  -a 
X Z Ml  0^  X kM  W 

X bj  UJ  > z q: 

O »-  o ^ O bJ 
ca  la  o ta  ca  o H 


w cy'  w 


(a  bj 

H O 

« z 


X z 
X o 

o o 


M II 

OJ  ^ — 

CM  -a  "a 

X a: 
o X bj 
e o X 


•M  3 O < 

4o! 

II 


O’  O X 
•a  e o 


O'  f_  o 
UJ  iij  z 
X z w 


Q£ 

a. 


oc 

a 


D 

O 

z 

flfi 

D 

(a 


132 


UJ 

o 

0. 


► X UJ 

X 

^ z ar 

z 

CD 

a UJ  X 

ro 

X 

U X 

X 

ro 

Z 

X M 

z 

X 

X » ► 

» 

* 

(y 

ox  u. 

X 

%o 

0. 

O > X 

4 

• 

X 

4 u.  X 

z 

o 

•■  * (O 

3 

c 

-.X  X 

z 

u 

ox  » 

A 

X 

4 U T 

a 

u. 

UI 

(0  W • X 

o 

> 

o:  X 0^  ^ 

in 

M 

UJ  CO  O X 

X 

X O 3 • 

X 

& 

UJ  4 • to 

X 

4 

Z » X to 

Q 

O 

4 A X UJ 

o 

4 

Of  o • oc 

4 ^ Z 0. 

0^ 

0.  w UJ  * 

o 

o 

X s o 

o 

u. 

a o o M n 

m 

u<  fi  to  cy  4 

CO 

X 4 M Uf  • 

X 

0£ 

4 * X O 

X 

O 

X «»  CM  * ^ 

o 

X 

c^  o A no  4 

m 

O 

O -t  4 X •• 

A 

UJ 

to  — z z to 

o 

to 

to  Z * UJ  13 

-o 

4 X X X • 

in 

o 

Q O J 

X 

U 4 4 (M  Z 

O' 

z 

O ► CD  X O 

UJ 

X 

X Z Ckf  o 

X 

to  o J m 

o 

UJ 

UJ  (0  CO  X e 

# 

q; 

3 0 0 » 4 

X 

X 

-J  4 4 X cy 

o 

II 

X 

4 » Z X » 

o 

X 

X 

XX  » O X 

tn 

a. 

z 

UJ 

o — •*  o 

«i' 

o 

M 

a 

UI  O o ^ 

(V 

Of 

X 

Z 4 ^ w X 

u> 

Q. 

> 

X V •*-  M X 

X 

4 

X 

Z -1  O X V 

V 

o 

IV  O to  »•  Ui 

a 

X 

UJ 

UJ  z z ro  O' 

z 

K 

' z 

o 

X o •>  X X 

QC 

S 

M 

X 

UJ  4 «»  X X 

IS 

X 

ht 

O K o X V 

V 

3 

UI 

z -a-  * z 

z 

C 

z 

O O <v  CM  O 

O CO 

Of 

M 

Z Z X tH  z 

z z 

03 

X 

< Z X M z 

z ty 

3 

3 

0 0X0 

O X 

X 

CO 

O 

^04  » O 

u » 

z 

ct 

UJ  X CM 

X 

X 

CO 

Q 

0^ 

3 

O 

Q. 

(0 

z 

< 

u 

C3 

o 

a 


3 

cn 


w 

X 


UJ 

o 


a 

o 


w UJ 


< 0^ 

Z 

««'  • X 

M QiT 
< N.  O 
O I—  U.  >H 

^ j a o 

<M  ♦ li.  a ^ «* 

«<l  « (/)  V « X z 

♦ (/J  O W — C3 

in  O « CJ  z M ^ 

o <*  z in  M ftr  II 

« Z » « (/D  M 

O'  Z II  !\J  o II  Z 
••  M btf  M ••  fcr  O O 

(/)  o (M  H'  e ^ 

o o « e a t/)  CO 

«<oD«oe^xo 
zzzz«x«xoa 


U.  CO 

••  li 


X b* 

^ (O 

o o 

< < 


-I  o 
« t-l  M ui  J 

' w w o U. 

: I I X * 

• ^ ^ V u 

I o a — V UI 

: < < M > o 

' ta*  w w X 

• CO  Z K X I 

o « CO  CH  X 

U X O b.  X 

M II  < »-l  M U) 

» M «.  3 

M M a <*  M z 

w M w *H 

X X or  X k. 

u x ••  o a z 

a o X to  o o 

«x  <r  X z «K  u 


a 

IH  O 


^ < 
I o 


’ ^ UJ  . 


oj  z i-i  »-i  »-i 

^ Q t ^ ^ 

Qt  < M X X X 

4 Z X X X «i« 

IJ  «>  X 3 O O M 

«,«  o O X 


r to  : 


: CO  1 


u.  (D  u.  X o < o 
UI  Z I to  o X •! 
^ II  H IIMII#. 

X o 

z X M X X 4 

o o (k* 

CSJ  X X X X I 
(O  X CO  o X II 
nf  o c *3  o o X 
O O < < «X  4 X 


X UI 

X o 
X 
X I 
X X 
U X 
X I 

I M 

II  •i* 


«■»  X 
O Q 
CO  n 


UJ 

s 

z z 

X Of 

X s 
Z X o 
o UJ  z 
O Z UJ 


133 


subroutine  otsts 
dynamic  simulation  subroutine 


o 

<r 


4- 

sO 


3 

o 

Q£ 

m 

3 

CO 


CO 


z ► <vi  m 


03 


to  J 

• V 
U » 

M 7 M 

* -0  X 

« T • A 
7 K*  r> 

V z ^ • 

C 3 • M 

a cv<  « I 
uj  ^ a. 

^ z -cr  b>  < 
< z cr  ^ 
to  *<11; 
UJ  cn:  ^ ^ z 
»-  lu  z •>  h- 
« *00 
z ^ M CkT  * 
•“•nzZNtf 

a;  CL  CD  ^ UJ 
O ► t-^  Z I 
O Of  ••ZK 

o lu  CD  V V 

►-  Z z UJ 


•■  CO  ^ 

a >■  V CO 

>.  • to  •* 

fU  ^ Ul  CO 

O CL  » U. 

• CO  fO  CO 

CD  >-  >-  • 

>•  • CO  ^ 

CM  K>  U.  Ki 

e>  0.  • to 

• to  <M  * 

CD  > > *0 

» u.  u 

O CM  • » 

• a.  ^ <\i 

CD  to  T «J 

> >>  U.  • 

• • • ^ 

^ ^ fo  z 

>.  _J  T IV  • 

Q UJ  u.  ^ ro 

O » Z Z 


, >•  c <f  X a 

O 4 > » *H 


z 

O UJ 


O > 
O >- 
< U. 

X 
a X 

■s-  u. 


» CM 

m u 

> UJ 

♦ o 

<\i  » 

> ^ 

o z 

• ^ 

CM  CM 

>>  o 


Z ^ X Z «r 

»-*  JX  M o ^ 

H*  0.  • M (/: 

3 X X.  V 

o z r z z 

CkT  O o O 

a -I  z z 

3 to  < Z Z 

to  >-  u o o 

3 ar  o o 

to  a 

M U- 

X 3 
K O 


3 ►-  ^ 


O - 
UJ  > 
CO  CM 

X o 
r • 

O CM 

^ z z 

Z Z 

>-  O M 
o u o 


U.  CM 
•>  O'  CD 
ro  Z 

> z » 
• u.  » z 

» ^ M 

CM  lO  » 

> X X 

u,  u.  z 

X X 

UJ  ^ u 

o »o  o 
a X X 

O U.  Ul 

u.  • > 
X to  X 
Z X z 
o u.  o 
r » z 
Z to  K 
QUO 

o u.  o 


> X or  CL  • *- 

CM  •>  < «X  X 

o CO  Ul  UJ  >• 

• to  DC  O'  3 

n UJ  3 X o 

> cx  » » * 

CM  e. 

CD  • CM  CD  O 

4 U c9  -a 

z o W UN  tfs 

CO  «r  a * «« 

• • UJ  o X 

»o  o o »-t  >• 

3 V • ••  rt 

CV  « IJ  X * 

• • Z -J  — 

CO  to  t-<  a o 

Z CD  M • O 

3 • U-  • U 

^ * o «- 

ro  Z 1x1  <-i  (V 

X C Z UJ 

• U t-«  ♦“  ►- 

M • »-  J O 

a a o CL  • 

a.  ar  «x  • z -i. 

< • DC  Ul  • O 

lO  «M  • T X 

k«  Z ^ K 

CO  3 CD  H J 

4 • X a 

a.  CM  ^ o z 

o ^ ^ M X 

W X X 


¥•  O' 

to  to  lD 

10  3 
CL  «r  » 


a.  o » 

^ z 
CO  «■'  UJ 
3 X X ^ 

• a i-i 

^ o to  ar 

X TO  UJ 

M • 4 K 

• CM  » 

W-  O »o 

X ^ < M 

z z at 


X X 

CM  CM 

• • 3 K> 

^J  ^J  ^ ^ 

3 to  XX 

» • 

CD  CO  • • 

z z 


CM  CM  CM 


X X 

CM  CU 

PM  CM 


O PM  CM  CM 

IK^  *0  3 to  K> 

Z » » ft  CM  » » 

O AJ  3 ^ CD  CD 

Z Z z X Z 7 


z a 
a *s 

X ft 
ft  K 

to  C 
I to 
a o 


a tn 
T *- 

M at 


CO  : 


u* « 


CO  ►- 

X o 

a o 

<i 
X X 
UJ  -I 


o • 

«1  CM 

OD  ft>« 

z a 

ft  UJ 

to  I- 

O »T-> 
< z o 
Z V o 


u • 


U»  3 3 3 3 3 K: 

^ ^ ^ ^ _J  -J  -J 


U’  P>J  PM  CM  PM  CM  CM 


CM 


3 


: X X 
I 1*7  (M 
. J «| 


XXX 

CM  3 (M 


a 

►-  K 

z z z z 

Z M Z «* 


a ft  » ft  ft  ft  ft 


CM  O 

_>  u. 


I CO  a z a 
ft  a M >■  X 
o to  P-  O X 

o X X X X 

ft  z z z z 
o o o o o 
U)  z z z z 
UJ  z z z z 
a o o o o 
to  CJ  o u o 


^ ft 
a o 

ft  X 

a z 

t-  c 
CO  z 

^ z 
a o 

ft  u 


X o 

a < 

X X 

z z 
to  o o 

z z z 
a z z 
> o o 

ft  o o 


^ ft4  a 

W ^ 3 

O X CD 
(O  ft  1-4 
Z 3 
tR-H  Z 
— MO 
O X M 
^ *^10 
R-'  PM  Z 
M -4  UJ 


• CM  ^ 
» X >* 
) CM  P*J 


: X I 

z z 


^ CM  ^ 
XXX 
CM  CM  tVj 
• • • 

CM  CM  CM 

3 to  1*7 

ft  ft  ft 

X X z 
z z z 


^ CD 
X ft- 

CM 


M to 

a M 

O UJ 


M tD  II  O 

a M a M 

l9  Ui  (0  Ul 

3 Z 3 X 


ft4  ft4  M CM  CM  3 3 


a 

o 


X Ul  I 
a I 


y y 

X X 
3 

a a c 


' ft  C?  3 O O -ft 

I UJ  UJ  UJ  UJ  « a I 

» • • • • a uJ  I 


' w w w M X : 


Ui  «*  II  M «•  M II 
ZfM3^(Mft7.#Mrij3^ft-J0 

»-4aaaftaa*^^^kft 

zp-»-p“toiotoatototooo 

axzz]fc^>->-M>>xjfe-oo 


aaaaaa  v* 

OOOOOO  •ujo 
«<««r<«ccMoa 
333333»»  II  «r 
II 

aaaaaaujaftC 

ft4Mft4MMML><> 


O X 

< >• 


< X 

o to 


UN  I 
m > 
•ft  UJ 

sv 

ft*  :► 

a Ul 

M U 


^ ft 
« > 


O UJ 

o u 


a 

o 

3 

I O «• 

> « a 
I a a 

ft  H 

>•  4> 
UJ  o 

u a 

N < 
O II 

-z  ^ 


134 


s •*  • 

^ <\i  <n 

X >-  X 

H * <\t 


• ► >- 
X X eu 
’H  ► -H  D 


■t  < 

r £ 

or  ft: 

o o (\J 


er*  o ► 
ui  » (Vi 
»n  ^ e 
(Ni  u V • 

•H  uj  ««  (/n 

•>  o X m 

PO  • ^ U. 

• CM  » » 

»0  ^ M I/) 

U(  • u. 

• o u. 

X ««  tu  «. 

C9  «—  o 

• ^ a «-4 

o «-<  o 

J •-<  CT'  • <T» 

• X » 

:x^'aK 
. w X ^ X 
. z w * 

w 

C « Ui  < UJ  ^ ^ 

• z fc-  r k-  « 4 
j O'  M •-»  -n  X « 
: o O'  o or  ••  t« 

> u.  z u.  z n X I 


O J b • • 
Or  < X (b  fvi 
O.  D W V N. 

Z ^ ^ 
-i  « o > z 


t U I-I  3 z 

^ b*  111 

-i  tO  « UJ  u 

< K I X 

O CO  ^ Ul 

O O o L>  o o 


C Z U.  b. 

.. 

X 

••  X 

e >>4  ^ Q,  (M 

z 

O *H  li  11 

oc 

II  «4 

Of  a. 

< •!  II  II 

b4 

Of  W «b  fb 

4 

to  ♦ 

< X 

V)  mi  lA 

K 

b»  UJ  b4  ^ ^ 

bJ 

z >-•  bJ 

Z l«  z ^ z 

z 

(A  to  O ^ T 

Of 

t.  M 

or  u. 

b^  3 0.  < a 

o 

b4  3 Z b.  u. 

^ X M 

X M 

^ U X 

o 

b>  ui 

1^  z 

« UJ  u 

^ bJ 

o 

K X X 

t 

UJ  bJ 

N. 

CA  ^ UI 

Of  tB 

<«4  bJ 

. K Of 

C » O «4 

! M 4 ^ 4 UJ  U\ 

l^*  H C/)  »-  3 

- X 0.  X z 

e b-  M e M o 

> M M ^ M II  b-  K 

I b ^ ^ ^ Z 

> U.  X c X X O o 

- •-•  bH  o -»  ^0  C3  e 


136 


137 


UJ 

o 

<x 

a 


fc- 

fiL 


O 

Of 


o 

« 

to 

3 

W -i 
Z 

O — 

M CVJ 


Of  «« 

O <VJ  ► 


-j  * r 

< «*  M 

u evf  ^ 

M H u. 

Of  ^ * 

e Ui 

r V r 

3 ► »-< 

r ^ K 

CM  O 

« w * 

k»  o-*  ui 


U.  ^ ^ 

M ^ K *»'  >■ 

O UJ  CM  >■  O > O 

^ o s.  ►-  « » 

tt;z«^ujrcM  <->  <MZ<M 

3>H<C3«4  «ii^MX 
UJOr  <£Z**>—  ♦-<'-*3**--»3*l 

z z>->-cM^M«->uja^^«0'H’ 

•-<aoc^a>«.~»'>oruj  m ► — uj  ii 

trUjMXVU;M>-o* 

30(/)ZZZ  IIIIZ 

ozzoo»-«c9«»<»nf*  o 

CD  z**o  '-'-'Mrj  '-'Mj 

sxwoo  II  ©►a'-’M<0«-i«i»<0 
COKOUUZaC/1>>V»«UO>>VUO 


o 

u. 


»<»«««»  3 


CM  w « M 


fi  o or 
zoo. 


» ^ 3 

» e >>  Z 
I >•  w 
M II  »> 
I CM  — Z 
' w O 
a w o 


Ul 

C/^  z 

d; 

-»  3 
^ fa-  o 
< Ul  z 
u UJ 


138 


u 

o 

4 

a 


3 

o 

s 

3 

CO 


a » 

c <VJ 

^ 3 

(/)  >- 

*■  » 

ro  ^a 


►H  t/)  •■ 

X >•  rvj 
a » X 
(VJ  fVJ  u. 

o h-  ► 

••  (O  ^ 

0 >-  V 

>.  •■  t/1 

w ^ u. 
a a •* 
• CO  ro 

01  >■  > 
> • CO 

cj  ro  u. 

o a •■ 


a*  k- 
4 • 

UJ  CO 

Qt  to 
“>  UJ 


3 

O 


* Li  * 
» O 

CO  ^i.  LS 
U.  O 
to  O 4 
I 4>  k* 
ro  -t  m 
-I  »o  X 

-»  to  K 

— C/l 

1*  •-'  » 

(O  I M 
IJ.  ^ 
U.  4 w 
4>  k-  » 
•>  iLl  tVJ 
4 X X 
^ U. 

lU  to  I 

X » «» 

td  4 


X I 

a a 

CO  to 


* » 

>0 


» > 

• to  CM 

• 

• 0£ 

♦ 

to 

o 

X 

X X 

kO  O 

CD  >>  V 

w 

•*.  a 

z 

» 

» 

♦ 

a 

to  to 

z * 

> • a 

a 

CM  • 

o 

• 

» 

4 

X 

a a 

3 o 

O (M  » 

• 

3 o 

• 

o 

X 

X 

* 1 

♦ > 

• a 3“ 

CM 

m 

to 

» 

« 

o 

a 

to 

M 

<M 

ro  CO  X 

a 

O'  4 

a 

♦ 

* 

* 

X 

» 

X 

Z M 

X X a 

a • 

» 

» 

X 

o 

a 

X 3 I 

• • • 

3- 

C o 

a: 

» 

to  » 

w 

u 

o 

u 

a * a 

j-  ^ ro 

^ X 

♦ X 

o 

» 

z * 

CO 

» 

t 

» 

^ m a 

X a X 

#v 

•> 

III  4 

a 

o * 

4 

4 

CO  ^ » 

o a a 

k- 

fO 

r 

«• 

M « 

# 

to 

X 

» 

► «.  w 

•*  Q • X X 

M to 

a 

* 

k-  » 

4 

m 

a> 

ro 

(O  CM  X 

4 • CM 

• 

• 

k-  CD 

X 

» 

4 * 

br 

to 

X 

X 

X a 

X kO  X kO 

CM 

a * 

a 

3 » 

a 

X 

to 

10 

a 

u. 


u. 

I 

K) 


u • 

M ^ < 

r T » 
4 » i 
z fo  : 


3 I-'  ffi  I 

or  M : 
lu  a * I 
a o 3 ' 


»■  -1  a -V  r 


i>  a » a : 
► X CO 


to 


to 


a k>  X 
3 4 X 
O 3 

a ^ 
(D  a 4 
3 a 
CO  a a 


X X 

o a 
u o 


» rvj  a 

ro  j P. 

V tir  w 


>■  jr  •» 

Cl  X tNJ 

•«  X > 

<\j  cu  a 

V 3 s 

V » a 

» kn  tS 
O X X 

a > O 
to  rj  a 
s.  o N. 
z • z 

O M O 

XXX 
Z X 

C M O 

u o o 


M • 

a ffl 

^ X a 

X * 4 

•.  X C-i 
4-  M kp 
fO  • CO 
>-  X • 

a X a 

• VO 
^ a h- 
r*’  o to 
XXX 

a a •> 

• > ts 
to  V o 
CO  z • 
a o a ( 

X a : 

CO  X a 

a o 0 
a c3  CO 


M CJ 
H-  • 
Q O 

• 4 

a a 
X • 
M W 
l~  G3 

V •* 
O H 
M a 
X N. 
4 Itl 

z a 

>►  M 

O k- 

V V 


o 

a o 

z a 

M O 

o o 

a 4 
_i  a 


X x 

o o 
u u 


to  z 
a U.I 
z z 

♦H  O 
N-  a 
3 x 

c o 

KO 


3 a — 


a a 
a a 


to 

a 


CO  -J 

u 

a M 
a X 
4 a 
o > 


^ K 
M »H 

a a 


a 

to  t3 

3 O 

to  OQ 


4 4 
O O 


w I W X 

CO  Q.  z a 

O 

0 to  CO  z 
H O 

4 O 4 CO 
k'  11  W M 
a M a M 

1 X X X 

k-  a w-  a 

O U to  CO 


Cf  « 

a » 


to  * 

> * 

to  • 


U»  f)  Z k 

[ a o a o a { 
a o td  a a c 
) o a r a CO  . 
I z a 4 > M i 
I 4 > o : 

) X T 

’ I I o o a I 

t o o ♦-  k-  > ! 

: |v  K M 4 • 

: i>4  »-•  a a UJ  i 

» a a X : 

a a 

a a k-«  M a 
a a c9  o a 
3 3 0 0 3 
X X 3 CD  X 


Z M 

a u 
a o 


a o I- 
a o z M 
to  a a (3 
k-i  a a o 
o > to  a 


1 1 _ 

k-'  a » 
o ro 
* •i'  a 
ro  * >— 
a CO  » 

CD  f<J 

» a >• 


V 4 ^ 


w '■'ax 
a 4 » X a 
aao>k-Xh»» 
aaaa  a»Hco<- 
aaxxaaxv*  4 
xx44aais  •o»- 
44  XXOCMOa- 
aa44w  II  ♦z* 
k>Kaa  M«»4k>-' 
azzooaix«-<Mk-toi 


CO  X X CM  X 

Mk-4a  xiMXfov 

xvoovxvx  • 

a ■r* 

'OCMCDO«40pHO|| 

> il  II  CO  ••  It  II  II  ••  •• 


aZZOOaQ£«H(Mk-tOI-rkOjk^U\tOIS.«Crv4 
aoooo44»i»'^’'a*  a-»-«»kO'-’»i'*-'*«'r''- 
XQ^Q^MkHaao^QfXkoxafQfxfkfQfCkfackrftf 
xaaxx'VQCaUJb’airaiaaaaaujauj 
««»«»aBu  I (300  • oooaoQ 

_ ..  . <HacM  ^ 


OOOOOUOUO(30C5UC30000 


139 


DETERHINE  VALUES  OF  VARIOUS  ACCELERATIONS 
OZTHETA  = DERI21 


UJ 

O 

<r 

0. 


D 

O 

Qi 


UJ  Qi 
O UJ 

M a 


«Vl 

Q O 


> Q 
r*5  * 
01  »-• 


UJ  *-• 
X I 

^ a 

OJ  O 

c » 
♦ 

« X 
^ fi. 
UJ  o 

X • 

^ <M 

ts  T 

♦ 

< o 

UJ  » 
X eo 
1A 

o ♦ 
» «*. 

< < 

UJ  UJ 


I i/i 

^ * 

«.  o (M  ^ 

« •-<  _J  UJ 
♦ X 

aracto»- 

UJ  UJ  UJ  >>  o 

o a a «sj  » 

II  M M a < z 

eg  ^ ^ 

X X X X UJ  a 

> >•  > > X b-  c 

(\J  (M  IM  (VI  W UJ  Z 

o a a Q o q:  UJ 


140 


«1 


* 

a 

a 

• 

o 

U) 

CO  -i- 

<;r  X 

a 

r* 

J X 

K J*  to 

CO 

X 

LU  V 

U b. 

«. 

o • 

I 00 

a 

X 

o ro 

»-  >-  K 

(O 

cy- 

z 

m 

r X 

CM  » X 

X 

t-it 

•, 

fO 

>■ 

O K)  oo 

•> 

h- 

X 

fO 

bJ  O 

li. 

<t 

M 

X 

ly 

to  •* 

•> 

» 

4 

►-<  >-  w 

X >•  CM 

CO 

a 

X 

4JD 

^ X 

a » X 

(O 

t-t 

< 

• 

•>  >> 

CM  CM  u_ 

4 

w 

z 

Lij  a.  » 

a ^ » 

to 

to 

-ft 

l.iJ  O CM 

• to  ^ 

CD 

4 

Z 

cy  h>  X 

L9  X X 

to 

X 00  V 

>-  ♦>  f/1 

4 

7 r 

o 

b. 

^ o 

CM  J-  U. 

CO 

r M 

o 

» * 

o a •> 

a 

4 _| 

to 

»-  ro  CM 

•*  CO  w 

to 

a 4- 

%■" 

to 

z 

00  -J  X 

cn  V >- 

< 4 

X 

T 

o 

ft/  • V 

V ♦ CO 

a to 

> 

a 

4H 

M CM 

CM  CO  a 

CO 

a to 

o 

-0 

H- 

b.  «J  o 

o a 4 

to 

1 a 

o 

4 

«< 

••  >■ 

4 to  CM 

4 

4 a 

4 

a 

> 

bj  fO  O 

CD  >-  > 

CO 

a 

A 

K 

a 

X z *• 

> 4 a 

CM  4 

o 

to 

a 

w -1  o 

O CM  • 

44  O 

o 

1 

a 

•>  V 

4 a ^ 

CM 

*■»  -ft 

(A 

X N • 

« to  T 

a < 

44 

K 

a 

O Z M 

X > a 

4 

X 

X 

o 

b.  -5  X 

4 4 » 

4- 

o o V 

M 

# A 

^ ^ 

7 

4 y' 

a 

4 

a 

fA  40  O 

V _l  X 

O' 

4 

111  «r 

«> 

h. 

to 

UJ  X » 

o ui  a 

ro 

X 4 

<Mi 

X 

z 

o • •H 

4 0 4 

X 

z 

M to 

a 

z 

< 

a CM  z 

4"  4 CM 

4 

4 

y-  CD 

-3 

4 

X 

o X a 

>•  ro  T 

40 

CM 

a 4 

tA 

CO 

o 

u.  • • 

•“  -i  a 

ft:  z 

4 a 

44 

X 

fo  a 4 

a 

4 

a z 

a 

a 

a 

iij  X *- 

>-  a ^ 

X CD 

a:  o 

a 

X 

o 

or  b 

C'  4 5». 

4 

M 

M o 

4 

o 


VH  * 

♦-  I- 

t/)  5^ 


V I-* 


M ro  I 
: »- 

m Cl 


fr 


• CM  u.  CJ  ^ • 
.j  * a;  CD  o Q 
KiLiro^zo.  •< 


y: 

* V 

«.  o 

M I 

>•  >• 

0 o 

1 —■ 

M * 

>-  K, 

O CD 

» 7 

^ M 

CD  ^ 
W ^ 4-  K 
X X ^ » 
> >>  ^ ^ 
— O Q 4 

<M  I t X •ii' 

X W ^ M * 

V >-  X -J  * 

O O Q I-  CM 

I w >-»  •*  M 

CM  4 * CM  W 

. c>  c» 


.J 

»!- 

4 

>— 

X 


ti 

a 

7 «r 

Al 

a 4 

X 

(0 

7 

4 

a 

to 

to 

f/> 

o 

CO  CO 

1- 

u 

< 

4 a 

a 

4 ^ 

a 

a 

a 

a 

lA 

CO 

a 

T 

X 

X 

CM  fO 

4 

* 4 

a 

a 

a 

cc  a 

a 

X 

CM  M 

4 

to 

a 

GQ 

44 

X 

z 

0. 

CL 

a 

CO  X X 

X 

* 

CO  > 

o 

z 

o 

a I 

4 

a 

> a 

X 

4 

V 

• 

O' 

a 

4 

X 

X 

X 

a a a 

a 

a 

a a 

M 

t3 

4 a 

CM 

CM 

a b. 

z 

a 

CJ 

a 

a 

a 

X 

X 

w 

w 

w 

O 1 1 

1 

CD 

a a 

o 

a 

CD 

> 

O 

N.  4 

V 

o 

a 

X 

a 

V 

Z 

X 

z 

X 

X 

X 

a CM 

4 

4- 

o o 

a 

< 

z a 

N, 

4 

a Jt 

_J 

a 

r 

N 

V 

a 

o 

* 

♦ 

a 

a 

a 

o >-  a 

a 

a a 

jf 

o 

4 a 

■> 

(3  M 

o to 

4 

bl 

a 

CO 

a 

bi 

<o 

to 

to 

o 

CD 

n 

a * 4- 

X 

4 4 

K 

a 

X <r 

(J 

X 

y X 

I 

I 

Z 

Of 

z 

4 

a 

cy 

X 

T 

X 

» 

* 

» 

CM  ro 

4 

4 

a a 

K 

z 

o 

a a 

a 

a 

o a 

a 

4 

a 

a 

X 

I 

4 

a a 

CL 

0. 

a 

a 

o 

a a a 

a 

CM 

X X 

.4' 

a 

a 

4 to 

to 

CM 

a 4 

a 

td 

a 

a 

tD 

a 

a 

a 

X X 

X 

a 

a 

a 

cy  (o  to  to 

a 

w ^ 

K 

a 

X N. 

V 

C 

V to 

a 

ts 

V 

V 

K 

N 

a 

a 

a 

a 

M a a 

a 

a 

4 4 

r> 

a 

X z 

z 

4 

Z CD 

z 

4 

z 

z 

Z 

Z 

U X 

O' 

O' 

cy 

Q 

a 

a 

a M II 

•• 

o 

• • 

c 

z 

o 

o 

CM 

o a 

o 

e o 

o 

o to  o 

a o 

a 

a 

a to 

lO 

CO 

CM  ro 

4 

a 

a CM 

X 

M 

J r 

z 

X 

Z 4 

z 

bJ 

z 

z 

X 

z 

Z 

4 

a 

a 

a 

II 

It 

M 

a a a 

a 

II 

II  II 

CD 

a 

4 Z 

X 

z 

a E CO  z 

a c 

z 

z 

X z 

z 

a 

II 

II 

•» 

CM 

to 

4- 

4 a a a 

a 

90  4 

O 

o a c > 

o 

CM 

o a 

o 

ft 

o 

o 

o 

a 

o 

a 

CM  r> 

a 

a 

a 

o o o o 

C3 

o o 

L> 

CO 

o 

a cj  o o cft 

u a 

a 

(O 

o 

o u 

4 

iJ 

4 

O' 

a 

a 

a 

e 

o 

p 

M a a 

a 

a 

a a 

a 

a 

a 

a 

a 

a 

X a 

a 

Z 

CO 

X o 

cy 

M 

o 

a 

z 

a 

o 

a 

to 

a 

3 

> 

L9 

• • * 

CM  fO  J- 

U.  b.  ij.  CM  ro 

UJ  UJ  •-»  LU  «•  O U 

o O u.  o U.  b.  b. 

I-  M t-  w I,  n 
»»  *i<'  H CM  W 

b.  b.  1/1  b.  CO  > > 
M M b-  M tr-  b.  b. 


U CM 
b.  O 

••  u. 
^ — 
> b. 
b.  »-4 


LS  U 

b,  u. 

b.  b. 


OC  0. 

o x: 

b.  »- 


Z CM 

o z 

M b. 


141 


FHJ=-FI/3*3YX  (KPHS*-NH» 
FHi»1-FV4*DYXIKPHS) 


.ua&ojTrNE  poi'^rc  74/74  o?tii  ftk  oi/ia/73  pagf 


142 


(J) 

a. 


I o '3  o ar  : 

: r M . 

: X K-  ' 

o o o * 

• CJ  c O ter  : 


a ^ I ^ < 

OH*  X ft  «-i 

f v<Lr»  ^ Vi : 

X^XUJS  ♦0«' 

^a(kr>'KZ  'T*  *t  lij^x  0^ 

X^U40HM(/)M  0(k:b.U.0M] 

V M It 

ilQ;MX^ZQ^«*-J  UMCVi  ( 

X0J->*UJ0HI^«0  :»XOUJ» 
>^>D^Ute>k’UOl!7U.U.c9<J» 
z a . 

K>  -J  O ^ O ! 

-J  « UN  M ! 


143 


SUBROUTINE  ADORG  T4/74  OPT  = l FTN  4.6  + 4I30  01/10/78  11.0a,'=4  PAGE 


tVJ 


X 

X 

U bJ  UJ 
• 3 3 


X O O UJ  z 
u.  o o o;  UJ 


144 


OJ 

O 

< 

0. 


Q£ 

er 

3 


<t  •* 
UJ  3 

r 


T*  U. 

» C/1 

I • K *.  X 


liJ 

UJ 


2r  o 
liJ  CO 

•>  c 

M -j  o: 
O or  i-U  o 
ijj  a <J 

h-  o a 
z r 

U UJ  o 
Of  U Ul  z 

M Q£ 


z u.  > 

O O W C 
»-  « 

(/I  H O. 
UJ  X < 
o o o . 

(y-  M 1 
O U)  ( 
U.  X (/)  < 
M : 

cn  e 

UJ  z o I 
z < UJ  : 

M (/)  » 
X UJ  =9  I 
Z -I  ( 
UJ  »-l  -I  • 
M u.  UJ  ( 
a.  UJ  o o • 


OK'  T 
D O 

o >- 
O k- 
«*  O 
CO 
Ul 

O'  UJ 

o I 


lu  a 

Or  Z 
< < 

lA  iA 
Z D 

o o 

w 3 


UJ  Z 
Ck'  O 


UJ  a 
I UN 
O 

t/l  X 


< ««  «>  C/N 


• O' 

I o 


• ^ 

o cu 
-I 

o • 

X -> 


(/I  < 
< in 


a m 

< M 

Z OC 

• J 

bC 


w ^ a o 


C/l  O' 
UJ  UJ 

Ul  ►- 

«x  o 

Ui  *■ 


: €i  «i  CM 


C9  ffi  t-rf  UJ 
►H  J X 
O'  O J- 


- z » 

\ • Oj 
> U)  »-  I 
I o •• 

: < 3 

• z X 

^ a 


^ -J 
< 

Z «•  «» 
U.  o O 
••  J-  J- 

I er  4« 


X a . I-  < 


w cna  — 

•■QeaJ^X  a ^ CL 


<x 

O O'  o 

o 

X 

ir. 

UJ 

U) 

u 

H- 

J- 

X 

a ^ 

< 

X 

V 

CM 

< » 

M 

3 

o 

a X 

:jJ 

a 

o 

s- 

O 

k- 

V 

W 

M 

s. 

O 

<x 

« 

o 

-J  4> 

«i^ 

« 

UJ 

a 

•oJ 

X 

z 

<c 

<J 

V 

-1 

o 

V 

-1 

^ a. 

-f 

Q 

or 

<r 

3 » 

o 

CM 

a.  UI  UJ 

tv 

u. 

o 

o 

z 

a 

lit 

cr 

t/1 

• 

o 

<A  u 

UI 

k4 

z 

• 

O * 

t/N 

• 

IS- 

UJ 

«a  X (V 

X 

C/3 

X 

c.? 

M 

z 

z 

O' 

X 

z 

m 

I 

X 

Q 

^ 00 

<x 

*■ 

o 

• 

» 

• 

-J  » 

Z 

a 

N. 

z 

O K M 

1- 

Ui  Z 

UI 

»- 

Z 

o 

nr 

M 

o 

•t 

UJ 

*■  Z 

«•>  o 

•f4 

II 

o 

o 

o 

M 

• 

< * 

* 

M 

<I  H> 

Ui 

I o 

I 

o 

z 

M 

M 

o 

S- 

«< 

t-t 

> 

o o 

O 

(L 

11 

•• 

O 

o 

• •* 

II 

• 

O 

o 

a 

o » 

a 

• 

z 

U- 

-J 

k-  M 

U' 

z 

o -> 

»- 

>1. 

V 

a 

V 

v 

a t-i 

o 

^ trv  o 

3 

X 

11 

II 

O '«» 

^ o 

•• 

II 

3 

» 

LS 

o 

3 

UJ  C UJ 

o tr 

UI 

UJ 

< 

z 

z 

z 

z 

•>  (/I 

• II 

II  ^ 

II 

z 

o ♦ 

3 O 

o 

Z X 

<I 

< 

-J 

UI 

-J 

o 

tn 

o 

o 

w 

CM 

o 

o z 

.J  N-l  V 

O 

o 

< 

y 

>«» 

a 

M 

M * 

M 

z 

O' 

»-i  z k- 

UJ 

-1 

o O' 

CD 

CM 

3 

X 

X 

X 

X 

ktf 

x 

UJ  UJ 

Of 

Ul  II  o 

z 

CM 

or 

M 

-J 

X 

Z » 

O 

3 

CO 

♦-  o 

u 

M 

M UJ 

II 

X 

z 

a X 

X 

X 

UJ  X 

-1 

a CM  < 

M 

_l 

-J 

CM 

— X 

UJ 

< 

w 

k-r 

a 

z 

K » 

rA  O 

3 

3 M • 

•>4 

< 

y o 

r 

o o 

r 

M 

o 

V 

o 

o 

O 

> 

o 

a k-* 

UJ 

o M CD 

X 

O UJ 

Ul 

J 

_J 

-J  a 

r 

o a 

V CD 

o 

a » 

II  z a 

U) 

O 1-  UJ 

o 

> 

UJ  M 

UI 

< 

UJ 

U) 

o 

► 

o 

o 

< 

o 

CO  a 

o 

» a I-* 

« 

O 

s o < 

«a 

4 < 

e 

M a 

a 

a 

u 

X * 

M 3 

M 

hr 

O'  M z 

X 

T tA 

I 

u- 

X 

w*  CM 

CM 

UJ 

o » 

a. 

CD  C/) 

13 

13  Z 

f 3 

z 

<A 

C3 

UJ 

J 

o 

liJ  » 

Ki 

< 

3 0-1 

t/l 

(O  O 

U5 

o 

<x 

LA 

X 

o 

« 

CM 

o » 

ro 

o 

cn  a 0 

M 

W O 

►- 

o 

X 

M 

K 

OOOOOOOOOOCJOOO 


(O  » 

z » 

o » 

l-»  * 


O tf> 

n I 
Z 


X ►-» 
> X 
o & 


«•  < «« 

H 

O X 
O M o. 


O X CM 

o X a.  1-4 
^ Z «x  s. 


» 

o o o 


UN 


145 


OErERNINE  VERTICAL  ANO  RADIAL  DEFLECTIONS 
DELTA  (I)=RAD»A0CH  U «-TER<  JNOO) -VI 
OOELCI)=«SFEEO*DVX( JN00)-DVHI/A9CT(I) 
IFdSCHEM.Ea.OJGO  TO  31 


ro  tn 
o o 


o o 


u Ul 
O ^ -J 

o * • 

^ CJ 
• MM 

o a.  CL 


• <VJ  fy 
. V N. 


X «. 

o 

n ^ 


o r 
c o 


a a 
< < 
I I 


X X 

0.  a 


o o 

L5  U5 


o o 

• • 

U M 


M O Q 

< 

in  m 

♦ ♦ 

0.  o o 

M 

1 

o O o 

O 

UJ  uJ 

♦ » 

1 < < 

•H 

o o 

•4- 

o 

o 

♦ ♦ 

•!^  ♦ 1 

lU 

M 

II 

V 's. 

• 

• 

♦ ♦ 

<«#•<» 

o 

w 

o 

M M 

O 

o 

o o 

♦ * 

M M M 

M 

w 

M 

'M 

90 

a o 

M 

M 

M 

z 

z 

M 

♦ ♦ 

^ ^ 

w 

♦ X 

w 

tn  in 

M 

< 

< 

a 

4 

» ♦ 

XXX 

< 

a. 

0. 

o 

o 

a a 

C 

O 

M 

o 

O 

IT 

• 

• 

M U. 

K» 

♦ ♦ 

h-  M M 

M 

M 

< 

r 

M 

M 

< <« 

o 

o 

w 

M 

a 

i£> 

M 

M 

M 4 M 

♦ ♦ 

a Q o -1 

w 

w 

<c 

» » 

M 

< 

vO 

K> 

«» 

o 

D 

M 

a 4 0 

O 

* ♦ 

«t  < < 

Ui 

z 

o o c 

X » « 

A 

U 

o 

M 

O 

M 

a 

o 

a. 

u •« 

M 

♦ • 

M 

1 1 1 

o 

{/» 

< 

MOO 

M ^ 

M 

M 

a 

o 

o 

M 

O 

w 

• 

• 

U. 

M o CO 

♦ ♦ 

K- 

A 

(«)  W M 

4 

o 

< 

X 

o 

o 

M 

n MO 

n 

♦ * 

1- 

M M M 

fO  < 

• 

a 

a 

a 

Z O' 

<r 

Qt 

w 

» 

a 

M 

M 

M 

O 

M 

• 

• <M  O 

1 u 

o 

♦ » 

a. 

^ 

M 

"•v 

A 

• 

• 

ijj  a^  ly 

U 

<VJ 

• 

w 

Z 

o 

o 

a 

M 

M 

w o V 

«» 

» ♦ 

o 

M M M 

w 

O «« 

M 

a 

o 

O' 

~i  ^ «i 

UJ 

UJ 

«J 

M 

a 

(y 

«» 

O 

o 

«« 

o 

tf 

_J 

•>  z 

z • o 

o 

♦ CO  • 

XXX 

tt"  M M 

w 

0 

• 

• 

u uj  a 

o 

w 

• 

-j 

M 

o 

o 

N 

0^ 

1 

• 

• 

M < 

4 <NJ  CO 

• 

♦ z ♦ 

(Lao. 

»- 

M 

UJ 

M 

MOO 

• 

• 

♦ 

^ < 

UJ 

uJ 

w 

«» 

• 

M 

«• 

«» 

4-  K 

M 'V  CO  X 

o 

♦ o ♦ 

< c <s 

UJ 

O X 

(/5 

O 

o 

O 1 t 

M 

M 

Mi  M 

u 

a 

X 

fVI 

a 

a 

• 

M 

M 

M 

* V 

V a X N 

• 

♦ M ♦ 

jt 

WWW 

O 

c5  n 

a 

• 

• 

« a a 

in 

M 

O -1 

• 

1 

M 

e 

z 

UI  Cf 

M 

w 

a. 

a. 

M •» 

4 M dk 

UJ 

♦ M ♦ 

K 

<A  in  tn 

«>  o 

< 

«* 

A 

M « 4 

• 

• 

O UJ 

o o 

3 

UJ 

M 

u. 

u 

4 M 

M Z • to 

-J 

♦ 4 ♦ 

>fc 

CD  CD  ca 

M < 

1 

M 

M 

M 

z QC 

H 

Z O 

M 

< 

«X 

z 

O o 

T 

• 

X M M 

0m,  W 

•••♦MX 

• 

» -J  ♦ 

< < < 

• 

• •• 

«« 

w 

w 

w 

o 

M 

M 

J u> 

II 

w 

Of 

II 

• 

oc: 

•• 

• UJ 

X 

M 

Bu 

a o M z 

»o  Z M U X 

u\  UJ  ♦ 3 ♦ 

K 

WWW 

C5 

a ** 

M 

0^ 

ac 

U II  •• 

w 

w 

UJ  •» 

« 

w 

«« 

M 

90 

o 

o x 

o 

M 

X r* 

Z <C 

w 

w 

0.  -J  4 -J  a • r- 

M m 3 ♦ U ♦ 

(O  (/)  t/)  UJ 

■At  *-* 

W 

M 

fVI 

M M 

-J 

II 

M 

O 

M 

H O 

« 

M 

o 

o 

II 

(/)  CO  u.  UJ 

UJ  a CO  II 

Z ♦ ~J  ♦ 

• • 

• 

UJ 

UJ 

UJ 

UJ  M M 

U 

-1 

o 

'd 

w 

M 

• 

r> 

Of 

w 

a 

X < o 

O M 

0(BM30a4X«» 

O M ♦ 4 • 

o o 

•s  « < 

t-4 

M or 

O 

o 

a 

a 

Z W. 

< 

< 

M 

M 

r m 

a 

M 

M 

M 

M 

O 0^ 

M M 

M 

O CD  M 

0- 

< < 

a M 

M II  ••  M M 

4 M 1-  ♦ O ♦ 

H H 

W w 

w 

u 

a 

w 

w 

w 

M M rj 

w 

w 

w 

M J 

w 

w 

X 

M 

w 

M 

w w 

w 

< 

z 

lu 

M 

w 

M (0 

(/)  CO  «•  w 

M 

Z • • 

X Jfe- 

u.  u.  u. 

li. 

u.  u 

< 

u. 

u. 

U. 

r ^ u 

U. 

J o 

X u- 

u. 

^ fr 

II 

u. 

O 

P 

UJ  -J  u.  u. 

a 

II 

o 

« 

M U. 

u x a X X u.  ^ 

U O O ♦ UJ  * 

X ^ 

M M M 

M 

M c 

w 

M 

M 

M 

< <x 

M 

M 

0X 

o 

o a 

M 

< 

< 

M 

M 

o 

UJ 

a 

« M 

M 

M M O 

M 

0 

M 

M 

M N a M M < 

M o u • u * 

u.  u. 

1- 

UJ 

(X 

in 

• 

M 

• Z ♦ 

a. 

90 

M 

> 

a.  M 

Ui 

ar  e 

fii 

M M 

ta  » n * 

<x 

to 

UJ 

a w 

to 

O 90 

< 

4- 

4 

4 4 

w ♦ u.  ♦ 

o 

a 

Z 

0. 

♦ # 

o o 

u o 

o 

o o 

u o 

o o u 

z> 

c 

o; 

m 

3 

iA 


146 


SU9P0UTINF  ADAPT  7A/7A  0PT!1  FTN  Cl/ia/78  11.03.5A  PAGE 


O 

4 


LJ 

O 


r Q 

I M Q.  < 
: -i  ^ t 


• 0£ 
O -J 
II  UJ 

> <«  o 
• M 4 

> — O 

I 0^ 


O 

o < 

• M 


I ^ X 

0^  <t  a 

-J  o < 

I UJ  ^ PX 

I O U.  (/I 

’ ^ » a 

4—0 
— M 4 

r — — 


— X X 

M & a. 

«x  «t  < 

X 'X  »x 
►-  Z 
O M O 
01  C ^ 

f 4 4 


M O Cl 
XII 

a u.  u. 

<t  ♦ 4. 


I LiJ  » 


M I-  « 
- ^ ( 
in  M liJ  PX  »-  U.  * 


in  »X 

uj  a 
. Of  u. 
4 a.  X. 


o as  ^ M I 


. I-  a w o 
I I o II 
'04CD  — 

• <I  U-  »-* 

f Of  a * 

. X.  Q ^ 

< M Z 

: in  4 px  u. 

lODin^  — 
I 4 OD  ^ o 

• H>  » U.  • 

> a a <x  a 

1 a < ♦ * 

«I+px  •j<ujl3<Q'-«.UJ 

• 4 4 M ^ 

• <1  O • 

i^coa  — 
’ < < li,  »-• 

r 4 H II  px 

i u.  — — o 

J M M M z 


O < M C/  UJ  £> 


• 4 
— O 

M ^ • 

PX  ^ I 

Or  t«  < 


I in  — 
■ CD  tH 

> 11  *x 

• — a 
r »-■  u. 


lOQ.U.U.Vl’U.tLtxfflZU. 


-XX 

0.  a 

i~i  < < 

<n 

CD  o M 
u.  o in 

« 4 4 


t-  i9  11 

z z 
u.  u.  u. 

w ♦ I UJ 

••  V X 3 

^ >-  X z z 

M U.  U,  t-t  ftf 

•X  II  II  t-  => 

O V X Z Q 

X X X O UJ  Z 

U.  U.  U.  O Of  UJ 


147 


z ♦ 

V » li. 

♦ W" 

CM  X 

ro  X 

ro  X 

to  ^ ^ 

a CM  •> 

(0  1H  O 

X » 

X » 

X * 

tlJ  • V 

» 0 

M 

Z 

2 ^ » 

» z 

» z 

* X 

O CM  •* 

cr  to  I 

to 

bJ  » <x 

z to 

Z CO 

Z to 

AT  Z l-l 

> >•  u 

& 03  h» 

to  » 

to  * 

to  * 

O *0  I 

* • » 

& 

ro 

to  to  IJ 

» 

* ^ 

» A 

Ul  ♦ o. 

^ wo 

^ T 

cc 

D ♦ r 

CM 

^ wo 

M ro 

ivt  O 

s.  »l  X 

ty  ♦ 

(V 

to  — w- 

CM  CD 

m tn 

ro  OD 

Z X ♦ 

S UJ  u. 

1-  CO 

* 

< >0 

z 

z 

z 

Z X 

Z K 

Z 0£ 

c ♦ •- 

♦ a ♦ 

X T 

ro 

U 1-  » 

t9 

o 

o 

-0  ♦ 

•o  ♦ 

"J  ♦ 

IH  CM  X 

^ • CM 

m ♦ 

z 

O U -i* 

M 

M 

M 

1 CM 

1 ro 

1 

to  X Ai 

>•  WO  X 

M M 

X fO 

to 

to 

to 

CM  U 

ro  Ul 

z ♦ • 

• u u. 

or  X 

• 

to  1-  to 

» 

» 

» 

^ Ul 

-J  w 

^ bJ 

-I  CD  Ul 

a.  »H  X 
(/?  u » H 
D -I  m N. 
(/)  «x  X UJ 


fc-  I Si 

S Ui  X z 
o z o 

»■  M _J  X 
OD  W < X 

:3  D ut  c 
£/l  O X O 
or 

CD 

o 

to 


^ Ul  • ^ » ro 

>>0^X0)  ^ 

o ' • > * i-f  • 

» eg  U.  cvj  • <\j 

fo  j • O'  a «iQ 

> uj  ^ ^ X a.  O' 
o > X • <r  • 

ru  * u.  * X i*) 

> ^ » 4"  M ►-  2 

n X CM  W • -5 

V > > X • • 

CM  CM  u,  ii.  X Q.  OJ 

> O > • O 

V • UJ  ^ K 

^ W O wo  t3  to  2 
o X Ck  X I X to 
4i  > o u.  iii  ••  a i 
^ (NJ  U.  » > O to  < 
•s  O X to  S.  o K 

Z • X flo  z • z < 

O CM  O U.  O O O 
^XXX  ♦XU'  XI 
XXXXCOXujXt 

>>o<MOu.caoi 

cuauu.tJtoo( 


I n ^ 

I to  to 

X > 

I Q O 

I I fM  CO  > 

t fO  ^ ^ ^ u 

rxxuiuiui- 
r >>  V O O O 

> o o » » * 

I It  II  CM  ro  n 1 
J CO  ^ X V J 
I u II  *•  II  « 
J UJ  UJ  CM  CO  I 

> O O X bC  » 

3 o a u.  u.  u.  ( 


Wo  UJ  O 
. U • td  • 

I < e • CM  I 

: w »•  CM  ^ 

. z CM  ^ UJ  f 

O & Ul  O I 

. u o B I I 


)u.criou.sa(ou.u.cou.u.<oh.b. 

MU.  MlkOU.»-tMU.MMU>MM 


148 


C***»ISCOUS  BAHPER  FORCE 
8PV*B2VJ 


SUBROUriNE  SUSP  74/71*  ORU 1 FTN  4.6*4333  01/10/78  11.00.64  PAGE 


CJ 


o ift 

%£)  iO 


149 


90v«  iS'cc'TT  H//?T>TQ  prf*>+q'‘i  Nid  T=idO  31308  3Minoa8ns 


V X 
* >- 

A » 

C 

h-  X 

1/1  V 

V a 

W CJ 


(/)  -J!  u 
UJ  » >• 
O »0  Q 

or  r • 
0*^0 
u.  • > 

fVJ  • 
ti  Z *-t 
O ^ X 
»*  » A 
</i  fA  n 
^ * 
o * »» 

Z C\4  X 
C ^ 0 


O » IL) 
O W X 
. CD  X 

* O 
U1  (Sj  * 
h-  I « 

J U CD  Ul 


) J CO  "S 
1 «;  X ui 


X 1/1 

►-  > r 

M • 5 

o ( 


T >•  ' 

a * 1 

CJ  M t 


M ^ I 

o a 

» CO  f 
X 1 


* tA 

\ u. 


O AJ 

*■  & . 

o (A  : 


»>  * • ^ 

- -I  T «V  • 
3 uj  tb  ro 
i-  G X X 
r • fvj  • • 

- W T fA  AJ 
^ U.  (X  Z 

> ui  • ♦ 

- o ^ X m 

3 • > • 

» (NJ  U.  (M  * 
5 j b Of  03 

- uj  ro  »-  r 

b G > X 


M 
fO 
a (X 
*s  • 

■ O CM 

c b-  Z «« 
► (/I  T < 


- ^ IM  t 

3 X (M  n ♦ ( 

» > > V X 

J(VjU.U.Za.^ib  UJ 

■ov  •voi^x  — I — 


CD 

CD 

U' 

T 

X 

T 

ll< 

t-l 

Q. 

* 

* 

I 

X 

O 

K> 

lO 

b- 

fi 

» 

X 

X 

CD 

CD 

w 

* 

» 

» 

X 

a 

0 

o 

CL 

o 

o 

o 

b- 

» 

O 

►“ 

1- 

UJ 

bD 

» 

CD 

fA 

« 

X 

X 

G 

X 

X 

1- 

w 

CD 

:d 

* 

Ul 

CD 

li. 

♦ 

A 

X 

» 

1 

u. 

G 

CD 

w 

o 

O 

C' 

O 

X 

•jj 

♦ 

4» 

» 

1 

0. 

X 

CD 

to 

A 

(D 

w 

A 

I 

X 

CD 

»- 

* 

O 

X 

•«» 

X 

w 

G 

UJ 

» 

X 

» 

» 

* 

» 

G 

X 

» 

< 

t 

Cl 

to 

m 

K* 

4> 

w 

K> 

w 

G 

0X0 

X 

(D 

o 

>- 

IJ 

IH 

» 

* 

» 

» 

_J 

» 

CD 

X 

N. 

0. 

0. 

0. 

0. 

(M 

u. 

w 

u. 

o o o 

o 

1 

o 

X 

bw 

1- 

b- 

» 

M 

» 

« 

CD  CD  tA 

CD 

CD 

» 

X 

< 

X 

X X 

X 

I 

O 

(M 

a. 

w 

X 

G 

• 

1 

« 

t 

u 

W 

UJ 

» 

T 

M 

•» 

M 

II 

* 

U. 

» 

I 

• 

V 

CD  CD  VD 

CD 

1 

G 

w 

(M 

u.  u. 

G 

G 

«* 

4 

CD  < 

o 

1 

rsj 

u. 

u. 

U. 

U. 

b* 

to 

» W 

» 

T 

♦ 

l«l 

X 

^ l>> 

T 

CD  CD  CD 

CD 

z 

u. 

W I 

ro 

S. 

Z 

(O 

X 

X 

Z 

Z 

u- 

» 

X X. 

CD  « 

O N 

• 

• 

• 

• 

« CD 

«» 

U.  CD 

♦ 

w 

X 

• 

ul 

UJ 

UJ  tm 

UJ 

CD  H>  » 

< 

• * 

< UJ  » 

(M 

e 

o 

-J  o 

o 

G UJ  O 

w 

«» 

o 

w 

z o 

♦ 

1 

1 

« 

• 

• 

• 

u.  z o 

U) 

«7 

X o 

UJ 

w 

X 

G 

0. 

a 

<*» 

•• 

I 

»- 

0.  4> 

X o » 

w 

«t 

« 

u. 

u 

G 

G 

CD  CD  4 

w 

UJ 

• 

CD  « 

b> 

• 

N 

o 

o 

o 

O 

O 

lb  • t“ 

CD 

X 

(M 

» ^ 

o 

< 

M 

Ul 

K 

w 

u.  ^ UJ 

» 

w 

V 

< UJ 

» 

w 

1 

X 

(/I 

tA 

CD 

CD  CD 

CD 

b o z 

ro 

o 

(M 

* X 

UJ  z 

» 

* 

G 

G 

G 

m 

(MOW 

* 

• 

Jt  W 

b»  X z 

z 

tm. 

a 

«X 

<r 

«x 

W ^ O 

• 

o 

fVJ 

X (A 

l>' 

b- 

X 

M 

< 

U.  oi*  » 

O bO 

CD  » 

z o < 

» 

I 

o 

e 

o 

o 

to  ^ 

w 

♦ 

* 

u.  ^ 

w 

* 

X 

X <M  « 


D I 

o : 

Of  1 
CD  I 
G : 
(A  < 


b- 

> r* 

(D 

>o 

() 

CD  Z 

X 

T 

UJ 

w 

• • 

• 

U. 

(D  Z 

0; 

X 

X 

I CD  CD 

0.  z 

a 

X 

X 

o o 

o 

o 

b 

t=- 

UJ  > 

C 

u. 

UJ 

• 0,  o 

w w 

w 

w 

a 

• • 

• 

• 

W (M 

CD 

CD  (VI  lb 

• 

> 

O CD  o CD  CD 

z 

UJ  UJ  UJ 

UJ 

CD  X 

to 

V 

>*  o 

N. 

CD  N 

O X M 

O II 

W 

CD 

» 

» 

U;  CD 

z 

z • 

Z 

G 

z 

• z < 

O «* 

CD 

bJ 

CD 

CD 

• 

• 

• * 

• 

• 

O lb 

z 

o 

O (M  O 

u.  o 

a o b- 

II  W 

•• 

O 

G 

U, 

o 

o 

u.  u. 

G 

G 

Of  tl 

Z 

* z 

z 

• 

z 

UJ  X iL 

W UJ 

•-« 

Of 

CD 

CD 

II 

II 

O J 

O 

O 

o ^ 

c 

z 

z >• 

Z CD 

z 

UJ  X z 

X X 

X 

O 

II 

•• 

CD  CD 

w 

w 

U,  tti 

II 

o 

V 

O (M 

O 

u. 

c a o w 

O 1- 

Q 

u. 

« 

tl. 

U. 

G 

u,  u 

u. 

u> 

X 

0 

CD 

O o o 

O U. 

o 

CD  o o 

ti  CD 

CD 

O 

O U.  lb 

W M 

M 

W 

Ul  U.  X 

w 

& 

M 

(/)  CD  ij 

U.  X 

C/>  ^ 
_ I C3 
h-  * 

Uj  — i «• 
X ««  G 
^ O 

u » •C* 
* CA  ^ 
fO  U.  W 
^ U.  (/) 
••  II 
CD  o ^ 
X X I 


< o 


X ^ 


o 

o 

X 


!-•  CA 

x ♦ 

Q.  < 
o ^ 

CD  CO  »H  • UJ 
» U.  X < Z 
(Nj4-a»»-> 
-i  l*>  (A  ^ u 

» (A  o CD  ro 
(M  U.  c/)  U.  ^ 
X w II  II  •• 
U.  o UJ  u.  o 
I X X X X 
I (M 


M I 
Z kr 
CL  CD 
CA  » 
» to 

m ^ 

CD  I 
II  II 
X M 
X X 


«>X  • (VI 

h-  G ^ 

UJ  Z X 
X V M 
►-  O V 
C3  X G 
» » X 
L9  • » 

CD  (VI  O 

♦ X z 

b->  *•  • Of 

X • 3 
OL  MO  (M  ^ a 
O UJ  z 

• U.  V Of  UJ 


150 


SUBROUTINE  PSTORE 


3 

O 

a 


O' 

o < 

U.  ! 


k*  4-  CO 
hi  t-  U. 
X tr  •• 
h-  >•  M 
OJ  X 

O ro  t/) 
► K li. 
M (/)  » 

X >•  <\J 
0,  * X 
(M  M U. 

n ^ 

•*  CO  4“ 
I-S  > V 
V-  » e/i 

cvj  > u. 
a a » 

• 

CD  V V 
► • (/J 
(M  ro  u. 
o a •> 
CO  <vj 
CO  V 

>-••14. 
O ^4  * 
• 2L  4- 
ff*  CO  X 
>■  V u. 


IX 
« «r 
Ui  Ui 
O'  q: 

•0  X 


o 
^ o 

w (.r\ 
(X  • 
i.i  o 

Q H 


CO 

X 

0. 


: X 

r 0. 


♦ fo  ^ • — 
> -j  -r  »v  I, 

O Ui  u.  K 3 

#0 

^ CV>  • ► 


a > 
UJ  c 


‘ .J  Uk  CX 

:s 


0>  o 

lU  T^ 
X 


» (M  U.  (\l  K -J 
W -»  •■  CX  C 0- 

> Lu  ro  ^ •>  z 


I ^ O'  If* 
• »-l  tf\  ■*•< 
I • -H  ^ 


m 

o 

> z 

UJ 

to 

a 

_J 

V 

• 

• 

PO 

CM  CM  CM  CM 

-to 

toto 

•i 

UJ 

to 

X 

X 

o 

to  A 

C3 

m 

V 

N.  V.  S.  N. 

to 

n 

fT' 

CO 

> 

4- 

»H 

0. 

If' 

(/I  O 

o 

« 

V 

CM 

PO  4 CM  kO  4- 

Z 

CO 

z 

to— 

Ul 

O 

z 

CM  fO 

u> 

X 

w 

X 4- 

« 

• 

o 

a ck:  a (X  a X 

a 

X 

I— 

a 

k- 

UJ 

I 

* 

> 

> V 

X 

X 

M 

Cki: 

a.  w 

» 

00 

»-  1-  k»  W 

a 

z 

k- 

4 

.J 

iV 

V- 

CJ 

CNJ 

u.  u. 

o 

z 

4 u: 

h- 

• 

<x 

« 

to4 

z 

z z z z z 

k- 

a 

a 

o 

V 

> 

a 

N.  to 

♦-H 

N 

k- 

N -J 

<x 

kto 

rH 

* 

PO 

V V V V V >to 

pE' 

4 

o 

CM 

4 

a 

-1 

a 

a 

UJ 

V 

to 

UJ  4 

z 

k- 

0. 

V. 

u)  a 

k- 

U' 

» 

CM 

z 

4” 

z 

CM 

PO  4 CM  to  4 

CM 

>- 

>■ 

_J 

toJ 

a a 

a 

V 

X 

fO 

»- 

s» 

»o 

CO  fn 

_l 

kto 

a 

(/i  ^ 

IjJ 

X 

T 

o 

CM 

z 

PO 

>> 

Z 

V 

> > X X X 

V 

CO 

CO 

CM 

pn 

4 

tr 

CM 

U; 

\l^ 

Ul 

fc- 

a 

z 

It 

II 

• 

<r 

L5 

z 

X X 

z 

0. 

to 

z 

< 

I 

k> 

I 

U 

a 

V 

z 

>• 

Z 

a 

>- 

u. 

a a a u.  u. 

a 

a 

a 

>- 

>- 

>- 

>- 

O 

n 

n 

o 

k 

-to 

-to 

•to 

a 

c 

U: 

>- 

O (4. 

>- 

>- 

^ a: 

X z 

o 

a 

U 

>- 

o 

>- 

O > 

•• 

II 

H 

II 

- N II  H •• 

ti 

II 

II 

It 

II 

11 

n 

II 

II 

II 

"0 

1- 

k- 

a 

u> 

(O 

CO 

M 

u » 

o 

X a 

o 

a.  o 

l> 

•• 

M 

II 

II 

•• 

II 

II 

II 

-to 

<ito 

-to 

-to 

«> 

-to 

* 

-to 

II 

to- 

X 

a 

X 

• 

V 

N*  Q 

N.  CO 

K 

to 

V 

V M 

•«to 

.to 

toto 

-to 

O 

^ CM  PO 

4 u\  a>  a 0 

o 

CM  PO 

4 

iO 

K 

CD 

a'  o 

n 

IH 

X 

X 

u; 

z 

Z 

to 

Z CD 

z 

z 

O' 

z 

Z CO 

CM 

ro 

4 

vD 

N 

«e 

O' 

toC 

•to* 

to4 

toH  toH  W toH  to^ 

w-t 

CM  CM 

CM 

CM 

CM 

CM 

CM 

CM 

PO 

PO 

.J 

to 

z 

to 

«tf 

z 

o 

O «\J 

o u. 

o 

O o 

CO 

o z 

to- 

w 

to- 

to- 

to- 

to- 

to- 

to- 

w 

to- 

to- 

to- 

to- 

to- 

to- 

to- 

to^ 

to- 

a 

-> 

k-t 

O' 

X 

M 

z 

4- 

z 

Z 

X » z 

X 

to 

r 

z 

Z UJ  »- 

U> 

k- 

K 

►“ 

hto 

k- 

k- 

1- 

kto 

»-  I-  k-  t* 

»- 

kto 

k-  fer 

k-  k- 

kto 

h- 

»- 

k* 

k- 

z 

to- 

kto 

w 

M 

k* 

z 

z 

z 

>• 

X CO  z 

z 

r 

O' 

z z 

.J 

-J 

4 

-J 

-J  u 

.4 

toJ 

4 

4 4 U -J  -4  ^ 

-1 

_l  .J  -J 

4 

toJ 

.J 

II 

u- 

Z 

k- 

to- 

3 

o 

>- 

o 

<\J 

o u 

o 

o 

0- 

o 

>-  o M a 

0. 

a 

a 

a 

a 

a 

a. 

a 

a 

a 

a 

a 

a a a a a 

a 

a a 

a 

a a 

a a a a 

a 

a 

a 

o 

a 

-t 

O 

-J 

a 

C 

CJ 

O CJ  O 

<J  b.  O 

o 

► 

CJ 

» o o >- 

>• 

>■ 

>•  >- 

>■ 

> 

> 

>- 

> 

>> 

^ IP-  >• 

>- 

>- 

>- 

>-  >- 

>• 

> V o 

r a 

u 

a 

k4 

o 

V) 


151 


SUSRflUTINE  I/PLOTS  ^4/74  DPTIl  FTN  4.6»433d  C1710/7a  11.0  0.54  t»AGF 


a 

— 

04  — 


a 

z 


^ z 
) uj  a *-i 

I Z Z -4 
♦-4  V » 
• ^ k»  Q. 
“3  -4  ^ 

o a 
: a >•  j 
I CD  X S. 

' 3 N.  » 

I (<0  Z O' 

o 


3 CL 


<X  ^ CL  M 


' *4  <\i  • 


O < J II  H I 


I «■'  LD  1-1 


*-*00.  Q.M-P30Z*-*. 
(/I  -J  3 3 


I < • II  *-  UJ  < U4 


or  O Z CO  Ul  II  ^ ( 

oa  z z .j  ir  z u . 
3 •-*  o a M a.  0.  I 
to  (r-  u *0^3* 


. M II  a M M 


a.  eg  ••  I-I 


I 3 X 
I z <r 
X 


I o . 
> (O  o < 


» O «>4 

• f ^ 

' CL 

► Z 3 

r J ^ 

a.  <f  -* 

► to  «.  «4 

4 • I 

T O 1 •• 

4 O'  • 

• X ► J 

• T ♦ (T* 

r « 

► ♦ # z 

• • «o  a. 

9 -3  • V 

► ►•  Ul  ♦ 

• • .J  ^ 
s O Of  ^ 

r •■»  <r  *0 

■)  CO  I-I  «•  » 

■*  *^  « in  o 

e X < • » 

t ^ ^ ^ ^ 

"r  m ^ to 

4 U O in  -f  I 

J -J  » rj  a. 

c 4 • • 3 • 

) O a o » o 

. ► I < ► 

■I  o *«*'•>  • 

■ • o O'  m <rt 

J cr  UJ  «- 

J ui  to  cc  UJ  »- 

• • «-•  z z o 

X X 3 l-l 

< < z ^ a z 

X or 

•-•  -»  -J  -I  J 3 

-i  ^ J ^ o 

U < <X  4 < UJ  z 

•-*  O U O O ftT  UJ 


t4  cn 


e tn 

<VJ  CM 


152 


o 

<r 


:>  II  ^ o o 

•<  M O 7 UJ  eo  Oj 

W Qi  Z 13 

roZ'-'Z»-^z 

} O w UJ  l-t  4 < 0^ 

: <i  z O'  i-  o 

)OuJt^OOOUJ2 
lOQ^ZOb.li.ft'UJ 


153 


Lj 

O 

<f 

Q 


iTv 


n 


-I  li‘ 
U *0 

< ^ 

• < 

►*  « 
< Q 
O •> 


r m * . 
r < ^ : 


•H  a:  X •H  3 


I UJ  (M  •• 
1 m <-i 
I «t  • X 
: -J  ^ 'O 

I » 

’ O «H 
' «3  ^ ^ 

I O o 


O 

UJ 

or 


o 

. 2 

: o • 

4 »-«  O 

: >.  oc 

- X T 
O «« 
J X o 
: X < 

: o tu 
u ot 


• «I  V 
H>  O <C 

X fo  r«9 


UJ  « < 

K r X D 
•-»  ar  a *-  Q 
ar  o o u z 
z u.  u.  a lij 

C3  O 


154 


SUBROUTINE  IOLAB2(JS) 

READ  AND  PRINT  WITH  A LINE  SKIP 

COHH0NiriONUH/NWT,RRO,LABEL(8  0) 


o «o 


<c  » 


■J  llJ 
UJ  CD 


<t  -J  V 

— C < 

(/>  -H  O O 

"7^3  — GO 

O <T  •• 

^ O'  » « X 


O z UJ  CO  CU 
2 

O Z w 2 h-  t-  Z 
C3  «-•  UJ  M < < tt 
C *-  h-  Z Z 3 

< M z or  or  a 

OUJO'OOOUJZ 

QCk^ZOU.U.(ZUJ 


155 


APPENDIX  B 


DETAILS  OF  THE  CONTROL  SYSTEM 


Details  of  the  control  system  software  included  in  this  appendix 
supplement  the  discussion  in  Chapter  3. 


B.  1 Principal  Nomenclature 


A/D 


Ascn 


Console  Device 


Core  Image 


CPU 


D/A 


DAM 


Data  Switches 
(Sense  Switches) 


DMAC 


EAI 


Analog  to  digital  converter.  A hardware  input  device 
which  converts  analog  voltages  to  digital  (binary) 
number  for  computer  utilization. 

American  standard  code  for  information  interchange. 

A binary  coding  standard  for  representing  characters 
and  symbols  in  digital  form. 

An  I/O  device  which  is  the  normal  operating  system 
communication  channel  between  the  user  and  the 
computer. 

The  binary  version  of  a program  in  a form  that  the 
computer  can  directly  execute. 

Central  Processing  Unit.  The  hardware  in  a computer 
that  controls  system  operation  and  performs  logical  and 
arithmetic  operations. 

Digital  to  Analog  converter.  A hardware  output  device 
which  converts  digital  (binary)  numbers  into  analog 
voltages. 

Digital  to  Analog  Multiplier.  A special  type  of  D/A 
which  multiplies  a D/A  output  by  an  analog  voltage 
reference. 

A group  of  eight  switches  on  the  computer  front  panel 
which  can  be  set  on  or  off  and  tested  by  the  computer 
for  logical  processes.  (See  Figure  11.) 

Direct  Memory  Access  Channel.  A hardware  device 
which  allows  high  speed  (one  million  words /second) 
transfer  of  data  to  and  from  the  computer  memory 
without  effecting  the  CPU. 

Electronic  Associates  Incorporated.  Company  which 
supplied  the  hardware  for  the  control  system. 
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EOF 

EOR 

FILE 

I/O 

K 

Logical  Unit 

Monitor 

Object  Code 

Octal 

Operating  System 
Record 

Source 


End  of  File.  A code  which  the  computer  recognizes  as 
the  end  of  a set  of  data. 

End  of  Record.  A code  which  the  computer  recognizes 
as  the  end  of  a discrete  record  in  a set  of  data. 

- A set  of  data,  blocked  into  records  and  terminated  by 
an  EOF,  stored  on  a disk.  The  file  is  identified  by  a 
file  name  and  may  contain  programs  or  data. 

- Input -output. 

- Symbol  for  binary  one  thousand,  1024. 

A logical  unit  number  (octal)  associated  with  a data 
set  or  an  I/O  device. 

Supervisor  program  which  controls  system  resources 
and  utility  program  operation.  Decodes  user  commands 
to  execute  operating  system  function. 

Binary  code  produced  by  the  FORTRAN  compiler  or 
ASSEMBLER.  The  object  code  is  input  by  the  loader 
and  linked  together  to  create  a core  image  file. 

A base  eight  number  system, 

(See  Monitor) 

A physical  or  logical  partition  of  data  for  input  or  out- 
put. Typical  block  of  data  read  or  written  during  one 
I/O  operation. 

Program  code  written  by  the  user  for  input  to  the  com- 
piler or  assembler. 

Carriage  return  character. 
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B. 2 Computer  Programs 


The  control  system  software  incorporates  two  programs:  (a)  file 
builder  program,  and  (b)  control  program.  These  programs  are  described 
in  the  following. 

B.2.1  File  Builder  Program 

This  program  reads  the  paper  tapes  containing  the  axle  dis- 
placement records,  packs  the  data  in  a buffer  and  writes  them  on  a disc  file, 
as  shown  in  Figure  9(a).  The  flow  chart  of  the  program  is  given  in  Figure 
B-1. 


B.  2 . 1 . 1 Listing 

The  file  builder  program  is  coded  in  EAI  FORTRAN,  and  a 
listing  of  the  program  appears  on  the  following  pages. 
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Figure  B-1.  File  Builder  Program  Flow  Chart 
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Figure  B-1.  File  Builder  Program  Flow  Chart  (Cont) 


Index,  Value 


Backspace 
File  1 Block 


Figure  B-1.  File  Builder  Program  Flow  Chart  (Cont) 
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C file  HUILOER  PROljRAM 

C INPUT  paper  tape  UATA  to  DISK/PRInt,AND  file  patcher 
C data  on  paper  tape  is  ten  time  steps  per  RECORo  on  disk 
c -1.0  IN  time  slot  Marks  end  op  data  set 

C PROGRAM  requires  LOGICAL  UNIT  ?1  OCTAL  BE  POSITIONED  TO  DATA  SET 
c data  switches  A*t)»C»D  USED  FOR  PROGRAM  CONTROL 
C A ON  prints  input  TAPE  AS  REAR 
C B ON  ENDS  PROGRAM, abnormal  E^IT 

C C ON  pause  PROORAM  before  tape  REcoRO  input, each  RECORD  » 1 STEP 
C D nN  SKIPS  TAPE  INdUT  AND  LINkS  Tr>  SCANNER  SECTION 
C DISK  I/O  DONE  VIA  QMONO  CALLS, 44  El  EPEnTS  PER  RECORD 
C LAST  LINE  oF  RECORD  fs  ZErO 
DIMENSION  aUFF(ljO) 

LOGICAL  SENSW 

NTNPsZ 

NTYP=1 

N0ISK=l7 

C octal  2l 

NPrNT=16 

NTApE=4 

WRlTE(NrYP,B300) 

9300  F0RMAT(‘iX,22H  FILE  CREAtIoN  PROGRAM,/, 

1 22H  set  data  switches  A-0»/f 

2 32h  hit  carriage  RETURN  TO  EXECUTE. »/> 

READ  (NINP, 93-50) 

9350  format (IX) 

REWIND  NDISK  » 

IF (SENSM (4) ) GO  To  2100  * 

00  lOOO  le'f *10000 

C clear  buffer  To  2tRO 
00  SO  jal,50 
0UFF(J)3O.O 
50  COmTInOE 

DO  100  J=l,10 
Kla ( J-1 ) *4*1 
K2sKi*3 
75  iao=i 

IF(SENSW(3))  pause 

rEaD  (NTAPE,9  0 0)  <BuFF  (IjK)  ',IJk»K1,K2) 

9000  FORMAT (F7,4,3F0,9) 

IF(SENSW(1)  )WRITE(NPRNT,95P0)  (nUFF(lJK)  ,IJK3K1»K2) 

9500  FURMAT{5X,AF15.5) 

IF(Bi)FF(K1)  .NE.-l.o)  GO  To  106 
WRITE(NrYP,90l0) 

9010  format (5X,31H  -1  INPUTjTTrE  0 TO  SKIP*ELSE  1,/) 
REAO(NIWP,9020) IGo 
9020  FORMAT(Il) 

IFdGO.EQ'.ol  GO  To  75 
100  continue 
C DUMP  DATA  buffer  TO  DISK 

CALL  QMOND  (16,NDISK,8UFF  (1 ) ,Bi|FF{44)  ) 

IF  (RIJFF  (K1  ) ,E0.-1 .0)  GO  To  2056 
IF(SFNSW(2))  go  io  3000 
1000  CONTINUE 
C end  file  ,CL0Sf  UNIT 
2000  ENDFILE  NDISK 
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c rewind  and  test  EILE 
DO  2025  iBltAA 
2025  Bi)FF(I)=0.0 

2050  re»jind  NDISK 

2100  CALL  QM0N0(]7»NDlSK.RUFF<1).B||FFtAA)) 

C TEST  data  to  LOOK  FOR  BAD  POINTS 
DO  2200  1=2,40,4 

IF((BUfF(I).(3T,1.o).OR,(8uFF<T).LT.-1.0))  (jO  TO  2300 
IF( (0UfF(I+1) .Oi.l.0> .OR. (MUFf(I*1) .LT.-1,o) ) GO  To  2300 
IF  ( {R:JFF  ( 1*2)  .GT.l  .0)  .OR.  (BUFf(I*2)  .LT.-1,?1)  ) GO  TO  2300 
lF(BijFF(I-l)  .Eo,-1.0)  Go  TO  3n00 

IF((nuFF(I-l)  .lt.-1.o1).OR;(8i)FFII-1).GT.2o.()))  60  TO  2300 
IF (SENSW (2) ) GO  TO  3ooo 
2200  continue 

GO  TO  2100 

2300  WRITE(nPRNt*D1c0I  (I,BUFF(n»Isl,40) 

9100  F0RMAT(4(2X,I^,2A,E12.5)) 

2400  WRiTE(nTYP, 922(1) 

9220  format (54H  INPUT  ID  «=  TO  BUFFER  TO  CHANGE  AND  NEW  VALUE , I 2El2. 5 
1 /,16H  ID  .Ge.41  STOPS  ,/) 

REAr)(NINP,92lu)  lD,flijEF(IO) 

9210  F0RmaT( I2,E12.5) 

IF(ID.GE.41)  go  To  2450 
GO  TO  ^‘^00 

2450  9ACKSPACE  nOISK 

CALL  QmOND(19,nDISK,BUFfIi  ) ,B||FF(44)  ) 

GO  TO  2100 
3000  CONTINiJE 

CALL  EXIT  ; 

END 


B.2.2 


Control  Program 


The  control  program  transfers  the  axle  displacement  records 
from  a disc  storage  file  to  the  digital  to  analog  converters  as  shown  in 
Figure  9(b).  The  records  are  synchronized  by  a crystal  controlled  real  time 
clock  prior  to  the  transmission.  Figure  B-2  shows  a flow  chart  of  the 
program. 


B.2.2. 1 Listing 

An  EAI  FORTRAN  listing  of  the  program  appears  in  the 
pages  following  the  flow  chart. 
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Start 


Initialize  Pointers 
Rewind  Disc  File 


Figure  B-2.  Control  Program  Flow  Chart 


C TRUCK  SHAKER  PROUHAM 

C INPUT  data  file  must  be  POSITtonEo  to  logical  unit  21  TO  RUN 
DIMENSION  XBUFFdoO)  tZEROjlO) 
logical  SENSW 

data  ZERO/0<C(0«0t0«O*0.O«0.0t0*0>0*0«0«0t5*0«0t0/ 
call  QSHYIN(itRRt580) 

N0iSK=l7 

c equals  octal  21 

NTYPal 

NIriPsZ 

WRITE  (NTYP, 9001 ) 

9001  FORMATdSH  PROGRAM  SHAKER*/. 2iH  LU  «21  FOR  DATA  FILE  ) 
REWImO  nDISK 

c SET  interval  Timer  to  30  seconds 
TlMEaSoOOO.O 
CALL  QwTImR(TImE.IERR) 

C ONE  sector  = 8a  is  B.^.T  WORDs«OR  ll  time  steps, but  lO/RECOPD  USEO 
c Input  one  sector  to  xbuff 

CALL  0M0ND(ir*N0lSK,XBUFF(n»xnUFF(44)) 

C initialize  pointers 
iPNTsl 
INC=4 
NSTEP=10 

C set  0/a  lines  TO  ZERO 

CALL  QwHOAR(Zt.RO,0,T»lE9H) 

CALL  OSTOA 

C DISABLE  ERROR  MESSAGES  AND  WAIT  Fqr  SW  1 TO! GO  ON, TRUE 
CALL  OltlET 

C call  LOW  SPEEO  wait 
PO  CALL  OSULY(l'OO) 

C start  transmission  when  SENSW  1 TRi.lEtON 
IF<SENSW(1) ) Go  TO  75 
Go  To  50 

75  CALL  QwTIMR(TIHE*IERR» 

C SYNC  TO  program  ANO  CLOCK 
100  DO  200  I=1,NSTEP 
C TEST  FOR  EOF  ■ 

IF(XBUFF(IPNT) .EQ,-1,0)  Go  TO  300 
C LOAD  0/A  duffer  WItH  DATA 

call  OWBDAR (XBUFF (IPNT*1) ,I ,3,IERR) 
c wait  for  step 

XTIME=XBUFF (IPNT) *1000.0 
150  call  ORTIMR(RTImE,IERR» 

IF ( TIME-RT I ME-XT I ME) 150,140,160 
c after  clock  sync  transmit  Data  out  d/a  lines 
160  CALL  QSTOA 

IPNT--IpNT*INC 

C TEST  SENSW  3 FOR  END.  PROGRAM 
IF(SENSW(2))  go  To  5000 

200  continue 

C GET  NEXT  RECORD  FROM  DISK 

210  call  QM0ND(17,ndI5K,XBUFF(1),xrUFF(44)) 

IPNTal 
GO  To  100 

c end  of  data, rewind  disk  and  Restart 
300  rewind  NDISK 
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CALL  QWTIMRaiMEflERR) 

IREC=0 
GO  To  2l0 

c program  end  on  sensw  I link  Back  to  monitor 

GOOO  WRITE(NTYP.'?  10) 

90i0  F0RMAT(12H  END  PROGRAM  ) 

CAll  exit 
END 
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B.  3 Program  Generation  Details 


How  to  create  a core  image  version  of  a program  is  explained  in  this 
section.  Several  details  of  the  operation  are  not  included  in  the  description 
because  they  are  readily  available  in  the  EAI  reference  manuals  (reference 
4,  5,  6 and  7).  The  core  image  is  generated  according  to  the  following  three 
steps: 

a.  Creation  of  source  program  file  with  BASIC  TEXT  EDITOR 
program 

b.  Compilation  of  source  program  using  the  FORTRAN  compiler 

c.  Generation  of  core  image  file  with  CORE  IMAGE  GENERATOR 
program. 

A diagram  of  the  process  is  shown  in  Figure  B-3. 

B.3.1  Source  File  Generation 


The  source  file  creation  step  requires  a number  of  operations 
by  the  user.  First,  the  source  input  is  punched  on  paper  tape  in  a form 
required  by  the  text  editor.  (See  reference  6.)  Then  the  disk  file  for  the 
source  storage  is  created  with  a monitor  command,  and  the  Basic  Text  Editor 
is  loaded  and  executed.  The  complete  procedure  is  outlined  below: 

#N,  SOURCE,  22,  0 

where 


#N 

SOURCE 

22 

0 


create  file  entry  command 
file  name 

logical  unit  nvimber  of  file  opened,  octal 
identifies  source  file  type. 
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Figure  B-3,  Program  Generation  Sequence 
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#L,  BTE,  21 


where 


#L  = load  file  command 

BTE  = file  name  of  basic  text  editor 

21  = logical  \init  number  for  loader. 


#R,  220,  422,  17777 


where 


#R 

220 

422 

17777 


set  hardware  register  command 

command  input  = console  and  list  output  = printer 
text  input  = paper  tape  and  text  output  = disc 
biaffer  memory  space  (see  reference  6). 


#G,  1000  ^ 


NNN, 

w 

M ^ 


input  NNN  lines  of  program 
write  btiffer  to  disc  file 
write  end  of  file  mark 
return  to  monitor. 


If  changes  in  the  source  code  from  the  paper  tape  version  are  required  they 
would  be  done  with  the  editor  before  the  write  (W^)  instruction.  (See 
reference  6.) 


B.3.2  Program  Compilation 

The  FORTRAN  compiler  inputs  code  from  the  previously 
created  source  file  and  creates  an  object  file  which  the  program  loader 
requires.  An  object  file  for  the  compiler  output  must  be  created  before 
exe  cution. 
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Example  #N,  OBJECT,  23,  o 
where 


#N 

OBJECT 

23 

1 


create  disc  file 
name  of  object  file 
logical  unit  for  compiler  output 
identifies  object  file  type. 


#P,  SOURCE,  22 


where 


#P  = open  disc  file 

SOURCE  = file  name  from  text  editor  output 
22  = logical  unit  number  for  source  input. 

#L,  FORTR,  21  (FORTR  = compiler  file  name) 
#R,  , 20702,  41162  ]^ 


where 

#R  = register  set  command 

20702  = non-standard  compiler  options 

41162  = source,  object,  list  logical  units.  (See 

reference  4) 

#G,  1000 

If  the  program  compiles  correctly,  the  object  file  OBJECT  is  ready  for  input 
by  the  core  image  generator. 
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B.3.3 


Core  Image  Generation 


The  core  image  generator  takes  the  user's  object  file  and 
links  it  with  the  required  system  subroutines.  However,  before  initiating 
the  process,  an  output  file  must  be  created  and  the  object  file  positioned: 


Example  #N,  CORIMG,  23,  2 


where 


#N  = file  creation  command 

CORIMG  = file  name  of  core  image 
23  = logical  unit  file  opened  on 

2 = identifier  for  Cl  file  type. 


#P,  OBJECT,  22 


where 


#P  = command  to  open  disc  file 

OBJECT  = name  of  object  input  file 
22  = logical  unit  number  of  file. 


#L,  CIG,  21 


where 


#L  = loader  command 

CIG  = file  name  of  Core  Image  Generator 
21  = logical  unit  niomber  for  loader. 
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#R,  , 101000,  23000 


where 


#R  = register  setting  command 

101000  = non-standard  options,  standard  execution  address 

23000  = logical  unit  output  = 23,  zone  zero  starts  at  000 


#G,  1000  exe  cute 

#F,  OBJECT,  22  force  load  object  file  to  output 

#L,  LINKN,  24  ^ library  load  (selective)  required  subroutines  from 

hybrid  linkage  library 

#L,  RTL,  24  library  load  (selective)  required  subroutines  from 

run  time  library 


#m;^ 


close  file  on  disc  and  return  to  monitor. 


The  program  can  now  be  loaded  from  disc  and  executed  with  a single  load 
command  as  described  in  3,  2. 
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B* Data  Storage  Formats 


Axle  displacement  records  are  obtained  by  the  system  from  paper  tapes 
punched  in  ASCII  standard  eight  level  tape  format.  Corresponding  to  each 
time  step,  there  are  four  elements  punched  on  the  tape,  i.e.,  the  time  value 
(seconds),  front  axle  displacement  (feet),  middle  axle  displacement  (feet), 
and  the  rear  axle  displacement  (feet).  These  data  are  in  FORTRAN  format 
F7.4,  3F8.5  (see  2.2.2. 3),  therefore,  typical  data  look  like: 

Example  01.2345,  Time 

+0.12345,  Front  axle 
-0.98765,  Middle  axle 
+0.34567,  Rear  axle 


Counting  the  characters  in  the  above  show  that  each  time  step  requires 

31  frames  of  paper  tape.  These  values  are  stored  in  memory  and  on  disc  as 

32  bit  real  numbers,  therefore,  data  for  each  time  step  require  eight  16-bit 
words  of  storage  (specific  details  on  storage  can  be  found  in  EAI  publications, 
i.e.,  references  4,  5,  6 and  7). 


When  these  data  are  transferred  to  a disc  storage-l*  file  they  are  for- 
matted such  that  each  sector  (largest  block  of  data  that  can  be  transferred 
in  one  operation)  is  88  16 -bit  words.  This  allows  up  to  11  time  steps  to  be 
packed  into  each  sector.  The  program  actually  stores  10  time  steps  or  40 
values,  with  four  zeros  at  the  end  to  fill  up  the  sector. 


Example  1 
2 


00.0000  + 0.12345  - 0.98765  + 0.34567 
00.0039  + 0.12201  - 0.98321  + 0.34921 


10  00.0451  + 0.10310  - 0.93415  + 0.39216 

11  00.0000  00.00000  00.00000  00.00000 

EOR 


One 
Se  cto  r 


*Input/output  with  disc  storage  is  carried  out  by  calls  to  the  monitor  by  a 
system  routine  QMOND,  details  of  which  are  available  in  reference  4. 
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The  full  data  file  then  consists  of  a chain  of  blocks,  10  steps  each,  with  the 
final  block  containing  -1.0  fill  up  elements  to  the  EOR. 


Example 


EOR 

18.0000  + 1.00000  - 0.01023  + 0.97510 


0.0  0.0  0.0  0.0 

EOR 

18.0040  + 0.95610  - 0.020501  + 0.96321 


-1.0000  - 1.00000  - 1.00000  - 1. 00000  I Marks  end 

-1.0000  - 1.00000  - 1.00000  - 1.00000  f of  file 

EOR 

EOF  (End  of  Disc  File) 


Employing  -1  indicators  allows  file  lengths  to  be  changed  without  altering  the 
programs,  as  long  as  the  end  of  the  file  is  marked  by  a series  of  -1.0. 
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